Abstract
Field performance of southern highbush blueberry (Vaccinium corymbosum L. interspecific hybrids) cultivars Emerald, Jewel, and Primadonna derived from softwood cuttings (SW) and tissue culture (TC) was evaluated in Citra and Haines City, FL, in 2010–12. Both fields were planted in Apr. 2010 on sandy soil amended with pine bark. Plant height and width were recorded at both locations, from which plant canopy volume was calculated. Additionally, whole plants were harvested at planting and after the first growing season, after the first fruit harvest, and after the second growing season. Average plant height and width, number of major canes, and total shoot number were determined at each sampling date. Dry weights for roots, crowns, canes, shoots, and leaves were obtained. Although propagation method affected plant canopy volume during the first season, no effects were observed by the end of the second growing season. At planting and after the first and second growing seasons, TC plants of the three cultivars had more major canes. Total shoot number per plant was greater for TC ‘Jewel’ at all dates but ‘Emerald’ TC plants had more shoots only at planting and after the first growing and harvest seasons. Tissue culture resulted in increased plant dry weights of ‘Jewel’ and ‘Emerald’ after the first and second growing seasons. There were no significant differences in total number of shoots or plant dry weight between TC and SW-derived ‘Primadonna’ plants at any point during the study.
Blueberries have traditionally been propagated by softwood, semihardwood, or hardwood cuttings. Although plants can be produced at low cost, cuttings from some cultivars have low or very low rooting percentages, particularly when propagated from older, mature plants (Lyrene, 1981; Miller et al., 2004). Therefore, it may take several years to propagate and commercialize newly released cultivars so that significant acreage can be planted (Lyrene, 1980). In contrast, micropropagation has the potential to produce large numbers of plants more quickly than by rooted cuttings, and tissue culture plants can be produced throughout the year (Isutsa and Pritts, 1994; Miller et al., 2004).
Although changes in the growth habit of micropropagated Vaccinium and other small fruit crops have been reported, field performance evaluations of TC-derived plants are limited and results are variable and contradictory (Debnath, 2007). During the first 2 years of measurement, TC plants of half-highbush ‘Northblue’ (Vaccinium corymbosum L. × Vaccinium angustifolium Ait.) had significantly more lateral branching than SW plants (Grout et al., 1986; Read et al., 1988, 1989), whereas Albert et al. (2009) reported that ‘Northblue’ SW plants were taller and produced significantly more shoots than TC plants during the first two growing seasons. ‘Herbert’ TC highbush blueberry plants grew more vigorously and more uniformly than SW-derived plants during the first 3 years in the field (Litwiñczuk et al., 2005). However, in another study, propagation method did not have a significant effect in the number of shoots or plant size of 1-year-old highbush ‘Bluecrop’ (Smolarz and Chlebowska, 1997).
Because previous research had been done in temperate climates with much shorter growing seasons than those of Florida, our objective was to compare various vegetative growth traits of southern highbush blueberry cultivars obtained from TC and SW under Florida field conditions. Three commonly planted cultivars were assessed in production regions of Florida that differ primarily by the amount of winter chill accumulation.
Materials and Methods
Plant material.
‘Emerald’, ‘Jewel’, and ‘Primadonna’ were chosen because they are commonly used in commercial operations in the two test locations. Plants derived from TC were purchased from a nursery in Oregon and plants propagated from SW were purchased from two nurseries in Florida. Plants varied from 11 to 50 cm tall on arrival across all cultivars and propagation types. The exact age of the plants was unknown, but based on the typical timing for rooting softwood cuttings in Florida and the apparent age and number of growth flushes on the tissue culture plants, the age of the plants was less than 1 year.
Locations and experimental design.
The study was conducted at two locations within the major southern highbush blueberry production regions in Florida. The two sites have different average chill hour accumulations. The north site, located at University of Florida Plant Science Research and Education Unit in Citra, FL, receives between 420 to 540 chill hours, whereas the southern location in Haines City (HC), FL, ranges from 110 to 210 h/year (hours below 7 °C until 10 Feb. in 75% of winters; Ferguson et al., 2008). ‘Jewel’, ‘Emerald’, and ‘Primadonna’ have estimated chill hour requirements of less than 200, 250, and 400, respectively (Fall Creek Farm and Nursery, 2014). Treatments were arranged in a three × two factorial (three cultivars and two propagation methods) with five replications arranged in a completely randomized block design with 20-plant plots. The field at Citra was planted on 15 and 16 Apr. 2010 on beds of Arredondo sand (loamy, siliceous well-drained soils), whereas the field at HC was planted on 23 Apr. 2010 in Candler Sand (sand or fine sand, excessively drained soils). Both fields had been amended with incorporated pine bark before planting. At both locations, distance between rows was 3 m and plant spacing within the row was 0.75 m, and beds were covered with black polypropylene woven fabric for weed control.
Management practices.
The field at Citra was irrigated with a single drip tube, which delivered ≈374.2 L·min−1·ha−1. Irrigation was run in three daily cycles that totaled 45 min to 1 h·d–1. A complete liquid fertilizer (10N–2P–5K with micronutrients) was injected into the irrigation water. The first year, nitrogen (N) rate was 134 kg·ha−1 from planting until the end of Sept. 2010. During the second year, 100.8 kg·ha−1 was applied from the end of January until the middle of September. Overhead irrigation was used for freeze protection as necessary, depending on the year. Standard commercial practices were followed for weed, pest, and disease control (Williamson et al., 2013). Plants were not pruned the first year but, following standard commercial practices, plants were pruned the second year of the study. In Citra, plants were pruned 1 week after the first harvest season was completed. Plants selected for the June 2011 destructive sampling experiment were flagged and not pruned.
In HC, plants were irrigated with two rows of drip tube, delivering 561.2 L·min−1·ha−1. Irrigation was applied in five to eight cycles that totaled between 75 to 120 min·d–1. A complete 6N–6P–3K (with micronutrients) liquid fertilizer was injected once a day in spring and fall and three times per day during summer growth flushes for a total of 224 kg N/ha/year. Plants were hand-pruned at the end of the commercial harvest season (before destructive sampling on June 2011), only removing the old flower spikes while the new vegetative flush was allowed to grow.
Vegetative field measurements.


