Agronomic Performance of Di ﬀ erent Open-Pollinated Beetroot Genotypes Grown Under Organic Farming Conditions

: Due to the increase of the organic cultivation of beetroot and its growing importance as a functional food, the potential advantages of open-pollinated genotypes as an alternative for F1 hybrid genotypes has been investigated. In this regard, six ﬁeld experiments were carried out in 2017 and 2018 in three di ﬀ erent locations to investigate the agronomic potential of new and existing open-pollinated genotypes of beetroot and their performance under the speciﬁc conditions of organic agriculture. Fifteen beetroot genotypes, including one F1 hybrid as a commercial control and one breeding line, were compared regarding their total and marketable yield, individual beet weight, diameter of beet, and leaves-growth-base width. Furthermore, ﬁve randomly selected beetroots meeting common marketable criteria were evaluated optically with regard to skin smoothness, corky surface, root tail, scab incidence, and uniformity. Results of this study indicated a signiﬁcant impact of genotype on eight of the assessed traits. The cylindrical-shaped genotype, Carillon RZ, demonstrated signiﬁcantly higher total and marketable yields, with 53.28 ± 3.34 t ha − 1 and 44.96 ± 3.50 t ha − 1 , respectively, compared to the yellow-colored genotype, Burpees Golden, which obtained the lowest total yield, 36.06 ± 3.38 t ha − 1 , and marketable yield, 27.92 ± 3.55 t ha − 1 . Moreover, the comparison of the open-pollinated genotypes with the F1 hybrid, Monty RZ F1, revealed that except for the traits yield, scab, and uniformity, the open-pollinated genotypes indicated desirable competitive outcomes and thus o ﬀ er suitable alternatives for organic cropping systems. Overall, the observed genetic variability can be beneﬁcial for breeding and food product development.


Introduction
With an average of 15% annual increase in cultivation area and a total of approximately 1900 ha (2018), beetroot is a commercially important vegetable in Germany [1]. The high demand for beetroot is also reflected in the quantities purchased per household in Germany, which amounted in 2018 to two kilograms [2]. Approximately 40% of the beetroot produced in Germany grows on organic farmlands, which makes this vegetable a typical organic product [3]. Beetroot is often offered fresh, as well as pre-cooked and vacuum-packed [1].
Beetroot (Beta vulgaris L. ssp. vulgaris) belongs to the Chenopodiaceae family. It is considered a classic winter vegetable and its pollination is naturally carried out by wind [4]. In the first year, the beet is formed by the thickening of the shoot axis section below the cotyledons. Depending on the genotype, the beet has a round or cylindrical shape. For seed production, beets are replanted in the second year. Recently, with the classification of beetroot as a functional food [5], the importance of this vegetable has further increased. The term functional food generally covers those foods showing health-promoting and disease-preventing benefits due to their nutritive value [6]. Moreover, relatively easy cultivation and good storability without the need for costly storage equipment has increased farmers' as well as the food industry's interest in beetroot production [7]. In addition, information on the morphology of the root crops, especially the shape and size of the storage organ, due to its impact on processing, consumer choice, and transportation, is necessary [8] and consequently it can lead to quality improvement of the crop products [9].
The demand for organic products is constantly increasing [10]. According to different studies which compared the quality of organic and conventional food products, a higher amount of health-benefiting compounds and lower amounts of harmful components were reported in organically produced foods [11][12][13][14]. For instance, the study of Bavec et al. [10] depicted that the amount of total sugar content, total phenol content, and total antioxidant activity was higher in organically grown beetroot (cv. Rote Kugel) compared to those which were cultivated under conventional conditions. Moreover, in terms of sensory quality, red beetroot grown under organic conditions had a better taste and smell [15]. In order to meet the requirements and comply with the values of organic agriculture, cultivars need to be adapted to the agronomic conditions of organic farming [16]. One of the main criteria for organic crops is their ability to reproduce. Hence, open-pollinating genotypes are of major interest as the farmers can multiply the seeds themselves instead of depending on external inputs [17,18]. The other important principle of organic plant breeding is that the characteristics of the selected plant can be passed on homogeneously to the next generations [19]. Despite the fact that the F1 hybrids are known to be homogeneous genotypes, the resulting offspring are typically heterogeneous [20]. F1 hybrids are referred to those seeds or plants of the first filial generation which result from a cross-mating of different parental types [21]. Furthermore, the F1 hybrid genotypes of beetroot show a narrower genetic base in comparison with open-pollinated cultivars [22]. Consequently, the open-pollinating cultivars represent more advantages for the use in organic farming compared to F1 hybrids.
The objectives of the present study were to assess the agronomic performance of 14 new and existing open-pollinated beetroot genotypes, including one breeding line, and one F1 hybrid under the specific conditions of organic farming. The genotypes were evaluated for their genetic yield potential and morphological traits and their suitability and stability in different years and locations.

