Next Article in Journal
Development and Quality Characteristics of Jangjorim Prepared Using Long-Arm Octopus (Octopus minor) as an Elderly-Friendly Food
Next Article in Special Issue
A Healthy Brazil Nut Beverage with Opuntia stricta var. dillenii Green Extract: Beverage Stability and Changes in Bioactives and Antioxidant Activity during Cold Storage
Previous Article in Journal
Effect of Yogurt Ice Cream on the Viability and Antidiabetic Potential of the Probiotics Lactobacillus acidophilus, Lacticaseibacillus rhamnosus, and Bifidobacterium animalis subsp. lactis after In Vitro Digestion
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Leaf Mustard (Brassica juncea) Germplasm Resources Showed Diverse Characteristics in Agro-Morphological Traits and Glucosinolate Levels

1
National Agrobiodiversity Center, National Institute of Agricultural Sciences, Rural Development Administration, Jeonju 54874, Republic of Korea
2
Department of Biotechnology, Debre Berhan University, Debre Berhan P.O. Box 445, Ethiopia
3
National Creative Research Laboratory for Ca2+ Signaling Network, Jeonbuk National University Medical School, Jeonju 54896, Republic of Korea
4
Agricultural Technology Center of Yeosu, Yeosu 59633, Republic of Korea
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Foods 2023, 12(23), 4374; https://doi.org/10.3390/foods12234374
Submission received: 26 October 2023 / Revised: 23 November 2023 / Accepted: 27 November 2023 / Published: 4 December 2023

Abstract

:
Leaf mustard, characterized by its purple/red/green leaves with a green/white midrib, is known for its thick, tender, and spicy leaves with a unique taste and flavor. There were only a few studies reported on leaf mustard for its morphological and biochemical traits from Korea. A total of 355 leaf mustard accessions stored at the GenBank of the National Agrobiodiversity Center were evaluated for 25 agro-morphological traits and seven intact glucosinolates (GSLs). The accessions showed a wide variation in terms of most of the traits. The quantitative agro-morphological traits varied from 16.0 (leaf length) to 48.7% (petiole width) of the coefficient of variation (CV). The highest variation was observed in glucoiberin (299.5%, CV), while the total GSL showed a CV of 66.1%. Sinigrin, followed by gluconapin and gluconasturtiin, was the most abundant GSL, accounting for as high as 75% of the total GSLs, while glucobrassicanapin and glucoiberin were the least abundant, contributing 0.7% and 0.1% on average, respectively. Sinigrin had a positive significant correlation with all GSLs but gluconasturtiin, while glucobarbarin and gluconasturtiin were highly positively correlated to each other, but least correlated with other GSLs. The leaf length was negatively correlated with sinigrin and glucoiberin. The width of the petiole showed a positive correlation with gluconapin, glucobrassicanapin, and glucobrassicin, while the length of the petiole had a negative correlation with sinigrin, glucobrassicanapin, glucoiberin, glucobrassicin, and the total GSLs. A higher width of the midrib was associated with higher contents of gluconapin, glucobrassicanapin, and glucobrassicin. A PCA analysis based on the agro-morphological traits showed that the first and second principal components accounted for 65.2% of the overall variability. Accessions that form a head tend to exhibit a longer leaf length, a larger plant weight, a thicker midrib, and higher widths of the midrib, petiole, and leaf. The GSLs showed inconsistent inter-and intra-leaf variation. Accessions that identified for various traits in their performance, such as, for example, Yeosu66 and IT259487 (highest total glucosinolates) and IT228984 (highest plant weight), would be promising lines for developing new varieties.

1. Introduction

Mustard (Brassica juncea) is an economically important vegetable globally used as an oil seed, vegetable, and condiment with a long history dating back to 3000 B.C. [1]. The most cultivated mustard species are brown mustard (B. juncea), black mustard (B. nigra), and white/yellow mustard (Sinapis alba) [1,2]. Leaf mustard (B. juncea, 2n = 36, AABB), ‘gat’ in Korean, is a diploid species of the Cruciferous family produced via hybridization and chromosome doubling of B. rapa (2n = 20, AA) and B. nigra (2n = 16, BB) [3]. It is characterized by its vigorously growing leaves of a deep purple/red/green color with a green/white midrib. Known for its high resistance to heat, moisture, and cold stress, depending on variety, it is commonly harvested 35 to 40 days after sowing in the summer and 60 to 70 days in the autumn–winter season. Korean mustard has thick, tender, and spicy leaves with a unique taste, aroma, and flavor; hence, it is not only used as a spice but also as the main ingredient in ‘gat-kimchi’ (leaf mustard kimchi) [4]. Mustard leaves help prolong the storage period of ‘gat-kimchi’ due to their slow fermentation speed and maintain a stable color during storage, which could be related to the glucosinolate degradation products it accumulated, which are known for their antimicrobial and anti-fungal effects [4,5].
Glucosinolates (GSLs) are well-known secondary metabolites found in the whole order of Brassicales. GSLs can be found in the roots, seeds, leaf, and stem of the Brassica plant, usually with the youngest tissues containing the highest amount [6]. Florets contain higher concentrations compared to stalks [7]. The content of glucosinolate in Brassica generally accounts for about 1% of the dry weight in vegetables and can exceed 10% in the seeds of some plants [8,9]. In general, reproductive tissues (florets/flowers/seeds) contain as much as 10 to 40 times higher glucosinolates than vegetative tissues [10]. The predominant glucosinolates include, but are not limited to, sinigrin (allyl-GS) in mustards and horseradish [11], glucoraphanin in broccoli and cabbage [12,13], glucobrassicin in Brussels sprouts [14,15], gluconapin in Chinese kale shoots [16], glucoiberin in cabbage leaves [14,17,18], glucobrassicanapin in pak choi leaves [19] and kimchi cabbage [20], glucomoringin in moringa leaves [8,21,22], glucoraphasatin in radish roots [23], and gluconasturtiin in watercress leaves [24].
Various biochemical and crop improvement research-related reports are available in the literature. One of the major findings of genetic variability studies of B. juncea is delineating the existence of two major gene pools: the Indian and the East European [25]. Pradhan and Pental (2011) [25] have summarized some of the findings before 2007 in genetic divergence using agro-morphology traits, phylogenetic analysis, random amplified polymorphic DNA (RAPD) markers, amplified fragment length polymorphism (AFLP) markers, and biochemicals. In recent years, several studies using agro-morphological characters, SSR molecular markers, and RAPD markers showed a considerable variability among the genotypes of Indian mustard [26,27,28,29,30,31,32,33]. Previous studies of genetic variability among B. juncea germplasm collections based on the levels of glucosinolates were focused on the seeds [34,35,36,37], except for Kim et al. (2016) [38] who studied the levels of desulfo-glucosinolates contents in 210 accessions of Korean leaf mustard (Brassica juncea var. integrifolia). Park et al. (2007) [3] also studied the sinigrin levels of 24 dolsan mustard cultivars and found that the contents varied between 2.07 and 4.38 mg/g. In earlier studies, the aims of breeding programs for Brassica crops were mostly to decrease the level of glucosinolates and increase the oil content [39]. However, in recent years there is a renewed interest in GSLs due to their biocidal and anticarcinogenic effect; hence, increasing the levels of GSLs should also be given attention. Studies on the diversity and relationship among Brassica juncea germplasm based on morphological characters, specifically on leaf mustard from Korea, are also not readily available in the literature. Park et al. (2007) [3] analyzed 21 morphological characters in 24 dolsan leaf mustard cultivars and reported a significant variation in the fresh weight, dry weight, leaf length, and the number of leaves.
This study aimed to examine similarities and differences regarding twenty-five agro-morphological characters and seven glucosinolate levels in the leaves of mustard (Brassica juncea) germplasm collections in the National Agrobiodiversity Center GenBank. We also studied the relationship between the agro-morphological characters and the individual glucosinolates among the leaves of the genetic resources. This study represented one of the largest sets of Brassica juncea samples in a single study and could provide important agro-morphological and biochemical data for further studies. Quantitative data on the glucosinolate contents of Korean mustard leaves using a large population of samples are elusive. Hence, our results provide critical data for mustard consumers and breeders developing new varieties with enhanced levels of bioactive compounds having health-beneficial properties.

2. Materials and Methods

2.1. Plant Materials

The seeds of the leaf mustard plants used in this study were collected from Korea and stored at the Korean GenBank located in the National Agrobiodiversity Center (NAC), Jeonju, Republic of Korea. They were grown in the research farm of NAC in Jeonju (35°49′18″ N 127°08′56″ E). The seeds were sown in plug trays, and the seedlings were grown inside a greenhouse. After 20 days, the healthy-looking seedlings (four to five leaves) were transplanted to an area of 60 × 40 cm per plant in an experimental field. Plant cultural practices were followed as per the recommendation of the Rural Development Administration (RDA) of South Korea. Fertilizers (N-K-P-Ca-B = 65-45-100-100-1.5 kg/10 a) were applied before transplanting the seedlings following RDA standards, and drip irrigation tape was used for watering. A total of 355 accessions were considered in this study. Each accession consisted of 25 plants. Plant growth was maintained using nutrient solutions throughout the growing season. Figure 1 shows some selected representative samples showing the whole plant and the inner, outer, and middle leaves of the Brassica juncea plant.

2.2. Reagents and Standard Chemicals

The chemicals and solvents used for extraction and analysis were of analytical grade and purchased from Fisher Scientific Korea Ltd. (Seoul, Republic of Korea) and Sigma-Aldrich (St. Louis, MO, USA). Standards of sinigrin, gluconapin, glucobrassicanapin, glucoiberin, gluconasturtiin, glucobarbarin, and glucobrassicin were purchased from Phytoplan Diehm & Neuberger GmBH (Heidelberg, Germany). All standards of glucosinolates were of purity greater or equal to 97%.

2.3. Agro-Morphological Characters