Additionally, number of major canes and total number of shoots per plant were recorded from the plants selected to evaluate reproductive growth. Major canes were defined as stems arising from the first 12 cm above the soil line with a minimum diameter of 8 mm. Shoots that were shorter than 5 cm and had a diameter of less than 2 mm were not counted. Data from HC were taken on 10 Dec. 2010 and at Citra on 12 and 14 Jan. 2011. Number of major canes and total shoot number of two plants per plot were recorded. For all measurements, data from the subsample plants were averaged and used for statistical analysis.
Destructive sampling.
Whole plants were harvested at four biologically significant points during the study: at the beginning of the experiment (Apr. 2010), after the first season of vegetative growth (Nov./Dec. 2010), after the first harvest (June 2010), and after the second season of growth (Nov. 2011). At planting, 10 representative plants of each cultivar were selected and stems arising 12 cm from the soil line with a minimum diameter of 8 mm were recorded as major canes. Lateral branches of each cane were also counted and plants were divided into roots, shoots, and leaves. Roots were washed free of soil and all plant parts were placed in paper bags and dried separately at 60 °C until constant weights were achieved.
After the first growing season (22 Nov. in HC and 3 Dec. in Citra), one plant from each plot was excavated. Roots were immediately separated from the tops about at the soil line and bagged separately. Roots were placed in cold storage at 4 °C until they were washed before drying. Major canes and total number of shoots with a minimum of 2 mm diameter or 5 cm in length were counted. All plant organs were separated, dried, and weighed as previously described except that crowns were dried and weighed separately.
Destructive harvests of five plants per treatment were done after the first fruit harvest on 8 June 2011 in Citra and on 14 June 2011 in HC and after the second growing season on 9 Nov. and 28 Nov. 2011 in Citra and HC, respectively. The plant organs were separated, dried, and weighed as described previously.
Statistical analysis.
Analysis of variance (ANOVA) was performed using PROC GLM with cultivar and propagation type as main effects (SAS Version 9.2; SAS Institute Inc., Cary, NC). Treatment means were separated using Tukey’s honestly significant difference test at α = 0.05 level of significance. PROC GLIMMIX was used to analyze number of branches and total shoot tips in the second destructive sampling because of unbalanced data. When data transformation was necessary, treatment means were separated based on the transformed data, but means of non-transformed data are presented. Whenever interactions involving location (location by cultivar by propagation and location by propagation) were not significant, means of both locations are reported.
In cases in which the overall effect of propagation was not significant but interactions between main effects were significant (cultivar by propagation, location by propagation, or location by cultivar by propagation), the slice statement was used. This statement produces a separate ANOVA with a P value for propagation method for each cultivar without having to analyze the data for each cultivar separately, which would reduce the df of the analysis (Littell et al., 2002).
Data that did not follow a normal distribution were transformed using the logarithmic function. This included plant volume data from June, July, and Nov. 2010; plant volume; number of canes per plant; and root, shoot, leaf and plant dry weights data before planting (Apr. 2010). From data recorded in June 2010, leaf, shoot, and cane dry weights and from Nov. 2011 leaf and root dry weight data were transformed using the logarithmic function to approximate the normal distribution. Data for leaf and shoot dry weights produced during the first year in the field and data for number of shoots produced between the end of the first and second years were transformed using the logarithmic function as well.
Results
Plant height and canopy volume.
During the first season, propagation effect on plant height was the same across locations. However, the cultivar by propagation interaction for the same period was significant as a result of increased plant height for TC ‘Jewel’ and reduced plant height for TC ‘Primadonna’ compared with SW plants (data not shown). Throughout the first season, plants in Citra were larger than in HC (data not shown). Significant interactions of propagation type and cultivar for plant canopy volume were observed throughout the season (Supplemental Table 1 shows main effects and hierarchy of interactions). Canopy volume of TC ‘Jewel’ plants was significantly larger than for SW plants, whereas the opposite was observed for ‘Primadonna’. ‘Emerald’ TC plants had significantly larger volume in June, July, and November with no significant differences between propagation types during the rest of the first growing season (Table 1).
Effect of cultivar and propagation method on plant volume (dm3) in 2010.