Field Experiments and Plant Materials
Field trials were carried out in 2017 and 2018, at three different locations. Two of the locations were the on-farm breeding locations coordinated by Kultursaat e.V.: De Beersche Hoeve (Oostelbeers, Netherlands), and Horticulture station Heinze (Bingenheim, Hessia, Germany). The third location was the research station for organic farming Kleinhohenheim (University of Hohenheim, Stuttgart, Baden-Wuerttemberg, Germany).
The on-farm breeding location De Beersche Hoeve is located at a latitude of 51 • 28 47 N, longitude of 5 • 15 35 E, and 17 m above sea level. In 2017, during the growth period (July-October), the mean temperature was 16 • C and the mean precipitation reached 47 mm. In 2018, the mean temperature and precipitation during the growth period (June to September) were 16.5 • C and 27.7 mm, respectively. The soil texture is fine sand with a low clay content.
The Horticulture station Heinze (50 • 22 28 N, 8 • 53 45 E) is located at an altitude of 129 m above sea level. In 2017, during the growth period (June to September), the mean temperature was 16.2 • C and the mean rainfall was 53.5 mm. In 2018, the mean temperature and the precipitation during June to September were 20.7 • C and 24.1 mm, respectively. The soil texture is clayey loam.
The research station Kleinhohenheim is located at a latitude of 48 • 44 14 N, longitude 9 • 12 01 E, and 430 m above the sea level. In 2017, during the months July to October, the mean precipitation was 77.42 mm and the mean temperature 17.6 • C. In 2018, the mean rainfall reached to 38.2 mm and the mean temperature was 19.0 • C between June and September. The soil texture is considered as loess to loamy clay.
Total monthly precipitation as well as monthly mean temperature of all three research stations during the experimental period in 2017 and 2018 are shown in Figure 1A-C.
Agronomy 2020, 10, x FOR PEER REVIEW 3 of 16 The Horticulture station Heinze (50°22′28 N, 8°53′45 E) is located at an altitude of 129 m above sea level. In 2017, during the growth period (June to September), the mean temperature was 16.2 °C and the mean rainfall was 53.5 mm. In 2018, the mean temperature and the precipitation during June to September were 20.7 °C and 24.1 mm, respectively. The soil texture is clayey loam.
The research station Kleinhohenheim is located at a latitude of 48°44′14 N, longitude 9°12′01 E, and 430 m above the sea level. In 2017, during the months July to October, the mean precipitation was 77.42 mm and the mean temperature 17.6 °C. In 2018, the mean rainfall reached to 38.2 mm and the mean temperature was 19.0 °C between June and September. The soil texture is considered as loess to loamy clay.
Total monthly precipitation as well as monthly mean temperature of all three research stations during the experimental period in 2017 and 2018 are shown in Figure 1A-C. All field trials were conducted under organic farming conditions in accordance with the European Commission's main regulations for organic farming. Briefly, no synthetic fertilizer was applied to the soil. Pests and disease control were done using biological controls. Frequent mechanical or manual weeding was conducted to suppress the weeds. The genotypes were allocated to the six field trials according to a non-resolvable row column design or a randomized complete block design (Table A1). In 2017, on each on-farm breeding location 30, and in Kleinhohenheim, 40 All field trials were conducted under organic farming conditions in accordance with the European Commission's main regulations for organic farming. Briefly, no synthetic fertilizer was applied to the soil. Pests and disease control were done using biological controls. Frequent mechanical or manual weeding was conducted to suppress the weeds. The genotypes were allocated to the six field trials according to a non-resolvable row column design or a randomized complete block design (Table A1). In 2017, on each on-farm breeding location 30, and in Kleinhohenheim, 40 open pollinating genotypes were cultivated in three field replicates. In 2018, 16 open pollinating genotypes were assessed on each on-farm breeding location, while 36 genotypes were tested in Kleinhohenheim. Over both years 15 genotypes were assessed on all locations forming the orthogonal core of the experiment. Note that data from all genotypes were analyzed, but the presented results were limited to the 15 genotypes which occurred in all location-by-year combinations. The 15 investigated genotypes were 13 red-colored beets: Akela, Betina, Bona, Bordo AS, BoRu1, Carillon RZ, Cervena Kulata, Detroit 3, Jawor, Monty RZ F1, Nobol, Nochowski, Ronjana, one white genotype, Sniezna Kula, and one yellow-colored beet, Burpees Golden. BoRu1 is a breeding line of the breeding organization, Kultursaat e.V., which was used in order to assess its market potential as a new open-pollinated genotype. Table 1 depicts the information regarding shape and color of the beets as well as the seed origin of the 15 selected genotypes. The seeds needed for all six trials were ordered in spring 2017 and stored in a dark and dry place until the sowing day. At De Beersche Hoeve, the field experiments were conducted as row-column design in 2017 and as randomized complete block design in 2018. In the first trial year, the field was fertilized with compost equivalent to 60 kg N ha −1 and in the second year, compost equivalent to 90 kg N ha −1 was incorporated into the soil. The preceding crops were clover grass and Romanesco broccoli in 2017 and 2018, respectively. The field was ploughed two months before each experiment to a depth of 20 cm. Harrowing was done twice before sowing to a depth of approximately 3 cm. The plot size was set as 9 m 2 with three cultivation rows per plot and an inter-row spacing of 50 cm. In the first field trial, sowing was conducted using a precision seed drill on 14 July 2017, and the harvest was done manually on 26 October 2017. In 2018, the beetroots were seeded on 12 June 2018 and harvested on 27 September 2018. In both trials, weed control was done manually and mechanically by hoeing.
At the Horticulture station Heinze, the field experiment was carried out as randomized complete block design in both 2017 and 2018 with three replicates. In 2017, the preceding crop was broccoli, while in 2018, the preceding crop was carrot. In both experimental years, no fertilizer was applied to the field. The field was prepared by ploughing in autumn to a depth of 25 cm and harrowed three months before sowing to a depth of 8 cm. The plot size was 8.75 m 2 with three rows per plot and an inter-row spacing of 42 cm. In 2017, the sowing took place on 08 June 2017 using a hand-sowing machine and the harvest was carried out manually on 30 September 2017. In 2018, the trial started on 09 June 2018 and the beetroots were harvested on 10 September 2018. Weeding was performed by hand and by a wheel hoe.
At the research station Kleinhohenheim, the experimental design was set up as row-column with three replicates in 2017 and as non-resolvable block design, with a block size of ten and a number of treatment of 36, in 2018, with three replicates for 32 genotypes and six replicates for four genotypes. The unequal replications were due to the point that the four genotypes with six replicates belonged to breeding lines, and their performance in the field was unknown. Therefore, in order to have sufficient samples, more replicates were considered for them. In 2017, the preceding crop was clover grass which was considered as green manure, and additionally, the field was fertilized with Vinasse, equivalent to 8.3 kg N ha −1 , one month after seeding. In 2018, no additional fertilizer was used due to the N fixation by the pre-crop clover grass. In January 2017, the field was ploughed to a depth of 20 cm and the seed bed was prepared using a rotary harrow to a depth of 6 cm. In the second year of the experiment, ploughing was carried out in January 2018 and harrowing was done about 40 days before sowing. The plot size was 14 m 2 with four rows per plot and an inter-row spacing of 35 cm. In the first experiment, sowing was done on 06 July 2017 and beetroots were harvested manually on 16 October 2017. The sowing and harvest date in the second field trial was on 06 June 2018 and 10 September 2018, respectively. In both years, sowing was conducted using a pneumatic single-seed drill. Inter-row weed control was performed manually and intra-row weeding was carried out mechanically by hoeing.