Twenty-five agro-morphological characters were investigated based on physical observations using a ruler or a digital scale according to the location, either in the field or in the laboratory. The characters were evaluated based on the International Union for the Protection of New Varieties of Plants (UPOV) descriptors for Brassica juncea L. Czern [40]. Leaf length, leaf width, midrib width, midrib thickness, petiole length, and petiole width were evaluated using five matured leaves of randomly selected different plants of the same accession and the average values were reported. The number of leaves and fresh plant weight are reported as the average of four in five plants of the same accession.

2.4. Sample Pretreatment and UPLC-MS/MS Analysis of Intact GSLs

Ten to fifteen healthy plants were used for the analysis of GSLs. Leaves were collected from the outer, middle, and inner locations of each plant and mixed. In each accession, three replicate samples were prepared. Great care was taken to prevent the thawing of the sample to minimize the enzymatic degradation of the GSLs. The samples were immediately frozen and all the equipment in contact with them was held at subzero temperatures until further processing. The spatial distribution and inter-leaf variation of the GSLs in Brassica juncea were studied using 15 randomly selected genetic resources. The inner, middle, and outer leaves were manually separated for the inter-leaf variation study. Each leaf was then dissected into the top, middle, bottom, and midrib parts, as required. Triplicate samples were prepared for each part.
The samples were processed and extraction was conducted based on a previously reported protocol on Brassica spp. [20], with a slight modification. Briefly, the harvested leaves were placed in vinyl and stored in an −80 °C freezer, followed by lyophilization for 48 h using an LP500 vacuum freeze-drier (Ilshinbiobase Co., Seoul, Republic of Korea). The freeze-dried samples were powdered and extracted as follows: 0.1 g ground sample was mixed with 1 mL of CH3OH/H2O (4:1, v/v) mixture in an e-tube, vortexed (10 min, 30 °C), and centrifuged (13,300 RPM, 4 °C, 10 min) using a VS-180 CFi centrifuge (Vision Scientific Co., Daejeon, Korea). The supernatant was transferred to an LC vial and the GSLs were analyzed immediately using an UPLC-MS/MS.
Qualitative and quantitative analyses of the GSLs were conducted using the ultra-performance liquid chromatography-tandem quadrupole mass spectrometry Acquity UPLC System (Waters, Milford, MA, USA) and the Xevo™ TQ-S system (Waters, MS Technologies, Manchester, UK). Separation was carried out using an Acquity UPLC BEH C18 (2.1 × 100 mm Id, 1.7 μm,) column (Waters Corp., Manchester, UK). The separation conditions were set as follows: a flow rate of 0.5 mL/min; a column temperature of 35 °C; an injection volume of 5 μL; and a sampling chamber kept at 6 °C. The mobile phase compositions were 0.1% of TFA in H2O (A) and 0.1% of TFA in MeOH (B). The elution timetable was set as follows: the initial condition set at 100% of A; 0.0–1.0 min, 100% of A; 1.0–7.0 min, from 100 to 80% of A; 7.0–10 min, from 80 to 0% of A; 10–11 min, from 0 to 100% of A; 11–15 min, 100% of A. The MS/MS analysis was conducted in a negative ion electrospray ionization (ESI-) and multiple reaction monitoring (MRM) mode. The ionization source parameters were set as follows: capillary and con voltages of 3000 V and 54 V, respectively; ion source and desolvation temperatures of 150 °C and 350 °C, respectively; and cone and desolvation gas flow rates of 150 L/h and 650 L/h, respectively. Data acquisition was performed using the MassLynx 4.1 software. The concentrations of individual GSLs were calculated using linear regression equations obtained from calibration curves constructed using authentic standards. The results were calculated from the peak area responses and presented as μmol/kg of the sample dry weight (DW). The MRM transitions and other analysis parameters are presented in Table 1.

2.5. Statistical Analysis

The samples were prepared in biological triplicates and all the data presented are means of the replicate measurements. Pearson correlation and Duncan post hoc tests were performed using SPSS V25.0 (Chicago, IL, USA). The quantitative agro-morphological traits and GSLs contents in Brassica juncea were analyzed using principal component analysis (PCA) with SIMCA v. 13.0.3 (Umetrics, Umea, Sweden) without preprocessing. The data were visualized using score and loading plots of the first and second principal components. Each solid triangle on the score plot and the loading plot represented an individual sample and the contribution of an individual trait to the score, respectively.

3. Results

3.1. Variations in Agro-Morphological Traits among 355 Brassica juncea Accessions

Qualitative agro-morphological characters including plant habit (posture), head formation, leaf attitude at the apical part, leaf type, leaf shape, leaf waxiness, lobation of leaf margins, density of incisions of margins, size of the terminal lobe, number of serrates, leaf anthocyanin coloration, blistering, pubescence, and midrib transection were analyzed for 355 accessions based on physical observations in the experimental field and laboratory. Detailed information about the qualitative morphological characters of the genetic resources is presented in Table 2. The characters were evaluated based on the International Union for the Protection of New Varieties of Plants (UPOV) descriptors for Brassica juncea L. Czern (UPOV, 2017).
Quantitative agro-morphological traits including the petiole (length and width), leaf (length and width), midrib (width and thickness), plant weight, and the number of leaves per plant were measured using a ruler and a digital scale accordingly. Wide variabilities were recorded among the 355 accessions. The leaf length and width ranged from 20.8 to 57.0 cm and from 8.6 to 31.8 cm, with mean values of 36.3 and 18.1 cm, respectively. The length of the petiole ranged between 0.9 and 12.5 cm with an average value of 3.1 cm, while its width ranged between 0.4 and 3.8 cm with an average value of 1.3 cm. The thickness and width of the midrib also exhibited a broad range of 0.27 to 1.05 cm and 0.6 to 3.8 cm, respectively. An individual plant weighed from 91.7 to 1854.3 g of fresh weight, with an average weight of 338.9 g. A single Brassica juncea plant considered in this study could bear anything between 6 to 28 leaves and 16.3 leaves on average. The presence of intra-specific variations is a desirable condition for breeding programs. A wide range of variabilities for different leaf agro-morphological characters in Brassica juncea has been reported in earlier studies [41,42]. Similar observations were also found in Dolsan leaf mustard cultivars from Korea [3]. An analysis of variance showed highly significant variations for all the agro-morphological characters studied, revealing the presence of considerable genetic variation among the genotypes (Table 3). According to the UPOV, the seed color, the density of the incisions of the margin of the leaves, the blistering of the leaves, and plant head formation are agreed as useful grouping characteristics for examining the distinctness of Brassica juncea varieties.

3.2. Identification and Quantification of Individual GSLs in Brassica juncea

Seven GSLs were identified in B. juncea genotypes and quantified using UPLC/MS/MS (Figure 2). A wide range of variations of the GSLs was recorded among the 355 accessions, in which the total glucosinolates content ranged from 320.4 to 8055.3 μmole/kg DW. The aliphatic glucosinolates (SIN, GNA, GBN, and GIB) were predominant throughout the collections, on average representing 67.3% (22.7~98.3%) of the total GSLs, followed by the phenylakyl (GNS and GBB) and indole (GBS) GSLs, contributing 28.9% and 4.0% to the total GSLs, respectively. Among individual glucosinolates, SIN, GNA, and GNS were the most dominant GSLs, with average contents of 635.9, 496.4, and 401.9 μmole/kg DW, respectively. GIB was the least dominant GSLs, representing between 0.0 and 8.2% of the total GSLs throughout the genetic resources investigated.
SIN was found to be the most dominant GSL in 148 (41.7%) of the accessions, and its content ranged from 13.0 to 4184.6 μmole/kg DW, accounting for as high as 75% of the total GSLs. Most of the accessions (296) contained less than 1000 μmole/kg DW, while 57 of them accumulated in the range between 1000 and 4000 μmole/kg DW (Figure 3). Two accessions, IT259503 and IT237840, had 4078.7 and 4184.6 μmole/kg DW of SIN contents, respectively. The SIN content varied highly among the genotypes, with a 92.2% coefficient of variation (Table 3). Another aliphatic glucosinolate that was detected in significant levels was GNA, which was the most dominant in 82 of the accessions. GNA ranged between 44.1 and 7708.7 μmole/kg DW, with an average content of 496.4 μmole/kg DW. It was the most varied dominant GSL, with a 126.6% coefficient of variation. However, about 97% of the resources contained less than 2000 μmole/kg DW, while only 2.8% of the accessions had between 2000 and 5000 μmole/kg DW. Yeosu66 had, outstandingly, the highest GNA content (7708.7 μmole/kg DW). The remaining aliphatic GSLs, GBN and GIB, were found in relatively low amounts and comprised 0.0 to 4.8% (0.7% on average) and 0.0 to 8.2% (0.1% on average) of the total glucosinolates, respectively. Although reports on the contents of intact glucosinolates in the leaves of large sets of Brassica juncea genetic resources are elusive, some studies from seeds showed SIN and GNA as the predominant GSLs [34,37]. Kim et al. (2016) [38] reported nine desulfo-glucosinolates from the leaves of Korean mustard and found SIN and glucoiberverin as the first and second dominant GSLs in 210 accessions.
GNS, which is considered one of the most widely distributed phenylakyl GSL and is highly dominant in watercress leaves [24], was also found abundantly in our samples, dominating in 125 (35.2%) of the accessions, with values ranging between 27.2 and 3393.5 μmole/kg DW. Most of the accessions (~81%) contained 100 to 600 μmole/kg DW of GNS, followed by 24 accessions accumulating between 600 and 1000 μmole/kg DW, and 19 accessions had between 1000 and 3000 μmole/kg DW of GNS. The highest GNS content (3393.5 μmole/kg DW) was recorded in IT236761. Like SIN, GNS also varied highly, with a standard deviation of 342.9 (coefficient of variation 85.3%). GBB was the other phenylakyl GSL detected in our samples in a relatively low amount, and it ranged from 0.7 to 33.9 μmole/kg DW. The only indole glucosinolate detected in our samples was the Trp-derived GBS. The GBS content varied from 17.3 to 351.9 μmole/kg DW, with an average concentration of 57.9 μmole/kg DW. GBS is known for its health-beneficial effects such as its anticancer activity [43].