For much of the second season, plant height was unaffected by propagation method (data not shown). Significant location by cultivar by propagation interaction (P values shown in Supplemental Table 2) from June to Aug. 2011 is evidenced by larger plant canopy volume of ‘Emerald’ and ‘Jewel’ TC plants and reduced plant volume of ‘Primadonna’ TC plants compared with SW plants in HC, whereas in Citra, TC ‘Jewel’ had larger plant volume only in June without any other effect of propagation method between June and August (Table 2). In Citra, propagation method did not have a significant effect on plant volume of ‘Primadonna’ at any point during the season. Averaged across locations, ‘Emerald’ and ‘Jewel’ plant volumes were greater for TC than SW plants in September and October (Table 3). Conversely, there was no significant effect of propagation type on plant volume of ‘Primadonna’ plants for the same months. In Nov. 2011, propagation effect on plant canopy volume, across cultivars and locations, was not significant (P = 0.0644).
Effect of cultivar and propagation method on plant volume (dm3) in 2011 at two locations.


Effect of cultivar and propagation method on plant volume in Sept., Oct., and Nov. 2011.


Number of major canes and total shoots/plant.
TC plants of all cultivars had more major canes at planting (Table 4) and at the end of the first and second growing seasons (Table 5). After the first harvest (June 2011), significant cultivar by propagation interaction is evidenced by increased number of canes for TC ‘Emerald’ and ‘Jewel’ plants compared with plants from SW cuttings with no significant effect of propagation method on ‘Primadonna’. TC ‘Emerald’ and ‘Jewel’ plants had more total shoots per plant at planting and on Dec. 2010 (Fig. 1) and June 2011 (Table 5). Additionally, TC ‘Jewel’ plants had more shoots per plant at the end of the second growing season. There was no significant difference in total number of shoots per plant between TC and SW ‘Primadonna’ plants at any point during the study.
Effect of cultivar and propagation method on initial plant size and dry weights before planting (Apr. 2010).


Effect of cultivar and propagation method on the number of major canes and total shoots/plant at three field sampling dates performed during the study.



Effect of propagation method on ‘Jewel’ plant architecture. Haines City, Dec. 2010. Similar results were observed with ‘Emerald’ plants from the two propagation methods. JSW = ‘Jewel’; softwood cuttings; JTC = ‘Jewel’, tissue culture.
Citation: HortScience horts 49, 5; 10.21273/HORTSCI.49.5.556

