Multivariate Based Variability within Diverse Indian Mustard ( Brassica juncea L . ) Genotypes

Indian mustard (Brassica juncea L.) germplasm consisting 167 accessions including one check cultivar was evaluated for qualitative and quantitative traits. The present study was conducted to investigate genetic diversity and correlation among studied genotypes of B. juncea L. based on agro-morphological at NARC, Islamabad, Pakistan. To investigate the genetic diversity based on morphological characters, data was recorded on 20 quantitative and 12 qualitative traits. The calculated data was analyzed through two complementary methods, i.e. PCA (Principal Component Analysis) and cluster analysis. Among all the studied cultivars, significant diversity was recorded for different agro-morphological characters. Among all the parameters, maximum variance was recorded for pod shattering (427.2) followed by plant height (345.6), days to 100% flowering (336.2) and main raceme length (210.0). Among all the characters, the greatest and highly significant association (0.99) was found between days to maturity 50% and days to maturity 100% followed by correlation (0.86) among days to flowering 50% and days to flowering 100%, correlation value (0.71) was calculated among leaf length and leaf width. Using cluster analysis all the genotypes were divided into five major groups. It was observed that 7 out of 20 principal components with an Eigen value of ≥1.0 calculated for 73.92% of the total diversity observed between 167 accessions of Indian mustard (B. juncea L.). The contribution of first three PCs in the total PCs was 23.25, 12.87 and 11.24, respectively. Among all the investigated accessions two genotypes 26,813 and 26,817 showed great How to cite this paper: Saleem, N., Jan, S.A., Atif, M.J., Khurshid, H., Khan, S.A., Abdullah, M., Jahanzaib, M., Ahmed, H., Ullah, S.F., Iqbal, A., Naqi, S., Ilyas, M., Ali1, N. and Rabbani, M.A. (2017) Multivariate Based Variability within Diverse Indian Mustard (Brassica juncea L.) Genotypes. Open Journal of Genetics, 7, 69-83. https://doi.org/10.4236/ojgen.2017.72007 Received: February 11, 2017 Accepted: May 7, 2017 Published: May 10, 2017 Copyright © 2017 by authors and Scientific Research Publishing Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/ Open Access

ysis.Among all the studied cultivars, significant diversity was recorded for different agro-morphological characters.Among all the parameters, maximum variance was recorded for pod shattering (427.2) followed by plant height (345.6), days to 100% flowering (336.2) and main raceme length (210.0).Among all the characters, the greatest and highly significant association (0.99) was found between days to maturity 50% and days to maturity 100% followed by correlation (0.86) among days to flowering 50% and days to flowering 100%, correlation value (0.71) was calculated among leaf length and leaf width.Using cluster analysis all the genotypes were divided into five major groups.It was observed that 7 out of 20 principal components with an Eigen value of ≥1.0 calculated for 73.92% of the total diversity observed between 167 accessions of Indian mustard (B.juncea L.).The contribution of first three PCs in the total PCs was 23.25, 12.87 and 11.24, respectively.Among all the investigated accessions two genotypes 26,813 and 26,817 showed great

Introduction
Brassica juncea L. commonly known as Indian mustard is globally used as oilseed, vegetable and condiments.The crop species is said to be one of the earliest domesticated species.B. juncea annually grows as wild and cultivated species over a large geographical range, across the Asia, Africa, Australia, America and Europe [1].In the genus Brassicaceae, Indian mustard is one of the highly significant crop species.In Pakistan after cotton and rape seed, Indian mustard is the third highly important source of oil [2] [3].
So for the prosperity of the country, it is essential to increase the production of these crops and for maximum production, exploitation of genetic diversity of these crops is essential.Pakistan being an agricultural country but unfortunately it is deficient in the oilseed production and a large amount ($ 1054.7 million) is spent on the import about 77% of edible oil.After petroleum, Pakistan pays the second largest import bill for edible oil among food items [4].
Diversity can be determined by three different levels, first one is genetic diversity, in which the variation observed at gene and germplasm level, the second one is diversity of species, in which richness of species at a specific location and the third one is the Eco-system diversity, in which crop species interact with their environment.Diversity is potent for different crops to fulfill the gap between demand and production across the world [5].Various morpho-biochemical and molecular methods are used to study genetic variability among local and exotic plant germplasm.The proper evaluation of important crop species helps in the identification and utilization of improved genotypes [6]- [11].Therefore, the present study was performed to estimate and characterize the genetic variation and relationships among Indian mustard genotypes through agro-morphological traits to identify and select promising germplasm for traits of economic significance.Row to row distance of 70 cm was maintained with one accession per row and the length of the row was 2.5 meter for each and every accession and genotype.