3.3. Correlation Analysis among GSLs and Agro-Morphological Characters

To investigate the patterns of accumulation of individual glucosinolates and their relationship with quantitative morphological characters, a Pearson correlation analysis was performed using the SPSS software (V25), and the results are presented in Table 4. Among the traits, significant and positive/negative associations were found for some pairs, and some others were uncorrelated. Individual glucosinolates were correlated to each other differently. The highest correlated glucosinolate was SIN, which was significantly (positively) correlated with all but GNS. The two phenylalkyl GSLs, GBB and GNS, were least correlated with all other GSLs but moderately, positively correlated to each other (r = 0.565**). This could be because of the same precursor amino acid (HomoPhe) they both share. As seen in Table 4, the aliphatic GSLs, namely SIN, GNA, GBN, and GIB, were significantly correlated to each other. This could also be attributed to their similarity in their biosynthetic pathway. The leaf length was negatively correlated with SIN and GIB. The length and width of the petiole were correlated with GSLs differently. The width of the petiole showed a positive correlation with GNA, GBN, and GBS, while the length of the petiole had a negative correlation with SIN, GBN, GIB, GBS, and the total GSLs. In addition, a higher width of the midrib was associated with higher contents of GNA, GBN, and GBS. This is in concordance with a previous report where the midrib (white section) accumulates a higher GSLs content than the green part of the leaf of Brassica rapa [20]. Correlation studies of the morphological traits and glucosinolates content in Brassica juncea are elusive in the literature. In a study comprising 99 accessions of Ethiopian mustard (Brassica carinata L.), the SIN content was found to be negatively correlated with most morphological characters investigated, namely, leaf area, leaf length, leaf width, number of primary branches, and plant height, and positively correlated with the length of the petiole [42].

3.4. Relationship between the Morphological Characters and the Levels of GSLs

The relationship between some morphological characters and GSLs were evaluated using Duncan’s test (Table 5). Among the characters studied, plant habit (posture), head formation, leaf type, lobation of leaf margins, density of incisions of the margins, leaf waxiness, pubescence on the lower side of the leaf, midrib transection, and the number of leaf serrates had shown significant variation between accessions in groups of distinct agro-morphological characters either on some individual glucosinolates or on the total GSLs. For example, in terms of plant habit, the horizontal types had significantly higher average total GSLs than the semi-erect mustard leaves. The mustard leaves of the division type were significantly lower in their mean total GSLs content than those of the entire (with broad midrib) type. On the other hand, leaf types such as lyrate vs. division, lyrate vs. entire, division vs. lyrate, and division vs. entire (narrow midrib) were not significantly different from each other in their total GSLs content. The average total GSLs content of entire type mustard leaves that had a broad midrib was the highest compared to any of the leaf types (Table 5). This indicates that the midrib of the mustard leaves accumulated higher levels of GSLs, and it agrees with a previous report on other Brassica species [20]. Regarding the midrib transection, the average contents of SIN, GNA, GBN, GBS, and total GSLs were significantly higher in the horizontal compared to the intermediates and semi-circled types. The lobation of the leaf margins had also affected the levels of total GSLs significantly, where the dissected leaves accumulated fewer GSLs compared to the parted leaves. In general, GBS was found to be the most affected chemical with respect to the morphology of the mustard leaves, which exhibited significant variations in eight of the characteristics.

3.5. Principal Component Analysis (PCA) Based on Some Agro-Morphological Traits

Unsupervised PCA analyses were conducted to identify the prevalence of clusters of genotypes and the most relevant agro-morphological traits. The PCA was carried out based on eight quantitative agro-morphological characters. The first three principal components with Eigenvalues of 3.7, 1.43, and 0.99 accounted for 47.3, 17.9, and 12.5% of the overall variability, respectively. All the traits, but the petiole length (−0.0555), contributed positively to PC1. On the other hand, petiole width (−0.3438), midrib width (−0.2427), and leaf width (−0.0884) contributed negatively to PC2, whereas all the other characters had positive contributions, as shown in Table 6. For PC3, the highest positive and negative contributions were due to the number of leaves (0.9022) and the length of the petiole (−0.2568), respectively. As seen in Figure 4a–e, some genetic resources tend to associate based on specific agro-morphological characteristics such as plant habit, head formation, leaf shape, and leaf type. For example, in terms of the plant growth habit, most erect mustard accessions were distributed on the positive side of PC1, while the horizontal types were on the negative side of PC1 (Figure 4b). The erect types were characterized by a higher petiole width, higher leaf and midrib width, a thicker midrib, a longer leaf length, and a larger plant weight. However, mustard leaves with longer petioles are related to semi-erect and horizontal plant habits (posture). Most of the mustard accessions with longer leaves, wider midrib, leaf, and petiole, and a thicker midrib tended to group at the positive side of the PC1 and were characterized by the entire leaf type as opposed to the lyrate and division types (Figure 4c). The spatulate and lanceolate mustard leaf shapes were characterized by shorter leaf width, midrib width, and petiole width, a thinner midrib, a shorter leaf length, as well as a smaller plant weight. Whereas, circular (compressed), broadly elliptic, and elliptic leaves are characterized by higher widths of the leaf, midrib, and petiole, a thicker midrib, a longer leaf length, and a larger plant weight (Figure 4d). Mustard accessions that form head tend to exhibit a longer leaf length, a larger plant weight, a thicker midrib, and higher widths of the midrib, petiole, and leaf (Figure 4e).

3.6. Intra- and Inter-Leaf Distribution of GSLs in Brassica juncea

Fifteen randomly selected light green, green, light purple, and dark purple colored accessions were used to study the spatial distribution of glucosinolates in mustard leaves. To study the distribution of GSLs based on their location/development, the leaves were sampled based on their location (development) as inner (young), middle, and outer (old). On the other hand, to study the spatial distribution pattern within the leaf, the outer part (older) leaves were separated into the top, middle, bottom, and midrib sections. The distribution of total glucosinolate levels based on the leaf sections and position in the plant are presented in Figure 5. There was a significant variation in the GSLs levels between different layers (inner, middle, outer) of leaf samples within the same plant and between the sections (top, middle, and bottom) within the outer leaves. The total GSL level was significantly higher in the outer layers of the leaves (Figure 5a) in 53% of the samples, while the inner and middle layers had significantly higher GSLs in 20% of the samples each. Contrary to this report, the younger leaves (inner leaves) of B. oleracea var. capitate and Raphanus sativus were shown to exhibit higher GSLs levels compared to the older leaves (outer leaves) [44,45]. On the other hand, the top section of the leaves had significantly high GSLs in 47% of the samples, whereas approximately 40% of the samples had significantly higher glucosinolates in their middle and bottom leaf sections (Figure 5b). However, the trend in the total GSLs content based on neither the leaf layers (inner, middle, outer) nor the sections (top, middle, and bottom) of the outer leaves was strictly consistent. A similar observation has been observed in a study employing three samples of kimchi cabbage [20].

4. Discussion