Total and Marketable Yield
After the manual harvest of the beetroots, the stems and leaves were cut by hand and immediately the beets' fresh-weight was measured. For the determination of the marketable yield, the harvested beetroots were sorted based on our defined diameter (between 5 and 13 cm) and quality (without any deformation and notable damage).

Diameter of Beet and Width of Leaves-Growth-Base
The average beet size of each genotype as well as the width of the leaves-growth-base was determined by using a Vernier caliper. Five randomly selected marketable beetroots of each plot were measured at the widest part of the beet and at the leaves-growth-base.

Individual Beet Weight
Five randomly selected marketable beets were weighed individually, directly after the harvest, to estimate the average weight of beets of the investigated genotypes and to understand if the yield consisted of a lot of small beets or of bigger and fewer beets per defined area.

Beet Evaluation
With respect to the remarkable influence of the beet morphology and disease resistance of different genotypes on the quality and customer acceptance, the following five traits were evaluated, smoothness of skin, corky surface, root tail, scab (Streptomyces scabies), and uniformity ( Table 2). In this regard, the marketable beetroots were examined optically using a 1-9 scale. Smoothness and corky surface of beet skin, as well as scab incidence, are among the traits affecting the attractiveness of the beets especially for the fresh market. Information on root tail and uniformity is particularly important for the use of the beet for processing purposes.