Effect of propagation method on ‘Jewel’ plant architecture. Haines City, Dec. 2010. Similar results were observed with ‘Emerald’ plants from the two propagation methods. JSW = ‘Jewel’; softwood cuttings; JTC = ‘Jewel’, tissue culture.
Citation: HortScience horts 49, 5; 10.21273/HORTSCI.49.5.556
Effect of propagation method on ‘Jewel’ plant architecture. Haines City, Dec. 2010. Similar results were observed with ‘Emerald’ plants from the two propagation methods. JSW = ‘Jewel’; softwood cuttings; JTC = ‘Jewel’, tissue culture.
Citation: HortScience horts 49, 5; 10.21273/HORTSCI.49.5.556
Destructive sampling.
Plant height and root and shoot dry weights were greater for TC plants than for SW plants for each cultivar at planting (Table 4). Initially, average width of ‘Emerald’ and ‘Jewel’ was greater for TC plants than for SW plants, whereas the opposite was observed for ‘Primadonna’. Leaf dry weight of TC plants was greater in ‘Jewel’ but reduced in ‘Primadonna’ plants compared with SW plants, whereas propagation type had no effect on ‘Emerald’ leaf dry weight. ‘Emerald’ and ‘Jewel’ TC had greater plant dry weight (data not shown) and plant canopy volume before planting with no effect of propagation type on ‘Primadonna’ plant dry weight or canopy volume (Table 4).
At the end of the first season (Nov./Dec. 2010), there was no significant effect of propagation type on plant height, plant canopy volume (data not shown), or root dry weight for any of the cultivars at either of the two locations (Table 6). However, significant location by propagation interactions were observed for leaf and cane dry weights. In Citra, TC plants had greater leaf dry weight compared with SW plants (SW = 110.68 g vs. TC = 145.39 g, P = 0.0413) and greater cane dry weight than SW plants (SW = 132.23 g vs. TC = 184.33 g, P = 0.0319) without a significant difference among propagation types for these variables in HC. ‘Emerald’ and ‘Jewel’ TC had more canes, lateral branches, and total number of shoots per plant (Table 5) as well as greater plant width and volume and shoot (data not shown) and plant dry weight compared with SW plants (Table 6). However, propagation type had no significant effect on total number of shoots, volume, or plant dry weight of ‘Primadonna’ at the end of the first growing season.
Effect of cultivar and propagation method on plant dry weights in Nov./Dec. 2010, June 2011, and Nov. 2011.z


In June 2011, ‘Emerald’ and ‘Jewel’ TC had significantly greater plant height and canopy volume (data not shown) and plant dry weights compared with SW plants without a significant effect on canopy size or dry weight of ‘Primadonna’ after ≈14 months of field growth (Table 6). Additionally, root dry weight was greater for ‘Emerald’ and ‘Jewel’ TC, whereas for ‘Primadonna’, greater root dry weight was recorded on SW plants. There was a significant location by propagation interaction for root:shoot ratio: TC plants had smaller root:shoot ratio in Citra (SW = 0.72 vs. TC = 0.45, P ≤ 0.0001) with no significant difference between propagation types in HC (SW = 0.33 vs. TC = 0.25, P = 0.0645). Although TC significantly increased leaf dry weight in both locations, treatment differences were larger in Citra (SW: 59.52 g vs. TC: 129.7 g; P ≤ 0.0001) than in HC (SW: 55.6 g vs. TC: 69.6 g, P = 0.0011).
After the second year in the field, plant height was not significantly affected by propagation method (data not shown). Conversely, TC plants had larger mean diameter (SW = 92.22 cm vs. TC = 99.12 cm, P = 0.0384) and greater plant volume (SW = 569.5 dm3 vs. TC = 664.2 dm3, P = 0.0399). Cultivar by propagation interaction is shown by increased cane, total shoot (data not shown), and plant dry weights of ‘Emerald’ and ‘Jewel’ TC plants but not of ‘Primadonna’ (Table 6). ‘Emerald’ and ‘Jewel’ root dry weights were unaffected by propagation type but TC ‘Primadonna’ plants had reduced root dry weight compared with SW plants. A significant location by cultivar by propagation type interaction was observed for leaf dry weight; there was no effect of propagation type on any cultivar in Citra, whereas in HC, ‘Emerald’ TC plants had higher leaf dry weight (ESW = 80.18 g vs. ETC = 134.98, P = 0.0033) and ‘Jewel’ TC plants had reduced leaf dry weight compared with SW plants (JSW = 73.42 g vs. JTC = 49.44 g, P = 0.0381) and there was no effect of propagation type on ‘Primadonna’.
Increase in cane and total shoot number and dry weights during the first 2 years of field growth.
To determine the effect of initial plant size on final vegetative measurements, ANOVA was performed on the difference in measurements taken after one season of field growth (Nov. to Dec. 2010) and measurements taken at planting (Apr. 2010) and on the difference between the end of the first and the end of the second seasons of field growth.
The effect of propagation on plant canopy volume (data not shown) and total number of shoots produced during the first year in the field were not significant (Table 7). Conversely, from planting until the end of the first growing season, TC plants produced more new canes than plants derived from SW (SW = 4.80 vs. TC = 5.97, P = 0.0241). Leaf, shoot, and plant dry weights produced by TC ‘Emerald’ and ‘Jewel’ plants during the same period were significantly greater than values recorded for SW plants of the same cultivars, whereas propagation effect on the same variables was not significant for ‘Primadonna’ (Table 7).
Effect of cultivar and propagation method on total number of shoots and leaf, shoot and plant dry weights produced during the first season of field growth.z