In this study, we employed 355 Brassica juncea accessions stored at the GenBank of the National Agrobiodiversity Center to study the variations on 25 agro-morphological characters and seven intact glucosinolates. The accessions showed wide variations in terms of both their agro-morphological traits and their GSLs levels. Significant correlations between the traits were also observed. We also performed an unsupervised PCA to identify the prevalence of clusters of accessions based on some selected agro-morphological traits. Finally, based on their quantitative trait performance, we identified potential candidate accessions for further studies.
The Brassica juncea accessions showed a huge variation regarding their agro-morphological characters. A different range of variability within and between the leaves has been observed for various agro-morphological characters including plant weight, number of leaf lobes, leaf length, leaf width, number of leaves per plant, anthocyanin coloration, petiole length, leaf shape, and number of serrates [3,29,42,46,47,48,49]. The variability in the agro-morphological characteristics observed in our study is in concordance with some of the previous observations. In addition to the agro-morphological traits of the leaves, previous reports had also explored the diversity in other parts of mustard plants [29,47,48,49], in their physiological traits, and in their biochemical traits such as their antioxidant activities and the total soluble sugars [46]. Rabbani et al. (1998) [49] investigated various morphological traits in the seedling, flowering, and maturity stages of 52 accessions of Brassica juncea collected from Pakistan and found narrow phenotypic variations amongst them. The authors highlighted that the reason could be the narrow genetic base of the germplasms that have undergone a high level of genetic erosion. In another multivariate agro-morphological study using the leaf, flower, siliqua, stem, seed, and phenological characters of Indian mustard, 62 accessions showed four distinct groups of varieties based on geographical location, with leaf, stem, and phenological characteristics accounting for much of the variability [48]. The diverse characteristics in the agro-morphology along with the biochemical contents of Brassica juncea plants offer opportunities for consumption purposes and the development of improved varieties.
Seven individual glucosinolates were detected in the leaves of mustard accessions. Sinigrin, followed by gluconapin and gluconasturtiin, was the most abundant glucosinolate in most of the accessions. Our results were consistent with previous studies on germplasm collections of leaves [38] and seeds [37] of Brassica juncea. The number of glucosinolates detected in our study was fewer compared to the reports of Kim et al. (2016) [38] on Brassica juncea leaves, with glucocheirolin and glucoiberverin not being detected in our study. This could be attributed to the difference in the method of extraction and the instrument used. These authors performed desulphonation during extraction and determined the desulphated GSLs, while, in this study, the glucosinolates were determined as intact glucosinolates. The indole glucosinolates, 4-hydroxyglucobrassicin, neoglucobrassicin, and 4-methoxyglucobrassicin, that have been reported in the seeds and roots of mustard [38,50] were also undetected in our study. Glucoraphasatin, the main glucosinolate of Raphanus sativus, has been reported in the roots of red mustard [50] but was undetected in our study. This could be due to the difference in the type of plant part used. The HomoMet-derived aliphatic glucosinolate, sinigrin, is highly abundant in Brassica juncea seeds with a concentration as high as 109.9 µmole/g DW and is showed to suppress nematode activity [11]. One of the major contributors to the total glucosinolates level in our samples, gluconapin, is also highly abundant in Chinese kale shoots (up to138.6 μmole/g DW) [16], pak choi (up to 70.67 μmole/g DW) [19], and leaf blade of Ezo-wasabi (up to 168.4 μmole/g DW) [51]. The overall glucosinolates levels in this study were lower compared to the previous reports discussed above. This could be related to a difference in sample preparation and recovery protocols. As observed previously in genotypes of mustard samples from Ethiopia [42,52], Korea [38], and different parts of the world [34,37], the glucosinolate concentrations showed a wide intra-species variation in the current study. The diversity of genetic resources reflects the frequency of important phytochemical, phenological, and morphological traits in germplasm collection, and breeding programs largely exploit the variation of these traits. The variation in the GSL content could determine important roles in the physiology, productivity, nutrition, and health benefits of mustard plants.
GSLs showed inconsistent inter- and intra-leaf variation. The spatial and stage-wise distribution of GSLs in the leaves of Brassica juncea could be used to explore how plants use GSLs in the defense management, ecological significance, and biosynthesis mechanisms of GSLs in Brassica plants. To the best of our knowledge, such studies are elusive in the literature. Previously, the levels of glucosinolates were found to be affected by leaf position [44] and leaf cross section [45]. Choi et al. (2014) [44] compared the level of glucosinolates between the inner and outer leaves of Brassica oleracea var. capitata in spring and fall sawing and found that the inner leaves contained higher individual glucosinolates in most of the samples. A study about the leaf spatial patterns of glucosinolate levels among different leaf regions of Raphanus sativus showed that proximal leaflets had significantly more glucosinolates compared to the main leaflet, edge, vein, and middle [45]. In a recent intra- and inter-leaf variation study on three kimchi cabbage varieties (Brassica rapa), Rhee et al. (2020) [20] reported that the proximal half of the leaves contained higher GSLs in most cases, that the middle layers had higher average total GSL levels compared to the inner and outer layer leaves, and suggested that the result was not strictly consistent. This experiment also found an inconsistent trend in the distribution of GSLs based on the leaf layers as well as the leaf sections.
This study is entrusted to identify potential candidates of genetic resources for different important traits for producing GSLs-dense lines and mustard leaves with high yields and other important leaf agro-morphological characters. Table 7 summarizes the promising accessions for various morphological and biochemical traits. Breeding strategies for improved glucosinolate/isothiocyanate profiles of Brassicaceae species are dependent on various factors, such as consumer preference, postharvest processing, biotic and abiotic responses, elicitor responses, and species selection [53]. The usual trend in many horticultural crops’ breeding has been mainly focused on increasing the yield, compromising the biochemical contents. Producing nutrient-dense Brassica juncea crops in terms of GSLs and adding them to our diet could help reduce the onset of chronic diseases such as cancer [5,54,55]. In addition, some GSLs could also contribute to producing varieties resistant to biotic and abiotic stress. Breeding crops that contain agro-ecologically relevant metabolic profiles requires investigation of the chemical diversity of large collections of germplasm. GSLs are known for their potential health benefits to humans [43,56] and as defense compounds for disease resistance in plants [57,58]. Producing lines of Brassica juncea with enhanced levels of GSLs requires assessment of diverse genetic resources in terms of biochemical components as well as agronomical traits. In this study, we have found that some of the genotypes highly accumulated specific GSLs. For example, the highest SIN content was recorded in IT259503 (4078.7 μmole/kg DW) and IT237840 (4184.6 μmole/kg DW), while the lowest SIN (13.0 μmole/kg DW) and the highest GNA (7708.6 μmole/kg DW) were in ‘Yeosu66’. The highest amounts of GNS and GBS were detected in IT236761 (3393.5 μmole/kg DW) and IT248036 (351.9 μmole/kg DW), respectively. The highest GIB content (93.6 μmole/kg DW), 3-fold higher than the second-highest in the entire domain, was detected in IT102894. In general, two accessions, ‘Yeosu66’ and IT259487, had relatively higher total glucosinolates (>6000 μmole/kg DW), while ‘Yeosu66’, IT237840, IT250121, IT259487, IT259503, and ‘Yeosu17’ all exhibited greater than 5000 μmole/kg DW total aliphatic glucosinolates. Given that each of the GSLs has its own significance for human health and in the plant defense mechanism, the use of these accessions in producing GSL-rich mustard lines would be highly beneficial.
Some accessions showed a superior performance on some of the agro-morphological traits. For example, IT228984 had the highest plant weight (1854.3 g), which was 1.69- and 5.47-fold higher than the second-highest and average plant weights in the whole domain. IT120115 and IT118972 were superior in leaf length, with 57 cm and 55.2 cm, respectively. Similarly, IT2288223, Yeosu42, Yeosu61, and Yeosu84 exhibited leaf widths greater than 30 cm and could be the best resources for producing wider-leaf mustards. As observed in the correlations analysis (Table 4), a wider petiole was correlated with higher levels of GNA, GBN, and GBS, while the tallness of the petiole was associated with higher levels of most GSLs. IT228223 and Yeosu80, which were also characterized by large leaf and midrib widths, could be the best germplasms in terms of width of the petiole. The superior genotypes identified for different traits in this study should also be investigated for their agro-climatic preference, disease resistance capability, and other biochemical characteristics, including their proximate compositions, antioxidant activity, and secondary metabolites, in order to be utilized in breeding programs for developing lines with enhanced yields and health benefits.

5. Conclusions

A total of 355 accessions collected from Korea and stored in the National Agrobiodiversity Center GenBank were evaluated for intact GSLs contents and agro-morphological characters. Four aliphatic glucosinolates (SIN, GNA, GBN, and GIB), two phenylalkyl glucosinolates, and one indole GSL were detected and determined using UPLC-MS/MS. Seventeen qualitative and eight quantitative agro-morphological characters were also evaluated. A wide variability was observed for most of the traits. This study provides enriching information for the scientific community and consumers about the agro-morphological traits and glucosinolate diversity in B. juncea and offers guidance for B. juncea selection and the development of high-quality cultivars. Promising potential candidate accessions in terms of their performance in various agro-morphological traits and GSLs levels are also identified. This report provides critical data specifically for breeders developing new varieties with enhanced levels of glucosinolates. Considering the anticancer properties of glucosinolates reported elsewhere, this report would be of interest to consumers, nutraceutical companies engaged in formulating anti-cancerous food supplements, and drug developers. Further investigation using molecular markers, disease resistance, and other agro-morphological and biochemical traits are also apt to develop high-quality leaf mustard cultivars. To the best of our knowledge, the number of accessions we used represented one of the largest sets in a single study to date.

Author Contributions

Conceptualization, A.D.A., S.-H.K. and J.-H.R.; methodology, A.D.A.; software, A.D.A.; validation, S.-H.K., B.-S.H. and H.C.K.; formal analysis, A.D.A. and Y.-J.C.; investigation, A.D.A., N.R. and P.S.; resources, Y.-H.L.; data curation, A.D.A. and J.-H.R.; writing—original draft preparation, A.D.A.; writing—review and editing, A.D.A., S.-H.K. and J.-H.R.; visualization, B.-S.H.; supervision, J.-H.R.; project administration, S.-H.K.; funding acquisition, S.-H.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was carried out with the support of the “Research Program for Agricultural Science & Technology Development (Project No. PJ01423001, RS-2019-RD008922)” of the National Institute of Agricultural Sciences, Rural Development Administration, Republic of Korea.

Data Availability Statement

