Caryologia. International Journal of Cytology, Cytosystematics and Cytogenetics 75(4): 133-148, 2022 Firenze University Press www.fupress.com/caryologia ISSN 0008-7114 (print) | ISSN 2165-5391 (online) | DOI: 10.36253/caryologia-1633 Caryologia International Journal of Cytology, Cytosystematics and Cytogenetics C i t a t i o n : S a e e d M o h a m m a d - pour, Ghasem Karimzadeh, Seyed Mahmood Ghaffari (2022). Karyomorphol- ogy, genome size, and variation of antioxidant in twelve berry species from Iran. Caryologia 75(4): 133-148. doi: 10.36253/caryologia-1633 Received: April 18, 2022 Accepted: December 04, 2022 Published: April 28, 2023 Copyright: © 2022 Saeed Mohammad- pour, Ghasem Karimzadeh, Seyed Mahmood Ghaffari. This is an open access, peer-reviewed article pub- lished by Firenze University Press (http://www.fupress.com/caryologia) and distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Data Availability Statement: All rel- evant data are within the paper and its Supporting Information files. Competing Interests: The Author(s) declare(s) no conflict of interest. Karyomorphology, genome size, and variation of antioxidant in twelve berry species from Iran Saeed Mohammadpour1, Ghasem Karimzadeh1,*, Seyed Mahmood Ghaffari2 1 Department of Plant Genetics and Breeding, Col lege of Agriculture, Tarbiat Modares University, P.O. Box 14115-336, Tehran, Iran 2 Institute of Biochemistry and Biophysics, University of Tehran, P.O. Box 13145-1384, Tehran, Iran *Corresponding author. E-mail: karimzadeh_g@modares.ac.ir Abstract. Twelve berry species, including Rubus fruticosus cv. Qaemshahr, Rubus occi- dentalis cv. Qaemshahr, and Morus alba cv. Mashhad, Morus rubra cv. Karaj, Fragaria vesca subsp. vesca, Ribes nigrum, Ribes rubrum, Ribes uva-crispa; Lycium barbarum, Lycium infaustum, Lycium ruthenicum, and Vaccinium corymbosum were surveyed for karyomorphology analysis, monoploid genome size, and antioxidant activity. Results indicated that all species were diploid (2n = 2x = 14, 16, 24, and 28). Among these spe- cies, chromosome counts and karyomorphology parameters of three cultivars includ- ing: R. fruticosus cv. Qaemshahr (2n = 2x = 14), R. occidentalis cv. Qaemshahr (2n = 2x = 14), and M. alba cv. Mashhad (2n = 2x = 28) are reported here for the first time. The flow cytometric mean monoploid 2Cx DNA of berry species was 2.35 pg, varied from 0.68 (Rubus occidentalis cv. Qaemshahr) to 5.15 pg (Lycium ruthenicum). The average antioxidant capacity of berry species was obtained 32.8 μmol g-1. Their average total phenol and flavonoid contents were 4.98 mg g-1 (3.08-8.61 mg g-1), and 5.18 (2.47- 10.63 mg g-1), respectively. Keywords: Blueberry, Raspberry, Cytogenetic, Antioxidant, Karyomorphology, Flow cytometry. INTRODUCTION Berries are economically important crop in many countries (Umdale et al. 2020). Interest in berries has recently increased because they are excellent sources of health-promoting vitamins, anti-oxidants, polyphenols, especially anthocyanins, micronutrients, fiber, and other valuable nutrients (Song and Sing, 2004). Berries are low in calories and are high in moisture and fiber (Basu et al. 2010). They contain natural antioxidants such as vitamins C and E, and micronutrients such as folic acid, calcium, selenium, alpha and beta carotene, and lutein (Basu 2019). However, berry production is deficient, even though soil and climate conditions in Iran are excellent for the intensive cul- tivation of berry fruits (Nestby et al. 2019). Some of the berry cultivars that 134 Saeed Mohammadpour, Ghasem Karimzadeh, Seyed Mahmood Ghaffari have recently been introduced to Iran are potentially profitable alternative and non-conventional fruit crops. Among those are cultivars of the blueberry (Vaccinium corymbosum L.), the goji berry (Lycium barbarum, L. infaustum, L. ruthenicum), and the raspberry (Rubus occidentalis) (Ahmadu and Ahmad 2021). In Iran, the production of blackberry (Rubus fruticosus) relies main- ly on the collection of the berries from stands of native varieties that grow wild in the mountains. Habitat destruction has led to a reduction in the supply of native blackberries. Although these native wild varieties are not as suitable for intensive culture as the new cultivars, they are still very valuable high dietary value and breed- ing potential. They can be a versatile raw material for the food processing and pharmaceutical industries. Recently, interest in growing various berry fruits has been increas- ing among either home gardeners, or small farmers. Ber- ry culture is increasing in the country’s mountainous, where soil and climate conditions are more favorable for intensive berry culture than in the lowland areas (Nest- by et al. 2019). The consumption of berry fruits and their contribu- tion to improving cardiovascular health is a significant issue (Blumberg et al. 2013). The commonly consumed berries in the United States, including blackberry, black raspberry, blueberry, goji berry, cranberry, red raspberry, and strawberries (Yang et al. 2019). Display quotations of over 40 words or as needed. Berry fruits are deciduous shrubbery that grows in different parts of the world (Donno et al. 2015). For instance, goji berry grows in China, Tibet, Argentina, Chile, Southern Africa, and other parts of Asia, and its fruits are 1-2 cm long, bright orange-red ellipsoid ber- ries. Goji berry is widely distributed in warm regions of the world, in particular, in the Mediterranean area and Southwest and Central Asia. It is also cultivated in North America and Australia as a hedge plant (Potter- at 2010). Mulberry (Morus alba) is native to China, but is also found worldwide on not native continents, such as Europe, Africa, North America, and South America. Currant (Ribes spp.), is widely cultivated across temper- ate Europe, Russia, New Zealand, parts of Asia, and to a lesser extent North America (Steffen et al. 2015). The main centers of diversity for blueberry are distributed in Europe (Poland, Germany, France, The Netherlands, Spain, and Sweden), New Zealand, Mexico, and North America. Karyotype analysis and chromosome observation are necessary to elucidate phylogenetic relationships, structure, function, organization, and evolution of plant genes and genomes (Amosova et al. 2019). For those rea- sons, many studies must be performed to observe plant chromosomes. Higher similarity in chromosome shapes suggests a closer phylogenetic relatedness. The base chromosome number in Rubus spp. (R. fruticosus, and R. occidentalis), Ribes spp. (R. nigrum, R. rubrum, and R. uva-crispa), Lycieae (L. barbarum, L. infaustum, and L. ruthenicum), Moraceae (M. alba, and M. rubra) is x = 7, 8, 12, and 14, respectively. The base chromosome num- ber in strawberry (Fragaria vesca) and blueberry (Vacci- num corymbosum) is x = 7, and 14, respectively (Zong et al. 2021). Most species are diploid with 2n = 24. Cytolog- ical information, including the number of chromosome and karyotypes, is available for many American and Asian Lycium and several species in Southern African (Bernardello et al. 2008). Monoploid genome size in the form of base-pair was calculated based on the formula proposed by Doležel et al. (2003), when 1 pg of DNA represents 978 mega base pairs (Mbp). Monoploid genome size was considered as the amount of DNA of one chromosome set, 1Cx- value, with chromosome base number x, and holoploid genome size as the amount of DNA of the whole chro- mosome complement, 1C-value, with chromosome num- ber n, irrespective of the degree of generative polyploidy (Greilhuber et al. 2005; Mahdavi and Karimzadeh 2010; Karimzadeh et al. 2011; Abedi et al. 2015; Tavan et al. 2015). In recent decades, significant progress has been made in the use of flow cytometry in various fields of botany. This growth is due to the advantages of flow cytometry such as high speed, ease of sample prepara- tion, and analysis of inactive mitotic tissues. Despite these advantages, the need for fresh plant materials may often prevent the further development of flow cytom- etry in field research (Suda and Trávníček 2006). Flow cytometry is known as the most powerful, reliable, and quick method for estimating the genome size (2C DNA) for a wide range of plant communities (e.g. Garcia et al. 2004; Doležel and Bartoš 2005; Doležel et al. 2007; Loureiro et al. 2007; Mahdavi and Karimzadeh 2010; Karimzadeh et al. 2010, 2011; Abedi et al. 2015; Tavan et al. 2015; Tarkesh Esfahani et al. 2016, 2020; Javadian et al. 2017; Hamidi et al. 2018; Zarabizadeh et al. 2022). The term “C-value” refers to the constant amount of DNA of an unreplicated haploid chromosome comple- ment (Swift 1950). Monoploid genome size was consid- ered as the amount of DNA of one chromosome set, 1Cx-value, with chromosome base number x (Greilhuber et al. 2005). The mean monoploid 2Cx DNA value was 2.35 pg. On the other hand, the monoploid genome size (1Cx DNA) varied from 420.54 Mbp (S3; Rubus occi- dentalis cv. Qaemshahr) to 2518.35 Mbp (S9; Lycium ruthenicum). Reports of 2Cx DNA in some berry spe- cies have been reported, including S3; Rubus occidenta- lis = 0.60 pg, S5; Ribes rubrum = 1.94 pg, and S6; Ribes 135Karyomorphology, genome size, and variation of antioxidant in twelve berry species from Iran uva-crispa = 1.88 pg (Meng and Finn 2002; Chiche et al. 2003). Antioxidants are compounds that effectively prevent oxidation in a variety of ways and, by slowing down the oxidation rate, significantly increase the oxidation peri- od (Dorman et al. 2003). From one perspective, antioxi- dants fall into two main categories: synthetic and natu- ral. Phenolic compounds are one of the best sources of natural antioxidants (Dorman et al. 2003). Because the use of such plant resources is effective in delaying oxi- dation and reducing chronic diseases, mutations, and cancer (Briskin 2000). In new research works in the field of natural antioxidants, much attention has been paid to essential oils and extracts of medicinal plants because a wide range of medicinal plants and their aromatic com- pounds as natural sources with antioxidant properties, (Moure et al. 2001; Tepe et al. 2006). It is generally rec- ognized that oxygen free radicals (ROS) are associated with many potential risks, including Parkinson’s disease, cancer, Alzheimer’s, as well as aging (Liochev 2013; Kim et al. 2016).Berry species are popular around the world as a ‘Super fruit’ due to their nutritional value, elevated levels of bioactive phenolic molecules, and excellent sen- sory evaluation (Kalt et al. 2020). In addition to these essential nutrients, berries contain a wide range of anti- oxidant phenolic molecules such as phenolic acids, fla- vonoids, flavonols, flavanols, anthocyanins, proantho- cyanidins, and ellagitannins (Prior et al. 1998; Piljac- Žegarac et al. 2009; Nile and Park 2014; Skrovankova et al. 2015). Cranberries ranked first in polyphenol content among the commonly consumed fruits in North Amer- ica which relate to high antioxidant activity (Vinson et al. 2001). The Consumption of lingonberries and bilber- ries has been proved in preventing human cancer which is ascribed to the high levels of phenolic and anthocya- nin compounds (Faria et al. 2010; Yang and Kortesnie- mi 2015). The average total antioxidant capacity, total phenol content, total flavonoid content of berries were 32.79, 4.98, and 5.18 μmol g-1, respectively. In the study of Islam et al. (2017) on genus Lycium, the mean of the antioxidant activity for red goji berry (S7; Lycium barba- rum) and black goji berry (S9; Lycium ruthenicum) were measured as 16.65 and 34.28 µmol g-1, respectively. In the report of Mustafa Ahmed et al. (2022) on genus Vac- cinium, the mean of total phenol content, total flavonoid content, and antioxidant capacity were reported as 3.73 mg g-1, 2.57 mg g-1, and 17.67 µmol g-1, respectively. In the present study, we provide a detailed karyo- type analysis for 12 berries species. Our goals were to confirm the number of chromosomes, karyotype deter- mination, estimate the nuclear genome size, assess the phenolic profile, antioxidant properties, and flavonoid content of these twelve species, and provide scientific insight into the phenolic and antioxidant functions of these species to consumers and nutraceutical industry. MATERIALS AND METHODS Plant materials Seeds were kindly provided by the Cultivation Development Company of Berries (Mazandaran, Iran). The culture medium used was MS (Murashige and Skoog 1962) with macroelements at 1/3 (Rache and Pacheco 2010), 1/8, and 1/16 of its original concentra- tion (NH4NO3: 1650 mg l-1, KNO3: 1900 mg l-1, KH2PO4: 170 mg l-1, CaCl2H2O: 440 mg l-1 and MgSO4.7H2O: 370 Table 1. Specifications of berry species collection sites (Berry spp.) used in this research. Species codes English name Scientific name Localities Longitude (E); Latitude (N) Altitude (m) S1 Strawberry Fragaria vesca subsp. vesca Iran, Mazandaran, sari 52°59’42.83”; 36°32’25.67” 83 S2 Blackberry Rubus fruticosus cv. Qaemshahr Iran, Qaemshahr 52°51’28.29”; 36°27’19.11” 282 S3 Black raspberry Rubus occidentalis cv. Qaemshahr Iran, Qaemshahr 52°51’28.29”; 36°27’19.11” 282 S4 Black currant Ribes nigrum USA, Arizona, Pinal 111°17’04.21”; 32°48’58.34” 660 S5 Red currant Ribes rubrum USA, Arizona, Pinal 111°17’04.21”; 32°48’58.34” 660 S6 Goose berry Ribes uva-crispa USA, Arizona, Pinal 111°17’04.21”; 32°48’58.34” 660 S7 Red goji berry Lycium barbarum USA, Arizona, Pinal 111°17’04.21”; 32°48’58.34” 660 S8 Goji berry Lycium infaustum USA, Arizona, Pinal 111°17’04.21”; 32°48’58.34” 660 S9 Black goji berry Lycium ruthenicum USA, Arizona, Pinal 111°17’04.21”; 32°48’58.34” 660 S10 Blueberry Vaccinium corymbosum USA, Maine 68°59’05.21”; 45°11’18.34” 426 S11 White mulberry Morus alba cv. Mashhad Iran, Karaj, Mohammadshahr 51°00’13.27”; 35°44’23.12” 1228 S12 Red mulberry Morus rubra cv. Karaj Iran, Karaj, Mohammadshahr 51°00’13.27”; 35°44’23.12” 1228 136 Saeed Mohammadpour, Ghasem Karimzadeh, Seyed Mahmood Ghaffari mg l-1). The medium was adjusted to pH = 5.8 and auto- claved for 20 min at 121ºC. These plants were grown in a greenhouse at the college of Agriculture, Tarbiat Modares University (TMU). Voucher specimens of the examined species are preserved in the herbarium of the Research Institute of Forests and Rangelands of Iran (TARI) and the Tarbiat Modares University of Iran (TMU). Also, the images of species are shown in Figure Figure 1. The picture of berry species that studied in this research. 137Karyomorphology, genome size, and variation of antioxidant in twelve berry species from Iran 1. All cultures were maintained in a growth chamber at 25 ºC, using a 16 h light/8 h dark photoperiod. The light was supplied, using white fluorescent lamps at an inten- sity of 50 µmol m-2 s-1. Karyomorphology The best method to study mitosis for preparing a karyotype is root tip meristem tissue, which was used for cytogenetic studies. For the cytological preparations, growing roots about 1 cm long were removed at the appropriate time (when the largest number of cells are in metaphase) and pretreated in 0.002 mol 1−1 8-hydrox- yquinoline for 4h at room temperature (RT) in the dark to induce metaphase arrest, followed by washing twice (each for 5 min) in distilled H2O and fixing in fresh 3:1 (v/v) absolute ethanol: glacial acetic acid for 24 h at 4 °C. The fixed roots were hydrolyzed in 1 M HCl for 7-10 min at 60°C washed two times (each for 5 min) in distilled H2O and stained in 1% (w/v) aceto-orcein for 3 h or in 4% (w/v) Hematoxylin for 4 h at RT. Five well-spread monolayer metaphase plates from different individuals were analyzed per species. High- resolution microscopic digital photographs were taken, using an Olympus BX50 microscope (Olympus Optical Co., Tokyo, Japan), equipped with an Olympus DP12 digital camera (Olympus Optical Co.). Eight chromo- somal parameters were either measured as long (L) and short (S) arm lengths or calculated as chromosome length (CL), arm ratio (AR; L/S), r-value (S/L), total chromosome volume (TCV = πr2 CL), where r = aver- age chromosome radius, percentage of total chromo- some form (F%), and centromeric index (CI = S/CL). Idi- ograms were drawn from the mean values, and chromo- some types were determined, using the formula of Levan et al. (1964). For karyotypic analysis, five parameters, including karyotype total form percentage (TF% = ƩS/ ƩCL × 100), coefficient of variation of total chromosome length (CVCL% = (total CL standard deviation/total CL mean) × 100), mean centromeric asymmetry (MCA = A × 100), where the degree of karyotype asymmetry (A = [Ʃ(L-S/L+S(/n] × 100(, Stebbins (1971) asymmetry cate- gories (ST), and Romero-Zarco (1986) indices: intrachro- mosomal asymmetry index (A1) and interchromosomal asymmetry index (A2), were measured. Flow cytometric genome size estimation Relative DNA content was determined, using PI- stained flow cytometry. Young leaves of the analyzed individuals and a reference standard (either Solanum lycopersicum; 2C = 1.96 pg., Petroselinum crispum; 2C = 4.45 pg., or Zea mays; 2C = 5.43 pg) were co-chopped with a razor blade in a glass petri dish, containing 1 ml of ice-cold WPB buffer (Woody Plant Buffer; 0.02 M Tris-HCl, 4 mM MgCl2.6H2O, 2 mM EDTA Na2.2H2O, 86 mM NaCl, 10 mM Sodium metabisulfite, 1% PVP- 10, and 1% (v/v) Triton X-100, pH 7.5) (Loureiro et al. 2007). The crude nuclei suspension was filtered through a 50-µm nylon mesh. RNase and propidium iodide (each 50 µg ml−1) was then added. After incubation for two min at RT, the relative fluorescence intensity of nuclei was analyzed. At least 5000 nuclei were analyzed in each measurement with CV% (coefficient of variation; %) val- ues below 5.0%. Subjects were randomly selected for the flow cytometric analysis. The DNA amount of a sample was calculated based on the values of the G1 peak means (Doležel et al. 2003, 2007; Doležel and Bartoš 2005) as follows: Sample 2Cx DNA (pg) = (Sample G1 peak mean/Stand- ard G1 peak mean) × Standard 2C DNA (pg). The obtained data were initially checked for the nor- mality test and analyzed, using SPSS (Version 22) and Minitab 17. The karyotypic, flow cytometric, and phyto- chemical content data were first tested for the normality and then analyzed based on a completely randomized design (CRD) with 5, 3, and 3 replications, respectively. The results were statistically evaluated by analysis of variance (ANOVA) and expressed as mean. Means com- parisons were performed, using the LSD test with SPSS; differences were considered statistically significant at P ≤ 0.01 and P ≤ 0.05. Linear regression analysis was carried out to find out the relationship between monoploid 2Cx DNA values and some traits, using Minitab 17. Chromatographic analysis The method used to determine total polyphenol con- tent (TPC) was based on Folin-Ciocalteu phenol reagent and spectrophotometric determination at 765 nm. The method used to determine the total flavonoid content (TFC) was based on aluminum chloride and spectropho- tometric determination at 420 nm. Antioxidant activity in the berries fruit pulp was evaluated by a 2,2-diphenyl- 1-picrylhydrazyl (DPPH) assay (Laczkó-Zöld et al. 2018). Total Antioxidant Capacity (TAC) The scavenging activity for DPPH radicals was determined using spectrophotometric analysis accord- 138 Saeed Mohammadpour, Ghasem Karimzadeh, Seyed Mahmood Ghaffari ing to a modified method of Lalegani et al. (2018). Brief- ly 100 µl sample solution was added to 900 µl of 0.1 M DPPH solution and mixed thoroughly at RT. Absorbance at 517 nm was determined after 30 min. Total Phenolic Content (TPC) The TPC was estimated, using the Folin- Ciocalteu colorimetric method described by Barreca et al. (2016). A 20 µl of methanolic extract was added to 2 ml of deionized water, and 100 µl of dilutions Folin- Ciocal- teu was then added. After 1-8 min (average 5 min), 300 µl of Sodium carbonate 7% (w/v) was added to it, and the absorbance was measured in the dark at 765 nm after incubation for 2 h at RT. Quantification was done based on the standard curve of gallic acid. Results were expressed as equivalent of the gallic acid (GAE), i.e., mg gallic acid g-1 DW. Total Flavonoid Content (TFC) The TFC was estimated, using the aluminum chlo- ride method described by Barreca et al. (2016). In this method, a 600 µl of methanolic extract was added to 600 µl of 2% (w/v) Aluminum chloride, and after 10 min, the absorbance was measured at 420 nm. Quantification was done based on the standard curve of quercetin. The TFC was calculated and expressed as quercetin equivalents, i.e., mg quercetin g-1 DW. RESULTS Karyomorphologic analysis In the current study, one chromosome type “m” formed karyotype formulas of “14m” for Fragaria vesca subsp. vesca (S1), Rubus fruticosus cv. Qaemshahr (S2), and Rubus occidentalis cv. Qaemshahr (S3), “16m” for Ribes nigrum (S4), Ribes rubrum (S5), and Ribes uva- crispa (S6), “24m” for Lycium barbarum (S7), Lycium infaustum (S8), Lycium ruthenicum (S9), and Vaccinium corymbosum (S10), and “28m” for S11 and S12 (Morus alba cv. Mashhad (S11), Morus rubra cv. Karaj (S12). The base chromosome number in Fragaria vesca subsp. vesca, Rubus fruticosus cv. Qaemshahr, Rubus occidenta- lis cv. Qaemshahr; is x = 7; Ribes nigrum, Ribes rubrum, Ribes uva-crispa x = 8; Lycium barbarum, Lycium infaustum, Lycium ruthenicum, Vaccinium corymbosum x = 12; Morus alba cv. Mashhad, Morus rubra cv. Karaj x = 14, and all of these species are diploid with 2n = 14, 16, 24, and 28. ANOVA shows significant (P < 0.01) differences between the berries species for the most studied chro- mosomal parameters (Table 2). Karyotypes of somatic complement and the idiograms of the haploid comple- Table 2. ANOVA for chromosomal parameters of berry species. S.O.V. DF MS S L CL AR r-value F% TCV CI Species 11 7.413** 8.825** 20.685** 0.128** 0.128** 23.661** 49.778** 0.199** Error 593 0.039 0.046 0.11 0.017 0.017 0.121 0.422 0.057 Total 604 **Significant difference at P < 0.01 Figure 2. Somatic chromosomes of 12 berry species. Scale bar = 5 μm. 139Karyomorphology, genome size, and variation of antioxidant in twelve berry species from Iran ment of studied berries are demonstrated in Figures 2 and 3, respectively. The mean value of chromosome length (CL) was 2.83 μm, varying from 1.06 μm (Morus rubra cv. Karaj; S12) to 3.75 μm (Ribes nigrum; S4). The mean TCV was 1.6 μm, ranging from 0.10 μm (S10) to 6.19 μm (S1). The mean CI of the complement var- ied from 46% (S12) to 51% (S5). The studied berry spe- cies showed different symmetrical groups according to various karyotypic symmetrical indices. For instance, the highest value of total form percentage of karyotype (TF%) was detected in S5 (50.51%; Table 3; the most symmetric), and the lowest value was identified in S12 (46.28%; the most asymmetric). The highest and the lowest values of coefficient of variation (CVCL%) were identified on S11 (18.03%; the most asymmetric) and S3 (11.98%; the most symmetric), respectively. On the other hand, karyotypes of all these species were classi- fied in 1A class of Stebbins classification (Stebbins 1971; Table 3). The basic chromosome number for these ber- ries is consistent with previously published data. The ploidy level agreed with prior and this study’s f low cytometric data. Despite rapid and convenient estima- tion of the ploidy level, using monoploid genome size measurements, chromosome counting is necessary to reliably infer the ploidy level when the DNA content of polyploids fails to match the ratio expected from closely related diploids. For more detailed studies of asymmetry, Romero-Zarco’s (1986) indices of A1 and A2 were also calculated (Table 3; Figure 4). Table 3 also shows two more asymmetry indi- ces of MCA and CVCL. The A1 and A2 indices illustrate four groups. The scatter plot of MCA vs CVCL illustrates three groups. To determine total variation in berries and parameters quota in total variation, the principal com- ponent analysis (PCA) of the karyotypic parameters was performed. The single dendrogram constructed based on karyotype similarities (Figure 3) shows four major clus- Figure 3. Idiograms of haploid chromosomes of 12 berry species. Table 3. Karyotypic symmetry formula and flow cytometric genome size of 12 berry species. Total form percentage of karyo- type (TF%), mean centromeric asymmetry (MCA%), coefficient of variation of chromosome length (CVCL%), karyotype formula (KF), Stebbins asymmetry categories (ST), intrachromosome asymmetry index (A1), interchromosome asymmetry index (A2), m: metacen- tric, M: metacentric. Means with the same symbol letter in the “2Cx DNA (pg)” column are not significantly different (P > 0.01), using LSD test. Species Asymmetry indices (Romero-Zarco, 1986) ST KF CVCL% MCA% TF% 2Cx DNA (pg) A1 A2 S1 0.46 0.13 1A 14m 13.38 4.37 47.79 2.10d S2 0.47 0.18 1A 14m 17.71 4.91 47.59 1.55g S3 0.47 0.12 1A 14m 11.98 4.73 47.66 0.86i S4 0.39 0.14 1A 16m 14.21 4.31 47.86 1.86e S5 0.33 0.16 1A 16m 16.17 4.54 50.51 1.94e S6 0.39 0.15 1A 16m 14.61 4.30 47.98 1.94e S7 0.10 0.16 1A 24m 15.82 5.35 47.38 3.93b S8 0.09 0.17 1A 24m 16.65 4.73 47.68 3.82b S9 0.09 0.16 1A 24m 16.51 4.75 47.60 5.15a S10 0.11 0.16 1A 24m 15.63 5.91 47.07 2.53c S11 0.00 0.18 1A 28m 18.03 7.31 46.37 1.65f S12 0.00 0.16 1A 28m 16.05 7.58 46.28 1.18h Figure 4. Dendrogram shows relationships of similarity among studied 12 berry species related to all chromosomal characteristics, constructed using Euclidean distance and Single method; cophenet- ic correlation r = 0.75. 140 Saeed Mohammadpour, Ghasem Karimzadeh, Seyed Mahmood Ghaff ari ters. Th e fi rst cluster consists of S11 and S12, the second cluster consists of S7-S10, and the third cluster consists of S2-S6, while S1 is separated within the fourth cluster. Nuclear genome size estimation Th e resultant fl ow cytometric data were fi rst tested for normality test and then analyzed according to a complete- ly randomized design (CRD) with 3 replicate cells. Th e monoploid genome size (2Cx DNA) of these twelve species (Fragaria vesca subsp. vesca, Rubes fruticosus cv. Qaem- shahr, Rubus occidentalis cv. Qaemshahr, Ribes nigrum, Ribes rubrum, Ribes uva-crispa, Lycium barbarum, Lycium infaustum, Lycium ruthenicum, Vaccinium corymbosum, Morus alba cv. Mashhad, Morus rubra cv. Karaj) were 2.1, 1.55, 0.68, 1.86, 1.94, 1.94, 3.93, 3.82, 5.15, 2.53, 1.65, and 1.18 pg, respectively (Table 3). Th e diff erences in nuclear DNA contents among these analyzed species were statisti- cally signifi cant (P < 0.01). Th e highest value was detected in Lycium ruthenicum (S9, 5.15 pg), while Rubus occiden- talis cv. Qaemshahr (S3) was the lowest (0.68 pg), and the mean total value was 2.36 pg (Table 3). Th e histograms obtained for analyzing nuclear DNA contained two peaks (Figure 6). Th e left peaks in S7 to S10 refer to the known samples reference standard, and the right peaks to the unknown samples. In other spe- cies, the peaks on the left refer to unknown specimens, and the peaks on the right refer to the reference stand- ard of known specimens. Young leaves of berries spe- Figure 5. Scatter diagram of 12 berry species: Romero-Zarco asymmetry indices (a), CVCL and MCA (b). Figure 6. Diagram resulted from principal components analysis 1 (highly related to the position of centromere) and 2 (strongly relat- ed to the length of the complements) of studied berry species. Figure 7. Histograms of relative fluorescence intensity were obtained, using WPB isolation buff er aft er simultaneous analysis of nuclei isolated from standard plant and berry species plants. Th e mean channel number and coeffi cient of variation (CV%) value of each peak are also given. 141Karyomorphology, genome size, and variation of antioxidant in twelve berry species from Iran cies contain many cytosolic compounds such as phe- nolic acids and flavonoids, which can interfere with the fluorescent staining of nuclear DNA. To compensate this problem, a new isolation buffer, WPB, was developed to analyze problematic or woody species because the Sodium metabisulfite and PVP-10 in WPB act as reduc- ing and binding agents to prevent the action of interfer- ing phenols and secondary metabolites. This buffer was found to be suitable for the analysis of Lycium species, as evidenced by the CV of DNA peak and flow cytomet- ric histograms of relative fluorescence intensity. Phyloge- netic analysis revealed that Lycium species are included in a monophyletic group, indicating a very close evolu- tionary relationship. This study discovered significant differences in the genome size among these three species (L. barbarum, L. ruthenicum, and L. infaustum). It was also shown in Ribes, Rubus, and Morus species. There is excellent genome size variation between species. It has been proposed that having a large genome has direct physiological consequences for plants, such as earlier flowering time, larger stomata size, and longer life cycles than small genomes. The flowering time may be partly affected by differences in DNA content. Phytochemical studies ANOVA shows significant (P < 0.01) differences between the berries species for the three studied phy- tochemical traits, including total antioxidant capacity (TAC), total phenol content (TPC), and total flavonoid content (TFC). The mean comparisons are shown in Table 4, using LSD test at 0.01 probabilty level. The aver- age total antioxidant capacity (TAC; Table 4) of berries was 32.79 μmol g-1, ranging from 6.92 μmol g-1 (S1; Fra- garia vesca subsp. vesca) to 63.84 μmol g-1 (S3; Rubus occidentalis cv. Qaemshahr). The average total phenol content (TPC; Table 4) of berries was 4.98 μmol g-1, ranging from 3.08 μmol g-1 (S12; Morus rubra cv. Karaj) to 8.61 μmol g-1 (S9; Lycium ruthenicum). The average total flavonoid content (TFC; Table 4) of berries was 5.18 μmol g-1, ranging from 2.43 μmol g-1 (S12; Morus rubra cv. Karaj) to 10.63 μmol g-1 (S9; Lycium ruthenicum). Black raspberry (S3; Rubus occidentalis cv. Qaemshahr) in terms of all three phytochemical traits, had the high- est score among all berry species. Interestingly, the level of antioxidant activity of black raspberry (S3; Rubus occidentalis cv. Qaemshahr) had statistically significant difference among all studied berry species in the present study at a probability level of 1%. For example, it had 6.6-fold increase over strawberry (S1; Fragaria vesca subsp. vesca) and a 24% increase over Black currant (S4; Ribes nigrum). In the next positions, Goose berry (S6; Ribes uva-crispa) and Black currant (S4; Ribes nigrum) had the highest values, respectively. Moreover, strawberry (S1; Fragaria vesca subsp. vesca) and Goji berry (S8: Lycium infaustum) had the lowest one. DISCUSSION Our results provided the basic genetic, genomic and phytochemical information for these species, which are helpful for the construction of genetic and physical maps and whole-genome sequencing in the future and provide scientific insight into the phenolic and antioxidant func- tions to consumers and the nutraceutical industry. For the first time, some of these species’ chromo- some number and karyotype were determined in Iran. Chromosome numbers 2n = 2x = 14, 16, 24, and 28 in these berry species are consistent with other reports. Differences in the karyotype formula of these species, which are geographically, climatically, and habitually different from each other, indicating the existence of chromosomal structural changes in the process of kary- otype evolution and species formation in these species. It varies according to geographical and climatic conditions. Structural changes in the speciation process of these spe- cies can be changed such as displacements, inversions, and other structural changes. According to Stebbins Table 4. Information obtained from phytochemical evaluation of studied 12 berry species. Total Antioxidant Capacity (TAC), Total Flavonoid Content (TFC), Total Phenol Content (TPC). In each column, means with the same symbol letter are not significantly different (P > 0.01), using LSD test. Species TAC (µmol g-1) TFC (mg g-1) TPC (mg g-1) S1 9.62h 2.47f 3.08f S2 38.40c 5.40c 5.75b S3 63.84a 9.60a 7.92a S4 51.57b 6.50b 5.82b S5 24.56e 3.43d 3.73d S6 52.19b 9.23a 7.68a S7 25.06e 3.45d 4.07c S8 16.38g 2.57ef 3.17ef S9 34.95c 10.63a 8.61a S10 27.50de 2.77e 3.28e S11 29.11d 3.70d 3.59d S12 20.32f 2.43f 3.08f Average (range) 32.79 (9.62-63.84) 5.18 (2.43-10.63) 4.98 (3.08-8.61) 142 Saeed Mohammadpour, Ghasem Karimzadeh, Seyed Mahmood Ghaffari (1971), cytogenetic studies can be used to understand better the relationships between species and different populations of a species and to orient the evolutionary tendencies of plants. Each of these plants shows a unique adaptation to the environment in which they grow. As this compromise increases, new varieties or species may appear in plant communities. This study developed an important tool to assess berries’ chemical composition and bioactivity, using different chromatographic meth- ods for its fruits’ comprehensive authentication and quality control. In the present study, all chromosom- al parameters were significantly different. Differences between species for all karyological parameters confirm the intraspecific chromosomal variations. These results are in agreement with those reported by Chen et al. (2013). In the current study, one chromosome type (“m”) formed karyotype formulas of “14m” for S1-S3 (Fragar- ia vesca subsp. vesca, Rubus fruticosus cv. Qaemshahr, Rubus occidentalis cv. Qaemshahr), “16m” for S4-S6 (Ribes nigrum, Ribes rubrum, Ribes uva-crispa), “24m” for S7-S10 (Lycium barbarum, Lycium infaustum, Lycium ruthenicum, Vaccinium corymbosum), and “28m” for S11 and S12 (Morus alba cv. Mashhad, Morus rubra cv. Karaj). Similar to our findings, Chen et al. (2013) report- ed the same one chromosome type in Lycium species. In the present study, the average total length of chromo- somes (CL) in the studied species is 2.83 μm. Karyotype symmetry was more indicative of evolutionary traits in the studied berries species, regardless of ploidy levels. To determine the karyotypic symmetry, several proper- ties were examined in 12 species. Based on TF%, S6 with the highest value (47.98%) has the most symmetric, and S12 with the lowest value (46.28%) has the most asym- metric karyotype. Gradual changes and changes in TF% values may be due to chromosomal abnormalities. Struc- tural changes in chromosome morphology probably due to chromosome duplication or translocation (exchange and displacement) between chromosomes (Das and Teng 1998). According to the CVCL% index, S3 (11.98%) has the most symmetrical karyotype and S11 (18.03%) has the most asymmetric karyotype among other spe- cies. Overall, the coefficient of variation in the sample was low due to the metacentric nature of most chromo- somes, indicating the symmetry and the uniformity of these karyotypes. In general, it can be concluded, species that are more evolutionarily advanced have variations in chromosome type and size and are therefore asymmetri- cally karyotypically. This asymmetry is exacerbated by the displacement of chromosomal fragments. Karyotypic dissimilarity in terms of karyotype symmetry measure- ment parameters can lead to failure in reproduction and inadequate seed production in the offspring while pre- venting successful intraspecific crosses. In other words, the results of such crosses may be somewhat sterile. In the study of Chen et al. (2013) on the genus Lycium, the chromosomal base number for red goji berry (S7; L. barbarum) and black goji berry (S9; L. ruthenicum) was reported to be 12, which is in complete agreement with the findings of the present study (Table 3). Costich et al. (1993) worked on the blueberry plant (S10; Vaccinium corymbosum) showed that the chromosomal base num- ber was 12, which is consistent with the findings of the current study. For the genus Ribes, the base chromosome number 8 was reported by Chiche et al. (2003), which is matched for the three different species of this genus (S4; R. nigrum, S5; Ribes rubrum, S6; R. uva-crispa) in the present study. On the other hand, karyotypes of all species were classified in 1A (Table 3) . For more detailed studies of asymmetry, Romero-Zarco (1986) indices of A1 and A2 were also calculated. The scatter diagram of these indi- ces (Figure 4) presents four species groups. According to the analyzed asymmetry indices, S3 was demonstrated as the most symmetric species, while S2 and S11 were among the most asymmetric species. It has been sug- gested that karyotypes with more asymmetry have a derived status in comparison to those with more sym- metrical morphology (Lakshmi et al. 1984). For exam- ple, differences in chromosome length (CL) may indicate the occurrence of cyclic changes in genome size during the diversification of the genus. Thus, the study of asym- metry indices and variation in genome size is a valuable means for the establishment of the evolutionary relation- ship between the species and the origin of diversification of the populations (Karimzadeh et al. 2011). To determine the total diversity in the popula- tion and the quota of the parameters in the total diver- sity, principal component analysis (PCA) of the karyo- typic parameters was performed. It shows that the first two main components make up 82% of the cumulative changes and they are predicted in a two-dimensional graph (Figure 5). The single dendrogram constructed on the basis of karyotype similarities (Figure 3) shows four major clusters. The first cluster consists of S11 and S12, the second cluster consists of S7-S10, and the third cluster consists of S2-S6, while S1 is separated within the fourth cluster. The arrangement of PCA species from this experiment is fully consistent with the analy- sis obtained by single grouping analysis. Therefore, the results of this study proposed that species within a clus- ter have the most homology in chromosomal variations. For this purpose, a cross between S7 and S8 or S9 is sug- gested because they have the most similarity in their chromosomal characteristics. 143Karyomorphology, genome size, and variation of antioxidant in twelve berry species from Iran The mean monoploid 2Cx DNA value was 2.35 pg in the studied species, verifying intraspecific genome size variation. On the other hand, the monoploid genome size (1Cx DNA) varied from 420.54 Mbp (S3; Rubus occidentalis cv. Qaemshahr) to 2518.35 Mbp (S9; Lycium ruthenicum). According to the results (Chen et al. 2013), red goji berry (S7; Lycium barbarum) and black goji berry (S9; Lycium ruthenicum) had 2Cx DNA values of 4.35 pg and 5.45 pg, respectively. In the pre- sent study, the genome size for the those was 3.93 pg and 5.15 pg, correspondingly, which is almost in agree- ment with the findings of the previous report. The existence of some deviation from this species may be due to the inter-species and intersex diversity. Know- ing genome size may be useful in genome research and studies of relationships between DNA content, physi- ology, and plant ecology (Thiem and Sliwinska 2003). Reports of 2Cx DNA in some berry species have been published, including S3; Rubus occidentalis = 0.60 pg, S5; Ribes rubrum = 1.94 pg, and S6; Ribes uva-crispa = 1.88 pg (Meng and Finn 2002; Chiche et al. 2003). In the present study (Table 3), the monoploid genome size (2Cx DNA) for the latter three species is 0.68 pg, 1.86 pg, and 1.94 pg, respectively, which is in com- plete agreement with the findings of the two previous reports. In the present study, the existence of a statis- tically significant difference in the monoploid genome size indicates the interspecific and intersex diversity among berry species. On the other hand, the correla- tion between 2Cx DNA with environmental conditions (longitude, latitude, and altitude), chromosome length, and the chromosome number showed that 2Cx DNA in the studied berry species had a significant correlation with the chromosome length and chromosome number, and there were no significant correlations between 2Cx DNA with environmental conditions. Hence, it is con- cluded that the genome size of berries species is inde- pendent of changes in environmental conditions. Antioxidant properties estimated by antioxidant assays (DPPH) showed significant differences among different concentrations. Evaluation of total polyphenol content in samples is a widely used method to determine the number of antioxidants in the samples. A rapid, simple, and inexpensive method to measure the antioxi- dant capacity of food involves the use of the free radical 2,2-Diphenyl-1-picrylhydrazyl (DPPH). DPPH radicals are frequently utilized in antioxidant studies. DPPH is widely used to test the ability of compounds to act as free radical scavengers or hydrogen donors and evaluate the antioxidant activity of foods. Antioxidants in a sam- ple can scavenge the DPPH radicals. A gradual reduction in absorbance is observed, indicating that DPPH radicals are being scavenged. Through adding samples, which are rich in antioxidants, to a DPPH solution. Therefore, the percentages we have presented pertain to DPPH radical scavenging capacity, which is directly proportional to antioxidant capacity. The total phenolic content of dry fruits showed a wide range, with values ranging from 3.08 mg GAE/g (S1; Fragaria vesca subsp. vesca) to 8.61 mg GAE/g (S9; Lycium ruthenicum) shown in Table 4. By comparing these data, it can be concluded that the S3 (Rubus occidentalis cv. Qaemshahr) has a high amount of antioxidants. The total flavonoid content (TFC) of dried fruits showed a wide range, with values ranging from 2.43 g QE/g (S12; Morus rubra cv. Karaj) to 10.63 g QE/g (S9; Lycium ruthenicum) as shown in Table 4. Based on the above studies and considering the diversity of these 12 species in terms of cytology and phytochemistry, it is concluded that some of these spe- cies are prone to enter the daily food basket. Black raspberry (S3; Rubus occidentalis cv. Qaemshahr) had the highest value in terms of all three phytochemical traits, followed by goose berry (S6; Ribes uva-crispa) and black currant (S4; Ribes nigrum), respectively. Straw- berry (S1; Fragaria vesca subsp. vesca) and red mulberry (S12; Morus rubra cv. Karaj) had the lowest value. We report that black raspberry (S3; Rubus occidentalis cv. Qaemshahr) is a remarkable source of antioxidant com- pounds compared to other fruits. Hence, research sup- ports deep-colored fruits as potent antioxidant sources. Berries and dried fruit compose a relatively small part of the average diet, but they are important antioxidant sources. Highly pigmented berries have the highest anti- oxidant activity. Such fruits are rich in antioxidant com- pounds that are known for their enhanced stability and bioaccessibility. Based on our findings and the cited lit- eratures, it can be suggested that black raspberries are among the fruits that provide antioxidants. In the study of Islam et al. (2017) on genus Lycium, the mean of total phenol content, total flavonoids content, and antioxidant activity for red goji berries (S7; Lycium barbarum) were measured as 3.16 mg g-1, 2.83 mg g-1, and 16.65 µmol g-1 respectively. These values were 8.33 mg g-1, 11.03 mg g-1, and 34.28 µmol g-1 respectively for black goji berry (S9; Lycium ruthenicum). In the present study (Table 4), these values were 4.07 mg g-1, 3.45 mg g-1, and 25.06 µmol g-1 respectively, for S7 with a significant increase of 29%, 22%, and 50% in all three cases, respectively, compared to the previous report of red goji berry (S7; Lycium barbarum). Also for S9, these values are reported 8.61 mg g-1, 10.63 mg g-1, and 34.95 µmol g-1, respective- ly, which is almost equal to the previous report for (S9; Lycium ruthenicum). In the study of Mustafa Ahmed et al. (2022) on genus Vaccinium, the mean of total phenol 144 Saeed Mohammadpour, Ghasem Karimzadeh, Seyed Mahmood Ghaffari content, total flavonoid content, and antioxidant capacity were reported as 3.73 mg g-1, 2.57 mg g-1, and 17.67 µmol g-1, respectively. In the present study (Table 4), these values were 3.28 mg g-1, 2.77 mg g-1, and 27.50 µmol g-1, respectively which are almost in agreement with those in the previous report. In the current study, the correlation was carried out between phytochemical traits with chromosome num- ber, chromosome length, the monoploid genome size (2Cx DNA), and environmental conditions (latitude and altitude). The results showed no significant correlation between phytochemical traits of species with environ- mental conditions, but there was a significant correlation between phytochemical traits with chromosome num- ber, antioxidant activity with genome size, total phe- nol with chromosome length, and total flavonoid with chromosome length (Table 5). From these results, it is inferred that phytochemical traits in the studied species are independent of changes in environmental conditions. For traits that had a significant correlation between them in the table above, linear regression analysis was performed (Figure 7). Figure 7 shows that the monoploid genome size (2Cx DNA) has a direct linear relationship with chro- mosome length and chromosome number of berry plants (Figures a, and b). Antioxidant has an inverse linear relationship with the genome size and chromo- some number (Figures c, and d). It means that berries with fewer chromosomes and smaller genome sizes pro- duce more antioxidants. Phenol and flavonoid have an inverse linear relationship with the chromosome num- ber (Figures e, and f ), but a direct relationship with the chromosome length (Figures g, and h). For further study, the correlation between genome size and chromosome length for four species (S7, S8, S9, S10) with the same chromosome number (2n = 2x = 24) was also calculated and was significant at the 1% prob- ability level (r = 0.89**). Therefore, Figure 8 shows the relationship of the direct linear regression between the monoploid genome size and chromosome length for the four species mentioned above. In another study, the effect of different ploidy levels on the quantity and quality of essential oils of differ- ent species of berry was investigated and their genetic modification was provided. It is suggested that these berry species should be compared with other species in Iran in terms of cytogenetics and phytochemistry. It is suggested that due to significant differences in genome size and morphology, further research should pin-point rDNA (5S and 45S) sites and additional repetitive DNA elements, using fluorescence in situ hybridization (FISH) to better uncover the processes involved in the chromo- some evolution of these twelve berry species. ACKNOWLEDGMENTS Thanks to Tarbiat Modares University for financial support, and to the Cultivation Development Company of Berries for sharing berry species from the USA. REFERENCES Abedi R, Babaei A, Karimzadeh G. 2015. Karyological and flow cytometric studies of Tulipa (Liliaceae) spe- cies from Iran. Plant Systematics and Evolution. 301: 1473-1484. Ahmadu T, Ahmad K. 2021. An introduction to Bioactive natural products and general applications. In Bioac- tive Natural Products for Pharmaceutical Applica- tions. Springer, Cham.140: 41-91. Amosova AV, Zoshchuk SA, Rodionov AV, Ghukasyan L, Samatadze TE, Punina EO, Loskutov IG, Yurkevich OY, Muravenko OV. 2019. Molecular cytogenetics of valuable Arctic and sub-Arctic pasture grass species from the Aveneae/Poeae tribe complex (Poaceae). BMC genetics. 20(1): 1-16. Barreca D, Laganà G, Leuzzi U, Smeriglio A, Trombetta D, Bellocco E. 2016. Evaluation of the nutraceutical, antioxidant and cytoprotective properties of ripe pis- tachio (Pistacia vera L., variety Bronte) hulls. Food Chemistry. 196: 493-502. Table 5. Correlation coefficients between phytochemical traits of 12 berry species with chromosome number, chromosome length, the mon- oploid genome size (2Cx DNA), and environmental conditions Phytochemical traits Latitude (N) Longitude (E) Altitude (m) Chromosome (No.) CL (µm) 2Cx DNA (pg) Total phenol content (mg g-1) -0.28ns 0.29ns -0.30ns -0.38* 0.40* 0.10ns Total flavonoid content (mg g-1) -0.29ns 0.29ns -0.25ns -0.34* 0.39* 0.09ns Antioxidant capacity (DPPH) (µmol g-1) -0.09ns 0.10ns -0.19ns -0.44* 0.08ns -0.34* ns, *, **: Non significant (P > 0.05), significant differences (P < 0.05) and (P < 0.01), respectively. 145Karyomorphology, genome size, and variation of antioxidant in twelve berry species from Iran Chromosome No. 2C x D N A (p g) Y = 0.577 + 0.088* (± 0.036) X Y = 58.26 - 1.263** (± 0.445) X A nt io xi da nt (µ m ol g- 1 ) Chromosome No.2Cx DNA (pg) A nt io xi da nt (µ m ol g- 1 ) Y = 42.84 - 4.26* (± 2.02) X CL (µm) 2C x D N A (p g) Y = 0.180 + 0.746** (± 0.184) X Ph en ol (m g g- 1 ) Y = 2.53 + 0.839* (± 0.333) X CL (µm) Y = 1.69 + 1.195* (± 0.49) X Fl av on oi d (m g g- 1 ) Y = 8.96 - 0.187* (± 0.088) X Chromosome No. Fl av on oi d (m g g- 1 ) Y = 7.85 - 0.142* (± 0.059) X Chromosome No. Ph en ol (m g g- 1 ) d e g c a f Figure 8. Diagrams of linear regression relationships between diff erent traits of studied berry plants 146 Saeed Mohammadpour, Ghasem Karimzadeh, Seyed Mahmood Ghaff ari Basu A, Rhone M, Lyons TJ. 2010. Berries: emerging impact on cardiovascular health. Nutrition reviews 68(3): 168-177. Basu A. 2019. Role of berry bioactive compounds on lipids and lipoproteins in diabetes and metabolic syn- drome. Nutrients. 11(9): 1983. Bernardello G, Stiefk ens L, Las Penas ML. 2008. Karyo- type studies in Grabowskia and Phrodus (Solanace- ae). Plant Systematic and Evolution. 275: 265-269. Blumberg JB, Camesano TA, Cassidy A, Kris-Etherton P, Howell A, Manach C, Ostertag LM, Sies H, Skulas- Ray A, Vita JA. 2013. Cranberries and their bioactive constituents in human health. Advances in Nutrition. 4(6): 618-632. Briskin DP. 2000. Medicinal plants and phytomedicines. Linking plant biochemistry and physiology to human health. Plant Physiology. 124(2): 507-514. Chen J, Liu X, Zhu L, Wang Y. 2013. Nuclear genome size estimation and karyotype analysis of Lycium species (Solanaceae). Scientia Horticulturae. 151: 46-50. Chiche J, Brown SC, Leclerc JC, Siljak-Yakovlev S. 2003. Genome size, heterochromatin organisation, and ribosomal gene mapping in four species of Ribes. Canadian Journal of Botany. 81(11): 1049-1057. Costich DE, Ortiz R, Meagher TR, Bruederle LP, Vorsa N. 1993. Determination of ploidy level and nuclear DNA content in blueberry by fl ow cytometry. Th eo- retical and Applied Genetics. 86(8): 1001-1006. Das, TK, Teng BS. 1998. Between trust and control: Developing confi dence in partner cooperation in alli- ances. Academy of Management Review. 23(3):491- 512. Doležel J, Bartoš J 2005. Plant DNA fl ow cytometry and estimation of nuclear genome size. Annals of Botany. 95: 99-110. Doležel J, Bartoš J, Voglmayr H, Greilhuber J 2003. Nuclear DNA content and genome size of trout and human. Cytometry. 51: 127-129. Doležel J, Greilhuber J, Suda J. 2007. Flow Cytometry with Plant Cells: Analysis of Genes, Chromosomes and Genomes. Wiley-VCH Verlag GmbH and Co., Weinheim, Germany. Donno D, Beccaro GL, Mellano MG, Cerutti AK, Bou- nous G. 2015. Goji berry fruit (Lycium spp.): anti- oxidant compound fi ngerprint and bioactivity evalua- tion. Journal of Functional Foods. 18: 1070-1085. Dorman HD, Koşar M, Kahlos K, Holm Y, Hiltunen R. 2003. Antioxidant properties and composition of aqueous extracts from Mentha species, hybrids, vari- eties, and cultivars. Journal of Agricultural and Food Chemistry. 51(16): 4563-4569. Faria A, Pestana D, Teixeira D, de Freitas V, Mateus N, Calhau C. 2010. Blueberry anthocyanins and pyruvic acid adducts: anticancer properties in breast cancer cell lines. Phytotherapy Research. 24(12): 1862-1869. Greilhuber J, Doležel J, Lysák MA, Bennett MD. 2005. The origin, evolution and proposed stabilization of the terms ‘genome size’ and ‘C-value’ to describe nuclear DNA contents. Annals of Botany. 95: 255- 260. Greilhuber J. 2007. Cytochemistry and C-values: the less- well-known world of nuclear DNA amounts. Annals of Botany. 101(6): 791-804. Hamidi F, Karimzadeh G, Rashidi Monfared S, Salehi, M. 2018. Assessment of Iranian endemic Artemisia khorassanica: karyological, genome size, and genes expression involved in artemisinin production. Turk- ish Journal of Biology. 42: 322-333. Javadian N, Karimzadeh G, Sharifi M, Moieni A, Behmanesh M. 2017. In vitro polyploidy induction: changes in morphology, podophyllotoxin biosyn- thesis, and expression of the related genes in Linum album (Linaceae). Planta. 245(6): 1165-1178. Kalt W, Cassidy A, Howard LR, Krikorian R, Stull AJ, Tremblay F, Zamora-Ros R. 2020. Recent research on the health benefi ts of blueberries and their anthocya- nins. Advances in Nutrition. 11(2): 224-236. Karimzadeh G, Danesh-Gilevaei M, Aghaalikhani M. 2011. Karyotypic and nuclear DNA variations in Lathyrus sativus (Fabaceae). Caryologia. 64(1): 42-54. Karimzadeh G, Mousavi SH, Jafarkhani-Kermani M, Jala- li-Javaran M. 2010. Karyological and nuclear DNA variation in Iranian endemic muskmelon (Cucumis melo var. Inodorus). Cytologia. 75(4): 451-461. Kim H, Lee GR, Kim J, Baek JY, Jo YJ, Hong SE, Kim HM. 2016. Sulfi redoxin inhibitor induces preferential death of cancer cells through reactive oxygen species- CL (µm) 2C x D N A (p g) Y = 0.573 + 0.983** (± 0.162) X Figure 9. Diagram of the linear regression relationship between genome size and chromosome length for four species (S7-S10) with the same chromosome number (2n = 2x = 24). 147Karyomorphology, genome size, and variation of antioxidant in twelve berry species from Iran mediated mitochondrial damage. Free Radical Biol- ogy and Medicine. 91: 264-274. Laczkó-Zöld E, Komlósi A, Ülkei T, Fogarasi E, Croitoru M, Fülöp I, Domokos E, Ştefănescu R, Varga E. 2018. Extractability of polyphenols from black currant, red currant and gooseberry and their antioxidant activity. Acta Biologica Hungarica. 69(2): 156-169. Lakshmi N, Venkateswara Rao T. 1984. Karyological and morphological investigations on some inbred strain of Trigonella. Genetica-Ibrica. 43: 187-200. Lalegani S, Gavlighi HA, Azizi MH, Sarteshnizi RA. 2018. Inhibitory activity of phenolic-rich pistachio green hull extract-enriched pasta on key type 2 diabetes relevant enzymes and glycemic index. Food Research International. 105: 94-101. Liochev SI. 2013. Reactive oxygen species and the free radical theory of aging. Free Radical Biology and Medicine. 60: 1-4. Loureiro J, Rodriguez E, Doležel J, Santos C. 2007. Two new nuclear isolation buffers for plant DNA flow cytometry: a test with 37 species. Annals of Botany. 100: 875-888. Lu Y, Guo S, Zhang F, Yan H, Qian DW, Wang HQ, Jin L, Duan JA. 2019. Comparison of functional compo- nents and antioxidant activity of Lycium barbarum L. fruits from different regions in China. Molecules. 24(12): 2228. Mahdavi S, Karimzadeh G. 2010. Karyological and nucle- ar DNA content variation in some Iranian endemic Thymus species (Lamiaceae). Journal of Agricultural Science and Technology (JAST). 12(4): 447-458. Meng R, Finn C. 2002. Determining ploidy level and nuclear DNA content in Rubus by flow cytometry. Journal of the American Society for Horticultural Science. 127(5): 767-775. Moure A, Cruz JM, Franco D, Domınguez JM, Sineiro J, Domınguez H, Núñez MJ, Parajó JC. 2001. Natural antioxidants from residual sources. Food Chemistry. 72(2): 145-171. Murashige T, Skoog F. 1962. A revised medium for rapid growth and bio-assays with tobacco tissue cultures. Physiologia Plantarum. 15(3): 437-497. Mustafa AM, Angeloni S, Abouelenein D, Acquaticci L, Xiao J, Sagratini G, Maggi F, Vittori S, Caprioli G. 2022. A new HPLC-MS/MS method for the simulta- neous determination of 36 polyphenols in blueberry, strawberry and their commercial products and deter- mination of antioxidant activity. Food Chemistry. 367: 130743. Nestby R, Hykkerud AL, Martinussen I. 2019. Review of botanical characterization, growth preferences, cli- matic adaptation and human health effects of Eri- caceae and Empetraceae wild dwarf shrub berries in boreal, alpine and arctic areas. Journal of Berry Research. 9(3): 515-547. Nile SH, Park SW. 2014. Edible berries: Bioactive com- ponents and their effect on human health. Nutrition. 30(2): 134-144. Nin S, Petrucci WA, Del Bubba M, Ancillotti C, Giorda- ni E. 2017. Effects of environmental factors on seed germination and seedling establishment in bilberry (Vaccinium myrtillus L.). Scientia Horticulturae. 226: 241-249. Piljac-Žegarac J, Belščak A, Piljac A. 2009. Antioxidant capacity and polyphenolic content of blueberry (Vac- cinium corymbosum L.) leaf infusions. Journal of Medicinal Food. 12(3): 608-614. Potterat O. 2010. Goji (Lycium barbarum and L. Chin- ense): phytochemistry, pharma-cology and safety in the perspective of traditional uses and recent popu- larity. Planta Medica. 76: 7-19. Prior R.L, Cao G, Martin A, Sofic E, McEwen J, O’Brien C, Lischner N, Ehlenfeldt M, Kalt W, Krewer G, Mainland CM. 1998. Antioxidant capacity as influ- enced by total phenolic and anthocyanin content, maturity, and variety of Vaccinium species. Journal of Agricultural and Food Chemistry. 46(7): 2686-2693. Rache L, Pacheco JC. 2010. Propagación in vitro de plan- tas adultas de Vaccinium meridionale (Ericaceae). Acta Botanica Brasilica. 24(4): 1086-1095. Skrovankova S, Sumczynski D, Mlcek J, Jurikova T, Sochor J. 2015. Bioactive compounds and antioxi- dant activity in different types of berries. Interna- tional Journal of Molecular Sciences. 16(10): 24673- 24706. Sliwinska E, Zielinska E, Jedrzejczyk I. 2005. Are seeds suitable for flow cytometric estimation of plant genome size. Cytometry Part A. 64(2): 72-79. Song Q, Sing KC. 2004. Agrobacterium tumefaciens-medi- ated transformation of blueberry (Vaccinium corym- bosum L.). Plant Cell Reports. 23(7): 475-484. Stebbins GL. 1971. Chromosome Evolution in Higher Plants. London: Edward Arnold, UK. Steffen K, Grousset F, Schrader G, Petter F, Suffert M. 2015. Identification of pests and pathogens recorded in Europe with relation to fruit imports. EPPO Bul- letin. 45(2): 223-239. Suda J, Trávníček P. 2006. Reliable DNA ploidy deter- mination in dehydrated tissues of vascular plants by DAPI flow cytometry-new prospects for plant research. Cytometry Part A. 69(4): 273-280. Swift H. 1950. The constancy of desoxyribose nucleic acid in plant nuclei. Proceedings of the National Academy of Sciences of the USA. 36(11): 643-654. 148 Saeed Mohammadpour, Ghasem Karimzadeh, Seyed Mahmood Ghaffari Tarkesh Esfahani S, Karimzadeh G, Naghavi MR. 2016. 2C DNA value of Persian poppy (Papaver bracteatum Lindl.) medicinal plant as revealed by flow cytometry analysis; a quick effective criteria for distinguishing unidentified Papaver species. International Journal of Advanced Biotechnology and Research. 7(2): 573- 578. Tarkesh Esfahani S, Karimzadeh G, Naghavi MR. 2020. In vitro polyploidy induction in Persian poppy (Papaver bracteatum Lindl.). Caryologia. 73(1): 133-144. Tavan M, Mirjalili MH, Karimzadeh G. 2015. In vitro polyploidy induction: changes in morphological, ana- tomical and phytochemical characteristics of Thymus persicus (Lamiaceae). Plant Cell, Tissue and Organ Culture. 122(3): 573-583. Tepe B, Sokmen M, Akpulat HA, Sokmen A. 2006. Screening of the antioxidant potentials of six Salvia species from Turkey. Food Chemistry. 95(2): 200-204. Thiem B, Sliwinska E. 2003. Flow cytometric analy- sis of nuclear DNA content in cloudberry (Rubus chamaemorus L.) in vitro cultures. Plant Science. 164: 126-134. Umdale S, Ahire M, Aiwale V, Jadhav A, Mundada P. 2020. Phytochemical investigation and antioxidant efficacy of wild, underutilized berries of economi- cally important Indian Sandalwood (Santalum album L.). Biocatalysis and Agricultural Biotechnology. 27: 101705. Vinson JA, Su X, Zubik L, Bose P. 2001. Phenol antioxi- dant quantity and quality in foods: fruits. Journal of Agricultural and Food Chemistry. 49(11): 5315-5321. Yang B, Kortesniemi M.. 2015. Clinical evidence on potential health benefits of berries. Current Opinion in Food Science. 2: 36-42. Yang H, Tian T, Wu D, Guo D, Lu J. 2019. Prevention and treatment effects of edible berries for three dead- ly diseases: Cardiovascular disease, cancer and diabe- tes. Critical Reviews in Food Science and Nutrition. 59(12): 1903-1912. Zarabizadeh H, Karimzadeh G, Rashidi Monfared S, Tarkesh Esfahani S. 2022. Karyomorphology, ploidy analysis, and flow cytometric genome size estimation of Medicago monantha populations. Turkish Journal of Botany. 46: 50-61. Zarco CR. 1986. A new method for estimating karyotype asymmetry. Taxon. 35(3): 526-530. Zong Y, Gu L, Shen Z, Kang H, Li Y, Liao F, Xu L, Guo W. 2021. Genome-wide identification and bioinfor- matics analysis of auxin response factor genes in highbush blueberry. Horticulturae. 7(10): 403. Caryologia International Journal of Cytology, Cytosystematics and Cytogenetics Volume 75, Issue 3 - 2022 Firenze University Press