Plant canopy volume produced between the end of the first and the second growing seasons was not significantly different between propagation types (P = 0.0937). However, TC plants grew more canes (SW = 1.8 vs. TC = 3.9, P = 0.0022) and total shoots per plant (SW = 121.7 vs. TC = 150.9, P = 0.0114). Additionally, TC plants produced more total shoot (SW = 534.66 g vs. TC = 700.92 g, P = 0.0003) and plant dry weight compared with SW plants (SW = 1079.12 g vs. TC = 1277.78 g, P = 0.0102).
Discussion
Differences in plant volume resulting from propagation method were larger at the beginning of the study but tended to decrease during the second year. At the end of the second growing season (Nov. 2011), plant volume was not significantly different between propagation types for any of the cultivars in the study (Table 2). Greater plant canopy volume of TC ‘Emerald’ and ‘Jewel’ observed in HC at the beginning of the second season could be attributed to greater plant width, because there was no significant effect on plant height for most of the second season. Similar results were obtained with lowbush blueberry, where TC plants had greater diameter at ground level than SW plants (Jamieson and Nickerson, 2003).
In addition to soil type and climate variations between our test locations, differences in management practices could have contributed to the differences in plant growth observed between the two locations. The reduced plant dry weights observed in HC during the second season could be attributed to insufficient chilling (Darnell and Davies, 1990) or the result of differences in climate and cultural practices between the two locations. Smaller plant dry weights recorded at the end of the second season for HC compared with Citra were mostly the result of differences in leaf and root dry weights. Smaller leaf dry weights in HC may be attributed to insufficient chilling but could also be the result of weather conditions favorable to leaf diseases (causing earlier and greater leaf drop), which were observed to be more prevalent in HC. Reduced root dry weight in HC compared with Citra may have resulted from differences in leaf number and health or differences in management practices such as irrigation and fertilization.
In June and July of the second year, plants at Citra had reduced plant volume compared with plants in HC (averaged across all levels of the other factors), and there was a significant location by cultivar by propagation interaction. One possible explanation for these differences is that only two of the data plants per plot were harvested for yield data in Citra and the rest of the plants in each plot were not picked almost until the end of the season, whereas in HC, the field was part of a commercial operation and all plants were harvested regularly. Moreover, in Citra, fruit yields were much higher than in HC, adding further to the stress brought on by young plants carrying a heavy cropload, which often results in reduced vegetative growth. Increased plant dry weight of young highbush blueberry cultivars when flower buds were removed has been reported (Strik and Buller, 2005; Williamson and NeSmith, 2007). The hypothesis that smaller plant volume in Citra compared with HC, in June and July 2011, was a result of the high cropload is also supported by the data from the destructive sampling after harvest (June 2011), where leaf, shoot, and whole plant dry weights were greater in HC, although plants at that location had been pruned before the experiment. Conversely, leaf, total shoot, and plant dry weights were larger in Citra than in HC at the first destructive sampling (end of 2010) and leaf, root, and plant dry weights were greater at the following sampling date (end of 2011).
Pruning differences between locations may have also contributed to the differences in plant volume between the two locations. Plants in HC were pruned earlier (4 weeks before measurements) than in Citra, where plants were pruned 2 weeks before starting measurements in June 2011. Moreover, differences in pruning style may have also played a role in the location effect on canopy volume. In Citra, manual topping and pruning included removal of canes that did not grow upright and/or originated too close to the soil line. However, in HC, pruning consisted mostly of topping the plants with much lighter cane removal. This can also explain the interaction between treatments and location observed at this point: in Citra, where plants had been pruned just 3 weeks before recording height and width in June 2011, there were less differences in plant volume between propagation types but in HC, where plants had time to grow after light summer pruning and before measurements were taken, we observed larger plant volume of TC ‘Emerald’ and ‘Jewel’ plants at the beginning of the second season.
Tissue culture resulted in greater plant dry weight at the end of the second year in the field with ‘Emerald’ and ‘Jewel’, which agrees with results obtained in TC lingonberry and blackberry plants (Gustavsson and Stanys, 2000; Swartz et al., 1983). Effect of propagation method on number of laterals also varied among different blackberry cultivars (Swartz et al., 1983). Here, during the first year after planting, TC ‘Emerald’ plants produced more shoots than plants from SW, whereas there was no difference in the other two cultivars. The greater number of shoots per plant observed in TC ‘Jewel’ after 1 year may have been the result of differences in plant age or size, between propagation types, at the beginning of the study. Although many of the previous studies done to evaluate field performance of TC and SW plants did not include initial measurements, Grout et al. (1986) reported that TC ‘Northblue’ plants increased 46% in size (centimeters per plant) compared with 24% increase for plants from cuttings under greenhouse growth conditions. At planting, for the present experiment, there were significant differences in plant height, width and canopy volume, and whole plant and certain organ dry weights between propagation methods. Because independence of the covariates (the initial values) with respect to the treatments (the plant itself represents the treatment) was not met (Cochran, 1957; Elashoff, 1969), we could not use analysis of covariance to determine the effect of initial plant size on final measurements. Therefore, ANOVA was performed on the difference in measurements taken after one season of field growth (Nov. to Dec. 2010) and measurements taken before planting (Apr. 2010) and on the difference between the end of the first season and the end of the second season of field growth. With this analysis, we observed results similar to previous work with significantly greater increase in whole plant dry weight for ‘Emerald’ and ‘Jewel’ TC plants compared with SW regardless of the initial plant dry weight.
We observed that 8- and 14-month-old ‘Emerald’ and ‘Jewel’ TC plants had greater total shoot number than SW plants and after two growing seasons, tissue culture plants of all the cultivars tested had more canes than plants propagated from rooted cuttings. Previous work with ‘Northblue’ blueberry resulted in TC having more basal branches (canes) and lateral branches or shoots than plants from rooted cuttings during their first 3 years of field growth (Litwiñczuk et al., 2005; Read et al., 1989). Whether the increased branching observed in blueberry plants derived from TC compared with SW lasts for 1, 2, or more years, it still poses a significant advantage during establishment years and can result in a greater number of flower buds per plant and greater early yields (Grout et al., 1986). The effect of propagation method on reproductive growth and yield and the relationship between shoot number and yield for the present trial is under evaluation and will be reported in a separate paper.
Despite differences in management practices between the two locations used in this study, the effect of propagation method on the three cultivars studied was, overall, similar in Citra than in HC. During the first season in the field, ‘Emerald’ and ‘Jewel’ TC plants tended to have greater plant canopy volume than SW plants, whereas the opposite was observed for ‘Primadonna’ (June to Nov. 2010). ‘Emerald’ and ‘Jewel’ TC plants had more major canes and plant dry weight after the first and second growing seasons with no significant effect of propagation type on plant dry weight or total shoot number for ‘Primadonna’ at any point during the study. Total shoots per plant were greater for TC ‘Emerald’ and ‘Jewel’ after one season in the field but only for ‘Jewel’ after two seasons compared with SW plants. The effect of propagation method on vegetative growth and field performance seems to be cultivar-dependent and more trials should be conducted in several locations with these and other southern highbush blueberry cultivars before general recommendations can be made.
Literature Cited
Albert, T., Starast, M., Karp, K., Kaldmae, H., Vool, E. & Paal, T. 2009 The influence of propagation method on growth of the half-highbush blueberry ‘Northblue’ Acta Hort. 812 141 146
Cochran, W.G. 1957 Analysis of covariance: Its nature and uses Biometrics 13 261 281
Darnell, R.L. & Davies, F.S. 1990 Chilling accumulation, budbreak, and fruit set of young rabbiteye blueberry plants HortScience 25 635 638
Debnath, S.C. 2007 Propagation of Vaccinium in vitro Intl. J. Fruit Sci. 6 47 71
Elashoff, J.D. 1969 Analysis of covariance: A delicate instrument Amer. Educ. Res. J. 6 383 401
Fall Creek Farm and Nursery 2014 Southern highbush. Low chill varieties. 2 Feb. 2014. <http://www.fallcreeknursery.com/commercial/variety/commercial_southern-highbush>
Ferguson, J.J., Andersen, P.C., Chaparro, J.X. & Williamson, J.G. 2008 Florida subtropical peaches: General concepts and cultivars recommended for grower trial. University of Florida Coop. Ext Serv., Gainesville. Publ. HS1125
Grout, J.M., Read, P.E. & Wildung, D.K. 1986 Influence of tissue culture and leaf-bud propagation on the growth habit of ‘Northblue’ blueberry J. Amer. Soc. Hort. Sci. 111 372 375
Gustavsson, B.A. & Stanys, V. 2000 Field performance of ‘Sanna’ lingonberry derived by micropropagation vs. stem cuttings HortScience 35 742 744
Jamieson, A.R. & Nickerson, N.L. 2003 Field performance of lowbush blueberry propagated by seed, stem cuttings and micropropagation Acta Hort. 626 431 436
Isutsa, D.K. & Pritts, M.P. 1994 Rapid propagation of blueberry plants using ex vitro rooting and controlled acclimatization of micropropagules HortScience 29 1124 1126
Littell, R.C., Stroup, W.W. & Freund, R.J. 2002 SAS for linear models. 4th Ed. SAS Institute. A two way factorial experiment. p. 74–80. 11 Nov. 2013. <http://books.google.com/books?id=JhxKC8ewcdAC&dq=slicing+proc+glm&source=gbs_navlinks_s>
Litwiñczuk, W., Szczerba, G. & Wrona, D. 2005 Field performance of highbush blueberries (Vaccinium corymbosum L.) cv. ‘Herbert’ propagated by cuttings and tissue culture Sci. Hort. 106 162 169
Lyrene, P.M. 1980 Micropropagation of rabbiteye blueberries HortScience 15 80 81
Lyrene, P.M. 1981 Juvenility and production of fast-rooting cuttings from blueberry shoot cultures J. Amer. Soc. Hort. Sci. 106 396 398
Miller, S.A., Rawnsley, E.K., George, J. & Patel, N. 2004 A comparison of blueberry propagation techniques used in New Zealand Acta Hort. 715 397 402
NASA 2010 National Aeronautics and Space Administration. 28 Jan. 2014. <http://www.grc.nasa.gov/WWW/BGH/volume.html>
Read, P.E., Hartley, C.A., Sandahl, J.G. & Wildung, D.K. 1988 Field performance of in vitro propagated blueberries Proc. Intl. Plant Prop. Soc. 37 450 452
Read, P.E., Wildung, D.K. & Hartley, C.A. 1989 Field performance of in vitro-propagated ‘Northblue’ blueberrries Acta Hort. 241 191 194
Smolarz, K. & Chlebowska, D. 1997 Growth, vigour and yielding of highbush blueberry cv. Bluecrop propagated from semi-woody cuttings and in vitro J. Fruit Ornment. Fruit Res. 5 53 60
Strik, B. & Buller, G. 2005 The Impact of early cropping on subsequent growth and yield of highbush blueberry in the establishment years at two planting densities is cultivar dependant HortScience 40 1998 2001
Swartz, H.J., Galletta, G.J. & Zimmerman, R.H. 1983 Field performance and phenotypic stability of tissue culture-propagated thornless blackberries J. Amer. Soc. Hort. Sci. 108 285 290
Williamson, J.G., Harmon, P.H., Liburd, O.E. & Dittmar, P. 2013 Florida blueberry integrated pest management guide. Univ. of Florida Coop. Ext. Serv. Publ. HS1156
Williamson, J.G. & NeSmith, D.S. 2007 Evaluation of flower bud removal treatments on growth of young blueberry plants HortScience 42 571 573
P values from analysis of variance for plant volume in 2010.


P values from analysis of variance for plant volume in 2011.