The datasets generated for this study are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Thomas, J.; Kuruvilla, K.M.; Hrideek, T.K. Mustard. In Handbook of Herbs and Spices, 2nd ed.; Peter, K.V., Ed.; Woodhead Publishing Limited: Sawston, UK, 2012; Volume 1, pp. 388–398. [Google Scholar]
  2. Lietzow, J. Biologically active compounds in mustard seeds: A toxicological perspective. Foods 2021, 10, 2089. [Google Scholar] [CrossRef] [PubMed]
  3. Park, H.J.; Lee, I.H.; Kim, H.T.; Yang, S.Y.; Nou, I.S. Analysis of genetic similarity of dolsan leaf mustard based on several horticultural traits and RAPD. Korean J. Hortic. Sci. Technol. 2007, 25, 305–310. [Google Scholar]
  4. Kang, S.K. Antimicrobial activities in the Korean traditional leaf mustard, Brassica juncea Coss. Korean J. Plant Resour. 2005, 8, 96–102. [Google Scholar]
  5. Vig, A.P.; Rampal, G.; Thind, T.S.; Arora, S. 2009. Bio-protective effects of glucosinolates—A review. LWT-Food Sci. Technol. 2009, 42, 1561–1572. [Google Scholar] [CrossRef]
  6. Blažević, I.; Mastelić, J. Glucosinolate degradation products and other bound and free volatiles in the leaves and roots of radish (Raphanus sativus L.). Food Chem. 2009, 113, 96–102. [Google Scholar] [CrossRef]
  7. Kushad, M.M.; Brown, A.F.; Kurilich, A.C.; Juvik, J.A.; Klein, B.P.; Wallig, M.A.; Jeffery, E.H. Variation of glucosinolates in vegetable crops of Brassica oleracea. J. Agric. Food Chem. 1999, 47, 1541–1548. [Google Scholar] [CrossRef]
  8. Bennett, R.N.; Mellon, F.A.; Foidl, N.; Pratt, J.H.; Dupont, M.S.; Perkins, L.; Kroon, P.A. Profiling glucosinolates and phenolics in vegetative and reproductive tissues of the multi-purpose trees Moringa oleifera L. (Horseradish tree) and Moringa stenopetala L. J. Agric. Food Chem. 2003, 51, 3546–3553. [Google Scholar] [CrossRef]
  9. Bennet, R.N.; Mellon, F.A.; Kroon, P.A. Screening Crucifer seeds as sources of specific intact glucosinolates using ion-pair high-performance liquid chromatography negative ion electrospray mass spectrometry. J. Agric. Food Chem. 2004, 52, 428–438. [Google Scholar] [CrossRef]
  10. Clarke, D.B. Glucosinolates, structures and analysis in food. Anal. Methods 2010, 2, 310–325. [Google Scholar] [CrossRef]
  11. Zasada, I.A.; Ferris, H. Nematode suppression with brassicaceous amendments: Application based upon glucosinolate profiles. Soil Biol. Biochem. 2004, 36, 1017–1024. [Google Scholar] [CrossRef]
  12. Jones, R.B.; Frisina, C.L.; Winkler, S.; Imsic, M.; Tomkins, R.B. Cooking method significantly effects glucosinolate content and sulforaphane production in broccoli florets. Food Chem. 2010, 123, 237–242. [Google Scholar] [CrossRef]
  13. Pérez-Balibrea, S.; Moreno, D.A.; García-Viguera, C. Glucosinolates in broccoli sprouts (Brassica oleracea var. italica) as conditioned by sulphate supply during germination. J. Food Sci. 2010, 75, 673–677. [Google Scholar]
  14. Charron, C.S.; Saxton, A.M.; Sams, C.E. Relationship of climate and genotype to seasonal variation in the glucosinolate-myrosinase system. I. Glucosinolate content in ten cultivars of Brassica oleracea grown in fall and spring seasons. J. Sci. Food Agric. 2005, 85, 671–681. [Google Scholar] [CrossRef]
  15. Carlson, D.G.; Daxenbichler, M.E.; Vanetten, C.H.; Kwolek, W.F.; Williams, P.H. Glucosinolates in Crucifer vegetables: Broccoli, Brussels sprouts, cauliflower, collards, kale, mustard greens, and kohlrabi. Am. Soc. Hortic. Sci. 1987, 112, 173–178. [Google Scholar] [CrossRef]
  16. Qian, H.; Liu, T.; Deng, M.; Miao, H.; Cai, C.; Shen, W.; Wang, Q. Effects of light quality on main health-promoting compounds and antioxidant capacity of Chinese kale sprouts. Food Chem. 2016, 196, 1232–1238. [Google Scholar] [CrossRef]
  17. Tiedink, H.G.M.; Davies, J.A.R.; van Broekhoven, L.W.; Van der Kamp, H.J.; Jongen, W.M.F. Formation of mutagenic N-nitroso compounds in vegetable extracts upon nitrite treatment: A comparison with the glucosinolate content. Food Chem. Toxicol. 1988, 26, 947–954. [Google Scholar] [CrossRef]
  18. Verkerk, R.; Dekker, M. Glucosinolates and myrosinase activity in red cabbage (Brassica oleracea L. var. Capatita f. rubra DC.) after various microwave treatments. J. Agric. Food Chem. 2004, 52, 7318–7323. [Google Scholar]
  19. Wiesner, M.; Zrenner, R.; Krumbein, A.; Glatt, H.; Schreiner, M. Genotypic variation of the glucosinolate profile in pak choi (Brassica rapa ssp. chinensis). J. Agric. Food Chem. 2013, 61, 1943–1953. [Google Scholar] [CrossRef]
  20. Rhee, J.H.; Choi, S.; Lee, J.E.; Hur, O.S.; Ro, N.Y.; Hwang, A.J.; Ko, H.; Chung, Y.; Noh, J.; Assefa, A.D. Glucosinolate content in Brassica genetic resources and their distribution pattern within and between inner, middle, and outer leaves. Plants 2020, 9, 1421. [Google Scholar] [CrossRef]
  21. Doerr, B.; Wade, K.L.; Stephenson, K.K.; Reed, S.B.; Fahey, J.W. Cultivar effect on Moringa oleifera glucosinolate content and taste: A pilot study. Ecol. Food Nutr. 2009, 48, 199–211. [Google Scholar] [CrossRef]
  22. Förster, N.; Ulrichs, C.; Schreiner, M.; Müller, C.T.; Mewis, I. Development of a reliable extraction and quantification method for glucosinolates in Moringa oleifera. Food Chem. 2015, 166, 456–464. [Google Scholar] [CrossRef] [PubMed]
  23. Han, N.; Súudi, M.; Kim, J. The major aliphatic glucosinolate content in Korean radish during vegetative and reproductive growth. Hortic. Environ. Biotechnol. 2015, 56, 152–158. [Google Scholar] [CrossRef]
  24. Giallourou, N.; Oruna-Concha, M.J.; Harbourne, N. Effects of domestic processing methods on the phytochemical content of watercress (Nasturtium officinale). Food Chem. 2016, 212, 411–419. [Google Scholar] [CrossRef] [PubMed]
  25. Pradhan, A.K.; Pental, D. Genetics of Brassica juncea. In Genetics and Genomics of the Brassicaceae; Plant Genetics and Genomics: Crops and Models; Schmidt, R., Bancroft, I., Eds.; Springer: New York, NY, USA, 2011; Volume 9, pp. 323–345. [Google Scholar]
  26. Devi, T.R.; Devshini Devi, N.; Vivekananda, Y.; Sharma, P.R. Genetic diversity analysis in Indian mustard (Brassica juncea L. Czern and Coss) genotypes using agro-morphological parameters. Electron. J. Plant Breed. 2017, 8, 749–753. [Google Scholar] [CrossRef]
  27. Khan, M.A.; Rabbani, M.A.; Munir, M.; Ajmal, S.K.; Malik, M.A. Assessment of genetic variation within Indian mustard (Brassica juncea) germplasm using random amplified polymorphic DNA markers. J. Integr. Plant Biol. 2008, 50, 385–392. [Google Scholar] [CrossRef] [PubMed]
  28. Khan, W.M.; Munir, I.; Arif, M.; Iqbal, A.; Ali, I.; Ahmad, D.; Ahmad, M.; Mian, A.; Bakht, J.; Swati, Z.A. Assessment of genetic diversity of Brassica juncea germplasm using randomly amplified polymorphic DNA (RAPD) markers. Afr. J. Biotechnol. 2011, 10, 3654–3658. [Google Scholar]
  29. Saleem, N.; Ahmad Jan, S.; Jawaad Atif, M.; Khurshid, H.; Ali Khan, S.; Abdullah, M.; Jahanzaib, M.; Ahmed, H.; Ullah, S.F.; Iqbal, A.; et al. Multivariate based variability within diverse Indian mustard (Brassica juncea L.) genotypes. Open J. Genet. 2017, 07, 69–83. [Google Scholar] [CrossRef]
  30. Sharma, H.K.; Kumar, A.; Singh, V.V.; Meena, H.S.; Priyamedha; Sharma, P.; Rai, P.K. Variability and genetic diversity study based on agro- morphological traits in a diverse set of Indian mustard [Brassica juncea (L.) Czern. & Coss.] germplasm. J. Environ. Biol. 2021, 42, 22438. [Google Scholar]
  31. Shekhawat, N.; Jadeja, G.C.; Singh, J. Genetic variability for yield and its components in Indian mustard (Brassica juncea L. Czern & Coss). Electron. J. Plant Breed. 2014, 5, 117–119. [Google Scholar]
  32. Tiwari, V.K. Morphological parameters in breeding for higher seed yield in Indian mustard [Brassica juncea (L.) Czern. & Coss.]. Electron. J. Plant Breed. 2019, 10, 187–195. [Google Scholar]
  33. Verma, K.; Tripathi, M.K.; Tiwari, S.; Tripathi, N. Analysis of genetic diversity among Brassica juncea genotypes using morpho-physiological and SSR Markers. Int. J. Curr. Microbiol. Appl. Sci. 2021, 10, 1108–1117. [Google Scholar] [CrossRef]
  34. Merah, O. Genetic variability in glucosinolates in seed of Brassica juncea: Interest in mustard condiment. J. Chem. 2015, 2015, 606142. [Google Scholar] [CrossRef]
  35. Saikia, S.L.; Rai, G.K.; Salgotra, R.; Rai, S.; Singh, M.; Rai, P. Erucic acid and glucosinolate variability in Brassica juncea L. Int. J. Chem. Stud. 2018, 6, 1223–1226. [Google Scholar]
  36. Velasco, L.; Becker, H.C. Variability for seed glucosinolates in a germplasm collection of the genus Brassica. Genet. Resour. Crop Evol. 2000, 47, 231–238. [Google Scholar] [CrossRef]
  37. Zhang, C.; Di, H.; Lin, P.; Wang, Y.; Li, Z.; Lai, Y.; Li, H.; Sun, B.; Zhang, F. Genotypic variation of glucosinolates and their breakdown products in mustard (Brassica juncea) seeds. Sci. Hortic. 2022, 294, 110765. [Google Scholar] [CrossRef]
  38. Kim, H.W.; Ko, H.C.; Baek, H.J.; Cho, S.M.; Jang, H.H.; Lee, Y.M.; Kim, J.B. Identification and quantification of glucosinolates in Korean leaf mustard germplasm (Brassica juncea var. integrifolia) by liquid chromatography–electrospray ionization/tandem mass spectrometry. Eur. Food Res. Technol. 2016, 242, 1479–1484. [Google Scholar]
  39. Sodhi, Y.S.; Mukhopadhyay, A.; Arumugam, N.; Verma, J.K.; Gupta, V.; Pental, D.; Pradhan, A.K. Genetic analysis of total glucosinolate in crosses involving a high glucosinolate Indian variety and a low glucosinolate line of Brassica juncea. Plant Breed. 2002, 121, 508–511. [Google Scholar] [CrossRef]
  40. International Union for the Protection of New Varieties of Plants (UPOV). Guidelines for the Conduct of Tests for Distinctness, Uniformity and Stability of Brassica juncea L. Czern. 2017. Available online: http://www.upov.int (accessed on 30 March 2022).
  41. Yadav, R.; Singh, R.; Kumar, S.; Prasad, T.V.; Bharadwaj, R.; Kaur, V.; Petapadi, A.; Kumar, A. Genetic diversity among indigenous germplasm of Brassica juncea (L.) Czern and Coss, using agro-morphological and phenological traits. Proc. Natl. Acad. Sci. India Sect. B Biol. Sci. 2017, 87, 1125–1131. [Google Scholar] [CrossRef]
  42. Teklehaymanot, T.; Wang, H.; Liang, J.; Wu, J.; Lin, R.; Zhou, Z.; Cai, X.; Wang, X. Variation in plant morphology and sinigrin content in Ethiopian mustard (Brassica carinata L.). Hortic. Plant J. 2019, 5, 205–212. [Google Scholar] [CrossRef]
  43. Possenti, M.; Baima, S.; Raffo, A.; Durazzo, A.; Giusti, A.M.; Natella, F. Glucosinolates in food. In Glucosinolates, Reference Series in Phytochemistry; Mérillon, J.M., Ramawat, K.G., Eds.; Springer: Cham, Switzerland, 2017; pp. 87–132. [Google Scholar]
  44. Choi, S.H.; Park, S.; Lim, Y.P.; Kim, S.J.; Park, J.T.; An, G. Metabolite profiles of glucosinolates in cabbage varieties (Brassica oleracea var. capitata) by season, color, and tissue position. Hortic. Environ. Biotechnol. 2014, 55, 237–247. [Google Scholar] [CrossRef]
  45. Shelton, A.L. Within-plant variation in glucosinolate concentrations of Raphanus sativus across multiple scales. J. Chem. Ecol. 2005, 31, 1711–1732. [Google Scholar] [CrossRef] [PubMed]
  46. Sharma, D.; Lal, J.N.; Singh, S.S.P.; Parmar, N. Genetic diversity in leafy mustard (Brassica juncea var. rugosa) as revealed by agro-morphological traits and SSR markers. Physiol. Mol. Biol. Plants 2020, 26, 2005–2018. [Google Scholar] [PubMed]
  47. Singh, K.H.; Shakya, R.; Nanjundan, J.; Thakur, A.; Singh, K.; Singh, K.K. Genetic diversity assessment and characterization of Indian mustard (Brassica juncea L.) varieties using agro-morphological traits. Indian J. Plant Genet. Resour. 2018, 31, 44–50. [Google Scholar] [CrossRef]
  48. Singh, K.H.; Shakya, R.; Mahawar, R.K. Genetic diversity and patterns of variation among Indian mustard (Brassica juncea (L.) Czernj. & Cosson) varieties. SABRAO J. Breed. Genet. 2014, 46, 329–339. [Google Scholar]
  49. Rabbani, M.A.; Iwabuchi, A.; Murakami, Y.; Suzuki, T.; Takayanagi, K. Phenotypic variation and the relationships among mustard (Brassica juncea L.) germplasm from Pakistan. Euphytica 1998, 101, 357–366. [Google Scholar] [CrossRef]
  50. Cuong, D.M.; Kim, J.K.; Bong, S.J.; Baek, S.A.; Jeon, J.; Park, J.S.; Park, S.U. Comparative analysis of glucosinolates and metabolite profiling of green and red mustard (Brassica juncea) hairy roots. 3 Biotech 2018, 8, 382. [Google Scholar] [CrossRef]
  51. Abe, K.; Kido, S.; Maeda, T.; Kami, D.; Matsuura, H.; Shimura, H.; Suzuki, T. Glucosinolate profiles in Cardamine fauriei and effect of light quality on glucosinolate concentration. Sci. Hortic. 2015, 189, 12–16. [Google Scholar] [CrossRef]
  52. Warwick, S.I.; Gugel, R.K.; Mcdonald, T.; Falk, K.C. Genetic variation of Ethiopian mustard (Brassica carinata A. Braun) germplasm in western Canada. Genet. Resour. Crop Evol. 2006, 53, 297–312. [Google Scholar] [CrossRef]
  53. Bell, L.; Wagstaff, C. Enhancement of glucosinolate and isothiocyanate profiles in Brassicaceae crops: Addressing challenges in breeding for cultivation, storage, and consumer-related traits. J. Agric. Food Chem. 2017, 65, 9379–9403. [Google Scholar] [CrossRef]
  54. Clarke, J.D.; Dashwood, R.H.; Ho, E. Multi-targeted prevention of cancer by sulforaphane. Cancer Lett. 2008, 269, 291–304. [Google Scholar] [CrossRef]
  55. Zhang, Y. Cancer-preventive isothiocyanates: Measurement of human exposure and mechanism of action. Mutat. Res. 2004, 555, 173–190. [Google Scholar] [CrossRef] [PubMed]
  56. Blažević, I.; Montaut, S.; Burčul, F.; Rollin, P. Glucosinolates: Novel sources and biological potential. In Glucosinolates, Reference Series in Phytochemistry; Mérillon, J.M., Ramawat, K.G., Eds.; Springer: Cham, Switzerland, 2017; pp. 3–60. [Google Scholar]
  57. Singh, A. Glucosinolates and plant defense. In Glucosinolates, Reference Series in Phytochemistry; Mérillon, J.M., Ramawat, K.G., Eds.; Springer: Cham, Switzerland, 2017; pp. 237–246. [Google Scholar]
  58. Liu, Z.; Wang, H.; Xie, J.; Lv, J.; Zhang, G.; Hu, L.; Luo, S.; Li, L.; Yu, J. The roles of Cruciferae glucosinolates in disease and pest resistance. Plants 2021, 10, 1097. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Representative photos of Brassica juncea plants (1)–(6) showing the whole plant (a); the abaxial (b–d) and adaxial (e–g) leaf surfaces; and the outer (b,e), middle (c,f), and inner (d,g) leaves.