Statistical Analysis
According to the design of the experiment, the data was analyzed by the following linear mixed model: where a j and l k are the fixed effects of the j-th year and the k-th location. b ljk , r mjk , and c njk are the random effect of the l-th block, m-th row, and n-th column within a year-by-location combination. Note that row and column effects were only fitted for row-column designs. Further note that block effects were only fitted if the design has blocks. Dummy variables were used to block out these effects for all other experiments. τ i is the fixed effect of the i-th genotype. (al) jk is the fixed interaction effect of the j-th year and the k-th location. (τa) ij ,(τl) ik , and (τal) ijk are the random interactions between the corresponding main effects. e ijklmn is the error of y ijklmn with year-by-location specific error variance. The assumptions of normality and homogeneous variances of residuals were checked graphically. If required, a logarithmic transformation of the data was used prior to analysis and results were back-transformed for presentation purpose only. Standard errors were back-transformed using the delta method. Depending on the significance of fixed effects, mean comparisons were done for the highest significant interaction term of all factors using Fisher's Least Significant Difference (LSD) test [23]. Means (or medians in case of back-transformed values) plus minus standard error of the mean (approximated standard error of the median) were presented. A letter display was used to present significant differences between means or medians [23]. Means followed by at least one identical letter were not significantly different from each other.

Results and Discussion
The statistical analysis showed that all investigated traits except the root tail and scab incidence were significantly influenced by the genotype (p < 0.05). This demonstrated a high variability between the genotypes regarding the agronomic traits. In order to assess the performance of different genotypes under different environmental and climatic conditions, the influence of location and experimental year on different agronomic traits was taken into account, as well. In accordance with the statistical analysis, all traits except scab, and leaves-growth-base width and individual beet weight were significantly influenced by the interaction of location × year (p < 0.05). Scab incidence was neither affected by year nor by location. In terms of leaves-growth-base, the results of the analysis of variance (ANOVA) showed a significant effect of year (p < 0.05) and individual beet weight was significantly impacted by both year and location. Outcomes of the ANOVA of the investigated agronomic traits are presented in Tables 3 and 4. Table 3. Mean values and analysis of variance (ANOVA) of results of total yield (t ha −1 ), marketable yield (t ha −1 ), beet diameter (mm), leaves-growth-base width (mm), and individual beet weight (g) of 15 different genotypes of beetroot grown in three research stations within the trial year 2017 and 2018. Results represent the mean values ± standard error. Means followed by at least one identical letter were not significantly different from each other.