Materials and Methods
During seed bed preparation the field was irrigated before sowing to provide conditions of optimum moisture while afterword no irrigations were given and the crop was grown under rain-fed conditions.The experimental germplasm were planted at a depth of 3 -4 cm using hand drill.After germination to maintain proper plant population thinning was carried out by hand.After 30 days of germination weeding was carried out and because of no serious attack of insect pest no insecticides and pesticides were sprayed throughout the plant growth stage.The crop was harvested when almost more than 75% of plants were turned yellowish in color.Data were recorded for twenty quantitative traits such as Days to 50% flowering, Days to 100% flowering, Days to 50% maturity, Days to 100% maturity, Leaf length, Leaf width, Leaf width/leaf length ratio, Leaves per plant, Plant height, Primary branches per plant, Main raceme length, Silique per main raceme, Silique length, Silique width, Silique length and silique width ratio, Stem thickness, Seed per silique, 1000-seed weight, Seed yield per plant, Shattering percentage (%).The description of calculating different parameters is given in Table 1.The recorded morphological traits data were averaged and analyzed for simple statistics i.e. mean, variance, range, frequency distribution, coefficient of variance and standard deviation using computer software (MS-Excel 2007).
Cluster and principal component analysis (PCA) was performed on the recorded data for quantitative traits.Before to cluster and PCA, mean of each parameter was standardized so that to avoiding scaling differences effects.For all the pairs of accessions Euclidean distance co-efficient were calculated.The Euclidean dissimilarity co-efficient matrices were utilized to estimate the association among the B. juncea germplasm with cluster analysis through complete linkage method (Statistica version 7.0 and NTSYS pc v 2.1).

Principal Component Analysis Based on Agronomic and Morphological Traits
Principal component analysis (PCA) was carried out based on twenty quantitative morphological characters.The 7 principal components account for 73.92% of the overall variability among the studied B. juncea L. accessions for the total phenotypic variations (Table 3, Figures 1-6).Among these seven principal components (PCs), the PCA-I was found to have 23.35% out of the total varia-
Our results were strengthened by many earlier plant breeders and scientists.
The highly significant and positive correlation was observed between the parameters of days to flowering 50% and 100% and days to maturity 50% and 100% physiological maturity.These results were supported by the results of Nasim et al. [13], Alemayehu and Becker [16], by reporting highly significant and positive correlation among days to 50% and 100% flowering and days to 50% and 100% maturity.The results obtained by Zare [17] were also strongly supported our finding, who reported positive and highly significant correlation among important traits like days to flowering, days to maturity and plant height of B. napus L. accessions.Similarly Azadgoleh et al. [18] also reported positive and highly significant association among days to physiological maturity and plant height.Similarly the results of Rabbani et al. [19] also strengthen our findings.In the current investigation main raceme length was highly significant positive correlated with the trait of number of silique per main raceme.The results obtained by Rabbani et al. [19] and Nasim et al. [13] had also confirmed our findings.Plant height was also found positive and highly significant correlated to primary branches per plant, main raceme length and silique per main raceme which was also close agreement by previous research of Yasir et al. [20] who studied strong positive and significant association of plant height with primary branches per plant, main raceme length, pod on main raceme, pod length and seeds per pod.
Numerous plant breeders and plant scientist in the earlier era had practiced and got remarkable results of diversity in agro-morphological parameters for various Brassica individuals through two complementary methods i.e. cluster analysis and principal component analysis (PCA).These techniques were successfully observed in Indian mustard (Brassica juncea L.) germplasm by Gupta et al. [21], in Ethiopian mustard by Alemayehu and Becker [16] and in white head cabbage Balkaya et al. [22].All the above mentioned plant breeders and plant scientist's results gives support to our present investigation that these 2 techniques are very supportive in estimated associations among individuals originated from different environments in a more clear approach.The classification of present studied Brassica juncea L. genotypes into various groups was not due to of its origination from diverse ecological zones of the country and world.
Amurrio et al. [23], who investigated qualitative and quantitative traits of Iberian pea genotypes, were agreement with our results.During the present investigation among different accessions of Brassica juncea L. considerable level of genetic differences was observed for different agronomic and morphological characters.These studied parameters can play significant and interesting role in Brassica juncea L. cultivars and varieties development programs in future.In present study the seven groups were recorded and PC1 showed highest level of variability.The same findings were also recorded by Dias et al. [24] in diverse germplasm of Brassica oleracea L. originated from Portugal.

Conclusion
Among all the characters, valuable genetic diversity was observed through different agro-morphological traits i.e. days to maturity 50% and days to maturity 100%, days to 50% flowering and days to flowering 100%, leaf length and leaf width, main raceme length siliquae per main raceme, leaf width, primary branches per plant, silique length.Most of the traits contributed highly significant association for the said trait during the present investigation.These genotypes show that there is great potential in the studied accessions, which can be used in future breeding program as a productive genotypes.
Plant material comprises 167 germplasm including 1 check cultivar of B. juncea L. for estimation of agro-morphological variation.The germplasm were provided by PGRI gene-bank, NARC, Islamabad, Pakistan.The present research work was conducted in the experimental area of NARC (33˚43'N and 73˚06'E)during the year 2012-2013 using augmented design.The average annual rainfall in this area ranges from 500 -900 mm with 70% in summer and 30% in winter.

Figure 1 .
Figure 1.Scatter diagram og PC-I and II for agro-morphological traits.

Figure 2 .
Figure 2. Scatter diagram of PC-I and III for agro-morphological traits.

Figure 3 .
Figure 3. Scatter diagram of PC-II and III for agro-morphological traits.

Figure 4 .
Figure 4. Contribution of morphological traits in the first two principal components.

Figure 5 .
Figure 5. Contribution of morphological traits in the second and third principal components.

Figure 6 .
Figure 6.Contribution of morphological traits in the first and third principal components.

Table 1 .
Detailed information for calculating 20 different quantitative traits of B. juncea L. genotypes.