Figure 1. Representative photos of Brassica juncea plants (1)–(6) showing the whole plant (a); the abaxial (b–d) and adaxial (e–g) leaf surfaces; and the outer (b,e), middle (c,f), and inner (d,g) leaves.
Foods 12 04374 g001
Figure 2. Structure of glucosinolates (top) and representative MRM chromatograms of glucosinolates (bottom) identified and quantified in Brassica juncea accessions.
Figure 2. Structure of glucosinolates (top) and representative MRM chromatograms of glucosinolates (bottom) identified and quantified in Brassica juncea accessions.
Foods 12 04374 g002aFoods 12 04374 g002b
Figure 3. Frequency distribution of individual and total GSLs levels in 355 accessions.
Figure 3. Frequency distribution of individual and total GSLs levels in 355 accessions.
Foods 12 04374 g003
Figure 4. Loading (a) and score plots (be) of B. juncea samples based on quantitative agro-morphological characters. The score plots are grouped using leaf type, plant posture (b), leaf type (c), leaf shape (d), and head formation (e). PC1 and PC2 explained 47.2% and 17.8% of the variations, respectively.
Figure 4. Loading (a) and score plots (be) of B. juncea samples based on quantitative agro-morphological characters. The score plots are grouped using leaf type, plant posture (b), leaf type (c), leaf shape (d), and head formation (e). PC1 and PC2 explained 47.2% and 17.8% of the variations, respectively.
Foods 12 04374 g004aFoods 12 04374 g004bFoods 12 04374 g004c
Figure 5. Inter-leaf (a) and intra-leaf (b) distribution of GSLs in 15 accessions of B. juncea. The letters in the graph indicate the differences between the different results.
Figure 5. Inter-leaf (a) and intra-leaf (b) distribution of GSLs in 15 accessions of B. juncea. The letters in the graph indicate the differences between the different results.
Foods 12 04374 g005
Table 1. Retention time (tR), MRM transitions, collision-induced dissociation values, and calibration equations of the individual GSLs identified in Brassica juncea leaves.
Table 1. Retention time (tR), MRM transitions, collision-induced dissociation values, and calibration equations of the individual GSLs identified in Brassica juncea leaves.
GlucosinolatestRMRM TransitionCIDCalibration Equation
Sinigrin1.08357.88 > 9620Y = 12.8988X + 16.8399, r2 = 0.999
Gluconapin2.58371.89 > 9620Y = 12.0411X − 13.1698, r2 = 0.998
Glucobrassicanapin2.90385.95 > 9625Y = 44.9459X − 37.6357, r2 = 0.997
Glucoiberin1.09421.91 > 9620Y = 8.48762X − 24.3857, r2 = 0.999
Gluconasturtiin3.23421.97 > 9625Y = 39.0356X − 44.9008, r2 = 0.998
Glucobarbarin2.91437.93 > 9620Y = 15.443X − 21.6364, r2 = 0.996
Glucobrassicin3.01446.95 > 9620Y = 28.0663X − 43.3473, r2 = 0.999
tR: Retention time; MRM: Multiple reaction monitoring; CID: Collision-induced dissociation.
Table 2. Modified International Union for the Protection of New Varieties of Plants (UPOV) descriptors used for the qualitative morphological characterization of Brassica juncea.
Table 2. Modified International Union for the Protection of New Varieties of Plants (UPOV) descriptors used for the qualitative morphological characterization of Brassica juncea.
S/NOTraitDescriptions *
1Plant habit (posture)Erect (73), semi-erect (162), horizontal (114), mixed (6)
2Plant: head formationAbsent (5), partially present (297), fully present (52), mixed (1)
3Leaf attitude at the apical partFlat (119), bend downward (215), bend upward (15), mixed (5)
4Leaf typeLyrate (257), division (17), entire and narrow midrib (46), entire and broad midrib (24), mixed (11)
5Leaf shapeLanceolate (17), oblanceolate (16), spatulate (186), ovate (67), obovate (35), elliptic (9), broad elliptic (11), compressed circular (3), mixed (11)
6Leaf waxinessWeak (52), medium (142), strong (148), mixed (13)
7Lobation of leaf marginsAbsent (10), lobed (103), cleft (181), parted (21), dissected (22), mixed (18)
7Leaf blade: density of incisions of marginAbsent (8), sparse (113), medium (164), dense (29), very dense (23), mixed (18)
8Leaf blade: size of the terminal lobeAbsent (16), small (196), medium (106), large (20), mixed (17)
9Leaf blade: number of serratesAbsent (17), few (152), medium (127), many (47), mixed (12)
11Leaf blade: anthocyanin colorationAbsent (27), very weak (121), weak (52), medium (65), strong (42), very strong (24), mixed (24)
12Leaf blade: blisteringAbsent/weak (80), medium (118), strong (143), mixed (14)
13Leaf blade: pubescence on the upper sideAbsent (81), present (263), mixed (11)
14Leaf blade: pubescence on the lower sideAbsent (54), few (101), medium (154), many (35), mixed (11)
15Stem hairAbsent (59), present (291), mixed (5)
16Midrib transectionHorizontal (41), intermediate (122), semi-circled (191)
17Petiole anthocyanin colorationAbsent (77), present (257), mixed (18)
* Values in parentheses indicate the number of accessions in each category.
Table 3. Variation of quantitative morphological characters and GSLs in Brassica juncea leaf samples from 355 accessions.
Table 3. Variation of quantitative morphological characters and GSLs in Brassica juncea leaf samples from 355 accessions.
GSLs and Agro-Morphological CharactersRangeAverage ± SDCoefficient of Variation (%)
SIN13.0–4184.6635.9 ± 586.492.2
GNA44.1–7708.7496.4 ± 628.5126.6
GBN0.1–180.212.8 ± 20.1156.9
GIB0.0–93.62.1 ± 6.3299.5
GBB0.7–33.94.7 ± 3.368.9
Total aliphatic GSLs133.7–7738.21151.8 ± 982.685.3
GNS27.2–3393.5401.9 ± 342.985.3
GBS17.3–351.957.9 ± 44.276.3
Total GSLs320.4–8055.31610.8 ± 1064.866.1
Leaf length (cm)20.8–57.036.3 ± 5.816.0
Leaf width (cm)8.6–31.818.1 ± 4.324.0
Midrib width (cm)0.6–3.81.8 ± 0.633.7
Midrib thickness (cm)0.3–1.10.6 ± 0.118.9
Petiole length (cm)0.9–12.53.1 ± 1.547.7
Petiole width (cm)0.4–3.81.3 ± 0.648.7
Number of leaves (ea)6.0–28.016.3 ± 636.8
Plant weight (g, FW)91.7–1854.3338.9 ± 152.745.1
The contents of glucosinolates are expressed in μmole/kg DW. Glucosinolate abbreviations: GSLs: Glucosinolates; SIN: Sinigrin; GNA: Gluconapin; GBN: Glucobrassicanapin; GIB: Glucoiberin; GBB: Glucobarbarin; GNS: Gluconasturtiin; GBS: Glucobrassicin; FW: Fresh weight; SD: Standard deviation.
Table 4. Pearson correlations among the glucosinolates and quantitative agro-morphological characters of 355 B. juncea accessions.
Table 4. Pearson correlations among the glucosinolates and quantitative agro-morphological characters of 355 B. juncea accessions.
Leaf LengthLeaf WidthNo of LeavesMidrib WidthMidrib ThicknessPetiole LengthPetiole WidthPlant WeightSINGNAGBNGIBGNSGBBGBS
Leaf width0.622 **
Number of leaves0.303 **0.167 **
Midrib width0.444 **0.854 **0.256 **
Midrib thickness0.475 **0.581 **0.161 **0.455 **
Petiole length0.210 **−0.128 *−0.004−0.229 **0.045
Petiole width0.320 **0.829 **0.0040.895 **0.429 **−0.268 **
Plant weight0.532 **0.451 **0.0850.348 **0.352 **0.122 *0.368 **
SIN−0.116 *0.008−0.0090.051−0.082−0.185 **0.055−0.138 **
GNA−0.0610.0810.0120.143 **0.119 *−0.0070.119 *−0.0980.272 **
GBN−0.0490.0650.0630.155 **0.041−0.166 **0.121 *−0.106 *0.346 **0.551 **
GIB−0.119 *−0.0050.0380.04−0.061−0.144 **0.061−0.0440.256 **0.159 **0.139 *
GNS−0.016−0.1020.092−0.053−0.087−0.000−0.101−0.0420.082−0.0850.097−0.058
GBB−0.005−0.0090.0120.048−0.01−0.0360.016−0.0240.139 **−0.0540.0980.0040.565 **
GBS−0.1010.065−0.0150.112 *0.006−0.124 *0.127 *−0.105 *0.463 **0.380 **0.173 **0.0640.110 *0.083
Total GSLs−0.111 *0.0240.0320.103−0.001−0.115*0.076−0.154 **0.766 **0.740 **0.575 **0.226 **0.325 **0.235 **0.560 **
** Correlations are significant at the 0.01 level; * Correlations are significant at the 0.05 level. Glucosinolate abbreviations: GSLs: Glucosinolates; SIN: Sinigrin; GNA: Gluconapin; GBN: Glucobrassicanapin; GIB: Glucoiberin; GBB: Glucobarbarin; GNS: Gluconasturtiin; GBS: Glucobrassicin.
Table 5. The correlations of agro-morphological characters with the levels of GSLs.
Table 5. The correlations of agro-morphological characters with the levels of GSLs.
Morphological CharactersDescriptions (N)SINGNAGBNGIBGNSGBBGBSTotal GSLs
Plant habit (posture)Erect (73)659.2 ab614.9 a12.9 a2.5 a370.4 a5.1 a71.1 b1735.9 ab
Semi erect (162)541.2 a412.3 a11.0 a1.5 a399.5 a4.5 a48.3 a1418.2 a
Horizontal (141)749.5 b537.7 a15.3 a2.7 a424.3 a4.8 a62.8 b1794.6 b
Leaf: attitude at the apical partFlat (119)680.1 a535.6 a12.2 a2.6 a412.7 a4.5 a58.8 a1703.9 a
Bend downward (215)607 a475.1 a12.9 a1.8 a397.8 a4.9 a58.5 a1557.9 a
Bend upward (16)674.9 a486.5 a16.8 a2.9 a406.6 a4.4 a43.9 a1635.3 a
Leaf typeLyrate (257)635.8 ab463.2 a12.3 a1.9 a427.2 a4.8 a55.7 ab1599.7 ab
Division (17)305.1 a431.9 a7.2 a0.3 a355.6 a4.6 a33.1 a1137.7 a
Entire, narrow midrib (46)678.5 ab631.0 a19.0 a3.2 a321.6 a4.2 a65.8 bc1723.1 ab
Entire, broad midrib (24)762.3 b668.8 a13.0 a3.9 a341.7 a5.0 a86.2 c1880.9 b
Leaf shapeLanceolate (17)788.8 a400.8 a8.6 a1.7 a327.7 a4.9 a63.1 ab1595.5 a
Oblanceolate (16)524 a343.1 a7.7 a1.3 a398.7 a4.0 a55.4 ab1334.1 a
Spatulate (186)633.6 a467.9 a11.5 a2.1 a441.0 a4.9 a56.3 ab1615.6 a
Ovate (67)665.5 a688.7 a20.5 a2.1 a378.8 a4.6 a57.4 ab1817.3 a
Obovate (35)510.0 a476.4 a13.9 a1.4 a286.0 a3.6 a64.5 ab1355.7 a
Elliptic (9)590.9 a383.8 a4.2 a2.0 a409.5 a5.0 a43.9 a1439.4 a
Broad elliptic (11)714.1 a261.0 a4.0 a5.3 a374.2 a5.7 a65.3 ab1429.4 a
Compressed circular (3)694.8 a676.6 a13.0 a1.5 a395.5 a7.1 a109 b1897.5 a
Lobation of leaf marginsAbsent (10)600.1 a571.5 a21.1 b2.1 a298.0 a4.0 a60.1 ab1556.6 ab
Lobed (103)718.9 a565.0 a16.6 ab2.4 a397.0 a4.9 a64.2 ab1768.9 ab
Cleft (181)592.5 a474.7 a11.0 ab2.0 a418.2 a4.8 a56.3 ab1557.6 ab
Parted (21)840.5 a668.7 a16.2 ab3.1 a387.5 a4.1 a72.8 b1993.0 b
Dissected (22)416.4 a199.5 a3.6 a0.6 a351.9 a3.6 a37.4 a1012.8 a
Leaf blade: density of incisions of marginsAbsent (8)614.3 a600.5 a24.4 b2.2 a335.1 a4.1 a67.2 a1647.8 a
Sparse (113)726.9 a543.3 a16.5 ab2.3 a403.4 a4.8 a65.1 a1762.1 a
Medium (164)591.7 a474.4 a10.9 ab2.2 a415.0 a4.9 a55.8 a1552.8 a
Dense (29)680.1 a557.5 a12.9 ab2.1 a377.0 a4.2 a64.9 a1698.7 a
Very dense (23)447.3 a351.5 a4.3 a0.8 a355.8 a3.7 a36.2 a1199.6 a
Leaf blade: size of terminal lobeAbsent (16)755.3 a560.7 a17.9 b1.9 a296.7 a4.4 a70.5 a1707.3 a
Small (196)623.0 a529.6 a14.2 ab2.3 a408.9 a4.6 a60.5 a1642.5 a
Medium (106)660.9 a477.4 a11.0 ab2.1 a385.7 a4.9 a52.9 a1592.7 a
Large (20)510.2 a253.0 a4.7 a0.7 a497.0 a5.0 a59.7 a1330.2 a
Leaf blade: number of serratesAbsent (17)906.1 a514.2 a17.4 a4.8 a373.4 a5.5 a77.0 b1898.3 a
Few (152)596.2 a490.8 a13.8 a1.6 a405.2 a4.7 a55.6 ab1566.2 a
Medium (127)670.2 a556.9 a12.6 a1.9 a394.3 a4.7 a62.4 ab1702.9 a
Many (47)595.6 a383.8 a9.8 a3.4 a423.6 a4.4 a49.4 a1468.3 a
Leaf blade: anthocyanin colorationAbsent (27)453.7 a296.3 a5.6 a1.5 a320.4 a3.9 a38.5 a1119.7 a
Very weak (121)649.0 a645.4 a16.7 a2.7 a362.0 a4.7 a64.9 a1744.8 a
Weak (52)704.1 a410.5 a13.2 a1.6 a513.3 a5.3 a63.9 a1707.3 a
Medium (65)742.5 a568.3 a15.4 a2.8 a367.8 a4.4 a59.3 a1760.3 a
Strong (42)510.0 a342.0 a6.3 a1.3 a420.5 a4.5 a50.9 a1335.2 a
Very strong (24)497.5 a342.2 a8.0 a0.5 a515.3 a5.8 a49.8 a1418.9 a
Leaf waxinessWeak (52)651.4 a367.6 a11.3 ab1.1 a449.4 a4.4 a63.6 b1548.6 ab
Medium (142)528.9 a417.0 ab9.6 a2.0 a377.2 a4.6 a45.5 a1382.4 a
Strong (148)721.9 a619.7 b16.6 b2.5 a410.1 a5.0 a68.9 b1844.6 b
Leaf blade: blisteringAbsent/weak (80)663.2 a465.3 a13.4 a1.9 a415.5 a4.6 a61.9 a1625.8 a
Medium (118)590.0 a518.3 a13.1 a1.7 a368.9 a4.9 a57.5 a1552.3 a
Strong (143)647.6 a496.0 a12.2 a2.5 a419.8 a4.6 a57.1 a1639.1 a
Leaf blade: pubescence on the upper sideAbsent (81)668.8653.720.12.8319.94.968.81738.9
Present (263)633.5449.310.81.9427.24.755.31581.5
Leaf blade: pubescence on the lower sideAbsent (54)677.9 b661.5 b20.5 b3.2 a310.7 a5.0 a74.1 b1752.8 b
Few (101)687.0 b483.3 ab14.2 ab1.7 a385.5 a4.9 a61.2 b1637.7 ab
Medium (154)659.5 b478.9 ab11.0 a2.4 a449.0 a4.8 a55.9 ab1659.3 ab
Many (35)378.1 a366.5 a7.2 a0.6 a383.0 a3.7 a37.9 a1177.0 a
Stem hairAbsent (59)670.9652.420.23.2321.25.073.81746.5
Present (291)630.8461.011.41.9417.94.755.01581.6
Midrib
transection
Horizontal (41)837.8 b684.4 b22.8 b2.4 a434.0 a5.2 a73.1 b2059.7 b
Intermediate (122)589.1 a541.8 ab12.8 a2.0 a420.8 a4.8 a55.3 a1626.6 a
Semi-circled (191)604.4 a422.8 a10.6 a2.0 a383.3 a4.6 a56.2 a1482.1 a
Petiole
anthocyanin coloration
Absent (77)468.6418.88.71.0385.14.343.21329.7
Present (259)683.5531.914.42.5411.94.862.71710.5
Plant
head formation
Absent (5)263.5 a216.4 a5.5 a0.2 a326.0 a3.5 a30.9 a846.0 a
Partially present (297)635.1 a483.5 a12.6 a2.0 a421.4 a4.8 a55.7 ab1613.9 a
Fully present (52)697.7 a619.4 a15.6 a3.1 a301.8 a4.5 a75.0 b1716.9 a
N = number of accessions. Glucosinolate abbreviations: GSLs: Glucosinolates; SIN: Sinigrin; GNA: Gluconapin; GBN: Glucobrassicanapin; GIB: Glucoiberin; GBB: Glucobarbarin; GNS: Gluconasturtiin; GBS: Glucobrassicin. Different letters within the same column and the same morphological characters denote significant differences (p < 0.05). Values in parentheses indicate the number of accessions in each category.
Table 6. Eigenvalues, proportions of variability, and agro-morphological traits contributed to the first five principal components of Brassica juncea genetic resources.
Table 6. Eigenvalues, proportions of variability, and agro-morphological traits contributed to the first five principal components of Brassica juncea genetic resources.
TraitsPC1PC2PC3PC4PC5
Leaf length0.360660.425130.050563−0.12065−0.014741
Leaf width0.48364−0.08839−0.0536830.0992520.13102
No of leaves0.140670.231870.90222−0.067980.070513
Midrib width0.45814−0.242690.0900630.0722380.3238
Midrib thickness0.350290.16403−0.0553670.54633−0.70066
Petiole length−0.0555340.67376−0.256850.372610.51955
Petiole width0.43236−0.34385−0.163080.0534910.27198
Plant weight0.305980.31576−0.27717−0.72497−0.19486
Eigenvalue3.780611.434250.9974870.5725430.398334
% Variance47.25817.92812.4697.15684.9792
Table 7. Top ten best accessions identified for different agro-morphological and biochemical traits.
Table 7. Top ten best accessions identified for different agro-morphological and biochemical traits.
S/NoTraitsTop Ten Accessions *
1Leaf length (≥50 cm)IT120115, IT118972, IT118979, IT203623, IT228987, IT102888, Yesou86, IT215792, and IT120114
2Leaf width (≥28 cm)Yeosu84, Yeosu 42, Yeosu61, IT228223, IT204152, Yeosu27, Yeosu 57, IT203623, Yeosu26, and Yeosu20
3No of leaves per plant (≥26 leaves)IT228222, IT204153, IT215792, IT218388, IT235416, IT216873, IT118948, IT235418, IT180994, IT235344, IT216872, IT218451, IT218458, IT118955, and IT218449
4Midrib width (≥3.15 cm)Yeosu80, IT228223, IT203623, Yeosu57, Yeosu15, IT204152, Yeosu84, IT228990, IT228988, and Yeosu61
5Midrib thickness (≥0.75 cm)IT218357, Yeosu42, Yeosu27, IT141424, Yeosu41, IT102888, IT208801, Yeosu6, Yeosu15, IT118974, Yeosu11, and IT259499
6Petiole length (≥6.3 cm)IT102941, IT102942, IT259517, IT100949, IT141423, IT228987, IT259499, IT100945, IT218383, and IT141422
7Petiole width (≥2.9 cm)Yeosu57, Yeosu80, IT228223, Yeosu61, IT191096, IT228988, IT228988, IT204152, Yeosu21, and Yeosu24
8Plant weight (≥650 g)IT228984, Yeosu45, IT118979, IT215792, IT228223, IT100949, IT248039, Yeosu42, Yeosu64, IT228987
9Sinigrin (≥2295)IT237840, IT259503, IT248036, IT259487, IT236762, IT248037, IT250121, IT247853, IT235418, and Yeosu17
10Gluconapin (≥2000)Yeosu66, Yeosu31, IT102941, Yeosu17, IT218460, IT191097, Yeosu41, Yeosu32, Yeosu2, IT250121, and IT259487
11Glucobrassicanapin (≥70)Yeosu31, Yeosu37, Yeosu53, Yeosu41, Yeosu23, Yeosu17, Yeosu32, Yeosu44, Yeosu52, and Yeosu75
12Glucoiberin (≥16)IT102894, Yeosu74, Yeosu2, IT259503, IT250121, Yeosu63, Yeosu52, IT259514, IT259492, and Yeosu38
13Gluconasturtiin (≥1340)IT236761, IT248035, IT259491, IT248614, IT237839, IT235624, IT235803, Yeosu16, IT259513, and Yeosu23
14Glucobarbarin (≥13)IT102937, IT102950, IT236761, IT109154, IT218453, IT248035, Yeosu43, IT235415, Yeosu80, and IT237839
15Glucobrassicin (≥180)IT248036, IT248037, Yeosu66, Yeosu20, Yeosu19, Yeosu80, Yeosu21, IT259487, Yeosu18, and IT259509
16Total Glucosinolates (≥4600)Yeosu66, IT259487, IT259503, IT237840, Yeosu17, IT248036, IT248035, IT250121, IT248037, and Yeosu31
* Samples are arranged in decreasing order of their values; the GSLs levels are reported in μmole/kg DW.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Assefa, A.D.; Kim, S.-H.; Ko, H.C.; Ro, N.; Subramanian, P.; Chung, Y.-J.; Lee, Y.-H.; Hahn, B.-S.; Rhee, J.-H. Leaf Mustard (Brassica juncea) Germplasm Resources Showed Diverse Characteristics in Agro-Morphological Traits and Glucosinolate Levels. Foods 2023, 12, 4374. https://doi.org/10.3390/foods12234374