Genotype
Total Yield (t ha −
The marketable yield of the studied beetroot genotypes varied between 27.92 ± 3.55 t ha −1 and 44.96 ± 3.50 t ha −1 for the genotypes Burpees Golden and Carillon RZ, respectively, and were statistically different from each other (Table 3). Although the second-highest total yield belonged to the genotype Monty RZ F1, the same order was not seen for the marketable yield. The genotype Bordo AS, which gained the fourth-highest total yield with 50.34 ± 3.97 t ha −1 , had the second-highest marketable yield with 42.29 ± 4.17 t ha −1 , which indicated its higher share of marketable yield compared to the F1 hybrid genotype, Monty RZ F1.
The results revealed that the highest total and marketable yield belonged to the cylindrical genotype, Carillon RZ. Although the spherical shape has been better accepted in the fresh market [24], an increasing use of cylindrical shape genotypes in beetroot processing, due to the higher percentage of utilizable uniform slices in canning and slicing practices, has been reported [25]. Nevertheless, this uniformity in the top, middle, and bottom of the cylindrical-shape beets can be affected by population density and genotype [25]. There are different factors which affect the yield, such as sowing [26] and harvesting [27,28] date, water availability [29], cultivation system [30], and plant density [26]. Feller and Fink [27] stated that the effect of sowing date on the beetroot yield is significant. Takács-Hájos and Rubóczki [31] studied the differences in beetroot yield, depending on sowing date (April, July, and August), and noted that the highest yield was reached at the earliest sowing date in April. Furthermore, Akhiyarov et al. [32] reported a total yield ranging from 35.30-53.70 t ha −1 for a beetroot cultivar sown at five different dates, between the beginning of April and end of May over three years. In the current research, despite the fact that the sowing dates were either in June or July, the total yield was in a similar range as in the study of Akhiyarov et al. [32]. The study of Kikkert et al. [26] depicted that the marketable yield increased significantly with delayed harvest dates. Among all three harvest dates that were investigated for one beetroot genotype within two years, beetroots with 120 days growing period indicated a marketable yield of 48.9 t ha −1 , while those with 100 and 80 growing days revealed a marketable yield of 44.2 and 28.7 t ha −1 , respectively. In the present study, the range of the growing period was between 93-107 days, the marketable yield for all investigated genotypes varied between 27.92 ± 3.55 t ha −1 and 44.96 ± 3.50 t ha −1 .
Considering the cultivation system as one of the yield influencing factors, it was reported that the yield of beetroots grown under organic cultivation systems was 75% of that for the conventional system [30]. Conventionally cultivation of seven different beetroot genotypes in Austria in 2013 showed a yield range of 25.5-32.5 t ha −1 [33]. However, in the present work, despite the fact that the beetroots were grown under organic conditions, total yields were considerably higher. Furthermore, according to Stagnari et al. [29], water stress can significantly decrease the yield of beetroot. In 2018, due to the extremely hot and dry summer in Germany, lower average yields in both locations was noted.

Beet Diameter and Leaves-Growth-Base Width
Except for the cylindrical genotype, Carillon RZ, which significantly had the lowest beet diameter with 48.54 ± 0.31 mm, the range of the beet diameter for all spherical genotypes varied between 71.25 ± 0.31 mm and 81.50 ± 0.31 mm belonging to Burpees Golden and Bordo AS, respectively, which were significantly different from each other (Table 3). Ijoyah et al. [34] compared the beet width of five different spherical-shape beetroot genotypes, and reported a diameter range of 3.64-7.55 cm. Moreover, the beet diameter of the cylindrical-shaped genotype, Regulski Cylinder, which was grown organically in Poland, was 5.43 cm in 2006 and 2008, and 5.0 cm in 2009 [7], which in all cases are slightly larger than the mean beet diameter of the studied cylindrical-shape genotype in the present study. This may be due to the difference in the investigated genotype and the cultivation period, which was longer in the study of Szopińska and Gawęda [7]. The genotypes Betina and Detroit 3 with 80.59 ± 0.31 mm and 80.55 ± 0.31 mm, respectively, had the second and third largest beet diameter. Following the yellow beetroot genotype, Burpees Golden, the cultivars Ronjana with 75.93 ± 0.31 mm and Nobol with 76.31 ± 0.31 mm had the smallest beet diameter ( Table 3). The size of the beet can be mainly influenced by nitrogen fertilization [35] and water availability [26]. In the present work, in accordance with the outcomes of ANOVA, beet size was significantly impacted by genotype and the interaction of location × year (p < 0.0001). Burpees Golden, which possessed the smallest beet diameter, was noticed as the genotype with the lowest total and marketable yields. However, the same relation was not seen for the genotype Monty RZ F1, which was categorized among the smallest-sized beets, but resulted in the second-highest total and marketable yields. Based on Ugrinovic [35], among the marketable beetroots (beet diameter between 6-15 cm), those with the beet diameter of 6-9 cm can be categorized as the first-class beets for the marketable purposes. Accordingly, all 14 spherical-shape genotypes in this study can be considered as first class beets.
The smallest leaves-growth-base width belonged to the cylindrical genotype, Carillon RZ, with 23.63 ± 1.84 mm, followed by the genotype Boru1 and Nobol with 25.22 ± 1.84 mm and 25.78 ± 1.84 mm, respectively. No significant difference was observed between them (Table 3). In contrast, the genotype Nochowski with 43.14 ± 1.84 mm had a significantly larger leaves-growth-base compared to all studied genotypes (Table 3), followed by Betina and Bordo AS with 37.78 ± 1.84 mm and 36.21 ± 1.84 mm, respectively. This trait is important in the processing of beetroot, as the leaves-growth-base is included in the part which should be cut during processing. Therefore, a smaller leaves-growth-base width is desirable in order to achieve less flesh loss. According to the outcomes of this study, as the genotype demonstrated a significant effect on the leaves-growth-base width (p < 0.0001), selection of the suitable genotype, especially for processing purposes, may lead to lower waste.

Individual Beet Weight
According to the ANOVA, genotype, experimental year, and location influenced the weight of the beetroot. The beets weight varied between 214.71 ± 17.42 g and 283.33 ± 17.42 g ( Table 3).
The fresh beet weight of 15 beetroot genotypes, which were sown in July in Hungary, ranged between 88.93-198.67 g, which is notably lower compared to the findings of the present study [31]. Furthermore, investigation of seven beetroot genotypes in Austria demonstrated a beet weight range of 214-366 g, in which the minimum value was similar to the results of this study but the maximum value was higher. Following the same position in beet diameter as well as the total yield, the yellow genotype, Burpees Golden, which owned the smallest beet diameter and the lowest yield, showed significantly lowest value regarding the beet weight with 214.71 ± 17.42 g. On the other hand, the genotype, Bordo AS, with the largest beet diameter, possessed the heaviest beets among all investigated genotypes with a mean beet weight of 283.33 ± 17.42 g. Specifications of the market regarding the size of the root crops have been increasing [36]. Although the agronomic factors affecting the size, such as plant density or harvest time, can be controlled, the size variation within a crop is still great [36]. Therefore, in order to achieve a comprehensive overview of the beet size, beside the beet diameter, the fresh weight of the individual beet was measured in this study. An important factor that determines beet size is the water availability during the growing period [26], as earlier sowing dates showed higher beet weights related to water supply in 15 different genotypes cultivated in Debrecen, Hungary [31]. In the current study, all the field trials were started either in June or July. In addition, the second year of the experiment was considered as an extreme climate condition due to the high temperature and low precipitation level. Therefore, considering the significant effect of the trial year on the beet weight in this study, it can be concluded that the abiotic factors impacted the beet weight.
Due to the importance of the beet size for the fresh market as well as for processing, information about the genetically determined beet size in order to select the suitable genotype for the specific utilization should be taken into account. For instance, the genotypes with beet sizes too large for use in the fresh market could be interesting for juice production or the use in the production of pre-cooked beetroot.

Beet Evaluation
Morphology of the beet can have a notable influence on the quality, and its attractiveness is a key factor for the fresh market [24]. It is reported [37] that besides high yield and resistance against diseases and pests, shape [38] and uniformity in color without white flecking [39] are the main criteria in the selection process of a beetroot genotype for breeding [40] and processing purposes. One of the factors which plays a key role in physiology and consequently for the quality of different plant species is drought stress [29]. Accordingly, due to the high temperature as well as low precipitation in the second year of this experiment, it can be assumed that the beet morphology has been influenced by this factor.

Smoothness and Corky Surface of Beet Skin
Evaluation of the beets revealed that the skin of the 15 studied genotypes differed significantly from rough to smooth with a range of 3.96 ± 0.40 to 6.35 ± 0.40 (Table 4). The results depicted that the white genotype, Sniezna Kula, with 6.35 ± 0.40, had the smoothest skin among all evaluated genotypes followed by Jawor and Carillon RZ with 6.06 ± 0.40 and 6.03 ± 0.40, respectively. There was no significant difference noted between them (Table 4). On the other hand, Detroit 3, Nobol, and Burpees Golden with 3.96 ± 0.40, 4.58 ± 0.40, and 4.59 ± 0.41, respectively, were evaluated as three genotypes with the lowest smoothness, although they were not significantly different from each other (Table 4). According to the results of the analysis of the variance showing a significant effect of genotype on skin smoothness (p = 0.0024), high variability of the beetroot genotypes regarding these characteristics can be concluded. Furthermore, the interaction of location × year significantly affected the smoothness of the beet skin as well (p < 0.0001).
Beetroots with lower corky skin have a more attractive appearance for the fresh market. Corky surface is resulted from the light brown stripes which are formed on the skin of the beet. Nobol, with 5.97 ± 0.34, had a significantly corkier surface compared to Cervena Kulata, with 5.04 ± 0.34, which was evaluated as the genotype with the second-corkiest surface. Contrarily, Sniezna Kula demonstrated the least corky surface with 2.85 ± 0.34 followed by the genotype Nochowski with 4.08 ± 0.34. The values were significantly different from each other (Table 4).
Comparing the results of the two traits, corky surface and beet skin smoothness, it can be seen that the genotype Sniezna Kula, which had the smoothest beet skin, was as well the genotype with the lowest corky surface which may be due to the point that the lower cork on the beet makes the skin smoother. Likewise, Nobol with the second-lowest skin smoothness was noted as the genotype with the highest corky surface. The study of Da Silva et al. (2006) indicated that an increase of the mean tuber weight in potato finally led to a smoother skin [41]. In the current study, the same correlation was not noted for beetroot. Furthermore, it has been stated that one of the factors which improves skin smoothness of carrots is to grow them on ridges [42]. It was reported that the incidence of carrot cavity spots is more common under the wet conditions [42] and due to the fact that the soil on ridges are often softer and dries faster [43] by reason of better aeration, carrots grown on the ridges showed less cavity spots and better root shape and skin smoothness. Similar cultivation methods may improve the skin smoothness of beetroot.

Root Tail
The beetroots were evaluated optically based on the degree of detachment of the root tail from the beet. The results showed that the genotype Akela, with 6.67 ± 0.42, was ranked as the genotype with the highest degree of detachment of the root tail from the beet followed by Ronjana and Detroit 3 with 6.56 ± 0.42 and 6.45 ± 0.42, respectively. The cylindrical genotype, Carillon RZ, with 5.0 ± 0.42, was ranked as the genotype with the lowest detachment of the root tail from the beet compared to the other studied genotypes. Although the value was not significantly different from that of the cultivars Burpees Golden and Betina, which were rated around the average level with 5.23 ± 0.42 and 5.32 ± 0.42, respectively (Table 4). There were no significant differences between genotypes analyzed (Table 4). However, the statistical analysis of the results demonstrated a significant influence of the interaction of location × year (p < 0.0001) on the root tail. This is in line with the findings of previous studies [44,45] which noted the impact of environmental conditions, especially soil type, on the size and shape of the root vegetables. Moreover, the shape of the beetroot as well as the root tail is considered as a major factor due to the loss which may occur during processing [24]. To be more precise, the higher the degree of detachment of the root tail from the beet, the lower loss will result.

Scab (Streptomyces Scabies)
According to Table 4, scab incidence ranged between 3.20 ± 0.34 and 4.65 ± 0.34, which based on the defined 1-9 scale in Table 2, it can be concluded that all evaluated genotypes demonstrated little to average scab formation. Sniezna Kula and Ronjana, with both 4.65 ± 0.34, depicted the highest incidence of scab among all 15 genotypes. The cylindrical genotype, Carillon RZ, with 3.20 ± 0.34 and the F1 hybrid, Monty RZ F1, with 3.43 ± 0.34, were the most resistant genotypes to scab occurrence. There were no significant differences between the genotypes analyzed (Table 4). Scab (Streptomyces scabies) is a bacterial disease which causes warty growth on the skin and produces corky spots on the surface of different root crops, including beetroot and potato [46,47]. Despite the fact that the scab arises in all types of beet (Beta vulgaris L.), the occurrence of this corky lesion in beetroot due to the use in the fresh market as well as pre-cooked product in transparent packaging (such as glass bottles), is more undesirable compared with sugar beet which is pulped [48]. Crop rotation and ensuring the soil pH close to neutral are two major ways to prevent scab emergence [47]. A study by Lapwood et al. [48] on 12 different beetroot genotypes, including ten spherical and two cylindrical, showed that the two cylindrical genotypes were among the most resistant ones against the scab incidence. In line with these results, our research has as well depicted that regarding scab disease, the open-pollinated cylindrical genotype, Carillon RZ, is as resistant as the F1 hybrid, Monty RZ F1, which is typically known as the genotype with the highest resistance to the biotic and abiotic stressor.

Uniformity
The study showed that, based on the defined ranking, using 1-9 scale (Table 2), all 15 genotypes had an average to high uniformity. The highest uniformity was noted for the genotypes Sniezna Kula, Monty RZ F1, and Betina with 6.65 ± 0.38, 6.61 ± 0.38, and 6.28 ± 0.38, respectively. Detroit 3 with 5.19 ± 0.38, Burpees Golden with 5.27 ± 0.38, and Bordo AS with 5.35 ± 0.38 were scored as the genotypes with the lowest uniformity, through which the results were not significantly different from each other ( Table 4). Strang et al. [49] reported a uniformity of the beet shape for Burpees Golden as 6.5 out of 10 (equivalent to 5.8 in our rating scale), which is close to the rate which has been noted in the present study. Due to the influencing factors such as growing space, and balanced and sufficient nutrients and water supply, uniformity of the beetroot in terms of shape and size is difficult to attain [26]. One of the influencing factors on the shape uniformity of the beets is uniform water supply in the soil [31]. Interim water deficiency can change the shape of the beet from the common spherical to oval or heart-shape, which is not desirable for the fresh market or processing purposes [24]. In the present study, due to the dry summer and low precipitation in 2018, it can be assumed that the lower water availability in the second year of the experiment affected the mean beet uniformity. Another major factor leading to non-uniformity in the maturity stage as well as the size of beets at harvest was reported to be the non-uniform seed germination [50] which occurs in the early stage of crop development [51,52]. Likewise, the cluster shape of the beetroot seed which results in up to three seedlings, causes the root size variability [53]. Furthermore, the distribution of the plants within the row and row width have been reported as two other factors affecting the uniformity of the plant [54]. In this regard, the study of Griepentrog et al. (2009) demonstrated the improvement of the spatial uniformity of sugar beet with the increase of row width from 32 to 50 cm [54]. Despite the point that the row distance width in the present study varied in all three locations (35,42, and 50 cm), the effect of the location on the beet uniformity was not significant (Table 4). However, the interaction of location × year indicated a significant impact on the beet uniformity.

Overall Outcomes
Comparing the open-pollinated genotypes with the F1 hybrid, Monty RZ F1, indicated that in all investigated traits except for yield, scab, and uniformity, the open-pollinated genotypes revealed desirable competitive results. Additionally, with the available gene pool and selection of genotypes with preferable traits, there would be a possibility for cross-breeding in order to incorporate the outstanding characteristics into a new bred line. Considering the positive aspects of the open-pollinated genotypes for the use in organic farming, including maintaining the plant integrity and passing the plant-specific features to the next generation as well as the self-reproductive ability, the outcomes of the present study revealed the significant influence of the genotype on most of the assessed traits. Therefore, by providing information regarding the agronomic performances of different open-pollinated genotypes, the study recommends the selection of suitable genotypes for various purposes. In terms of beet uniformity as well as the persistence against the common beetroot disease, Streptomyces scabies, all studied genotypes indicated a high stability regarding the environmental conditions. Nevertheless, with the significant impact of location on total and marketable yield (p < 0.0001), it can be noted that the investigated genotypes may not be recommended for a broad range of environments if yield would be the only trait of interest.
The traits, high total and marketable yield and high degree of uniformity of the breeding line, BoRu1, showed a potential for the cultivation in organic farming. However, other traits like skin smoothness and scab formation still need further improvement by further breeding efforts.

Conclusions
With investigating the performance of various beetroot genotypes in different locations, it can be concluded that due to the significant genetic variability in beetroot, selection of suitable genotypes should be aligned with the intended final utilization. The cylindrical genotype, Carillon RZ, with the highest total and marketable yield, smallest leaves-growth-base width and one of the most resistance genotype against the common beet disease, scab, can be suggested as a potential genotype with outstanding characteristics to farmers for cultivation in organic farming. Finally, this study only indicated the agronomic traits of the genotypes and the related impact of environmental conditions, therefore, for the application of the investigated genotypes in different sections, the sensory quality as well as the amount of the health-promoting compounds, should be additionally taken into account.