AMA Style

Assefa AD, Kim S-H, Ko HC, Ro N, Subramanian P, Chung Y-J, Lee Y-H, Hahn B-S, Rhee J-H. Leaf Mustard (Brassica juncea) Germplasm Resources Showed Diverse Characteristics in Agro-Morphological Traits and Glucosinolate Levels. Foods. 2023; 12(23):4374. https://doi.org/10.3390/foods12234374

Chicago/Turabian Style

Assefa, Awraris Derbie, Seong-Hoon Kim, Ho Chul Ko, Nayoung Ro, Parthiban Subramanian, Yun-Jo Chung, Yong-Hyuk Lee, Bum-Soo Hahn, and Ju-Hee Rhee. 2023. "Leaf Mustard (Brassica juncea) Germplasm Resources Showed Diverse Characteristics in Agro-Morphological Traits and Glucosinolate Levels" Foods 12, no. 23: 4374. https://doi.org/10.3390/foods12234374

APA Style

Assefa, A. D., Kim, S. -H., Ko, H. C., Ro, N., Subramanian, P., Chung, Y. -J., Lee, Y. -H., Hahn, B. -S., & Rhee, J. -H. (2023). Leaf Mustard (Brassica juncea) Germplasm Resources Showed Diverse Characteristics in Agro-Morphological Traits and Glucosinolate Levels. Foods, 12(23), 4374. https://doi.org/10.3390/foods12234374

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop