Open Journal of Genetics
Vol.07 No.02(2017), Article ID:76076,15 pages
10.4236/ojgen.2017.72007

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

Noor Saleem1, Sohail Ahmad Jan2,3*, Muahmmad Jawaad Atif1, Haris Khurshid4, Shahid Ali Khan2, Mohib Abdullah5, Muhammad Jahanzaib4, Humair Ahmed6, Syed Farman Ullah7, Azhar Iqbal7, Sameena Naqi7, Muhammad Ilyas7, Naushad Ali1, Malik Ashiq Rabbani2

1The University of Haripur, KPK, Pakistan

2Plant Genetic Resources Institute (PGRI), National Agricultural Research Centre (NARC), Islamabad, Pakistan

3Department of Biotechnology, Quaid-i-Azam University, Islamabad, Pakistan

4Oilseeds Research Program, NARC, Islamabad, Pakistan

5School of Crop Production Technology, Suranaree University of Technology, Nakhon Ratchasima, Thailand

6Land Resources Research Institute, NARC, Islamabad, Pakistan

7Department of Plant Breeding and Genetics, PMAS-Arid Agriculture University, Rawalpindi, Pakistan

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/

Received: February 11, 2017; Accepted: May 7, 2017; Published: May 10, 2017

ABSTRACT

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-morpho- logical 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 potential for seeds/silique, 1000-seed weight and seed yield/plant, respectively, so these genotypes are recommended for future breeding programs for achieving promising results.

Keywords:

Genetic Diversity, Agro-Morphology, Brassica juncea L., PCA Analysis

1. 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-biochem- ical 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-mor- phological traits to identify and select promising germplasm for traits of economic significance.

2. Materials and Methods

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. 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. 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).

3. Results

3.1. Correlation Study Based on Agronomic and Morphological Traits

During the present research work days to 50% flowering exhibited highly significant positive correlation with the parameter i.e. days to 100% flowering (0.86**), leaf length (0.28**) and leaves per plant (0.21**), similarly days to 50% flowering, had significant positive correlation with days to maturity 100% and leaf length/leaf width ratio (0.18*), plant height (0.17*), days to maturity 50%, leaf width and primary branches per plant (0.16*). While days to flowering 50% had significant negative correlation was observed with silique width (−0.18*). Days to 100% flowering had highly positive significant correlation with leaf length (0.37**), days to maturity 100% (0.29**), days to maturity 50%, leaves per plant and stem thickness (0.27**) and primary branches per plant (0.24**). While the trait like plant height (0.19*) and leaf width (0.18*) had found significant positive correlation with days to 100% flowering (Table 2).

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

Table 2. Correlation coefficient among agro-morphological characters in B. juncea L. accessions.

Note: *significant at 0.05 (0.15 to 0.19 value); **highly significant at 0.20 & above probability level (more than 0.00 level); The green, yellow, white, red colors shows strong positive, weak positive, moderate positive, and strong negative correlations among different traits respectively.

3.2. 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 variability. Silique length (0.073), silique length/width ratio (0.224) and shattering percentage (0.506) contributed positively to first principal component. In contrast, days to flowering 50% and 100% (−0.385 and −0.494), days to maturity 50% and 100% (−0.680 and −0.679), leaf length (−0.711), leaf width (−0.569), leaf length/leaf width ratio (−0.112), leaves per plant (−0.616), stem thickness (−0.606), plant height (−0.513), primary branches per plant (−0.645), main raceme length (−0.391), siliqua per main raceme (−0.653), silique width (−0.193), seeds per silique (−0.230), 1000-seed weight (−0.203) and seed yield per plant (−0.251) contributed negatively.

The contribution of PC-II in total differences was 12.87% and was positively associated with days to maturity 50 and 100% (0.520 and 0.540), leaf length/ width ratio (0.540), silique on main raceme (0.034), silique length (0.234), silique width (0.435), seed per silique (0.530), 1000-seed weight (0.270) and seed yield per plant (0.394), whereas days to flowering 50 and 100% (−0.176 and −0.133), leaf length (−0.149), leaf width (−0.506), leaves per plant (−0.187), stem thickness (−0.153), plant height (−0.580), primary branches per plant (−0.280), main raceme length (−0.243), silique length/width ratio (−0.245) and shattering percentage (−0.215) were negative associated with PC-II.

PC-III showed 11.24% of the total agro-morphological variation and positively associated with leaf width (0.161), leaves per plant (0.033), plant height (0.144), primary branches per plant (−0.015), main raceme length (0.365), siliqua per main raceme (0.266), seeds per silique (0.205), 1000-seed weight (0.232) and silique width (0.578), while negatively associated with days to flowering 50%

Table 3. Principal components (PCs) of agronomic traits among B. juncea L. genotypes.

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

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

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

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

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

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

and 100% (−0.610 and −0.587), days to maturity 50% and 100% (−0.183 and −0.207), leaf length (−0.066), leaf length/leaf width ratio (−0.353), stem thickness (−0.029), silique length (−0.290), silique length/width ratio (−0.696), seed yield per plant (−0.248) and shattering percentage (−0.119).

The PC-IV contributed 8.07% in the total variability and were positively contributed to days to maturity 50% and 100% (0.273 and 0.252), leaf width (0.093), stem thickness (0.193), plant height (0.133), primary branches per plant (0.243), main raceme length (0.095), silique per main raceme (0.040), silique length (0.592), silique length/width ratio (0.490) and seed yield per plant (0.236). Similarly the traits i.e. days to flowering 50% and 100% (−0.463 and −0.443), leaf length (−0.186), leaf length/width ratio (−0.423), leaves per plant (−0.180), silique width (−0.108), seed per silique (−0.103), 1000-seed weight (−0.116) and shattering percentage (−0.069) contributed negatively.

In the total variability the contribution of PC-V was observed 6.89% and the parameter i.e. leaf length (0.329), leaf width (0.103), leaf length/width ratio (0.274), main raceme length (0.246), siliqua per main raceme (0.319), silique length (0.228), silique length/width ratio (0.349), seeds per silique (0.332) were positively contributed. Whereas the characters i.e. days to flowering 50 and 100% (−0.054 and −0.135), days to maturity 50% and 100% (−0.151 and −0.147), leaves per plant (−0.026), stem thickness (−0.492), plant height (−0.118), primary branches per plant (−0.033), silique width (−0.253), 1000-seed weight (−0.007), seed yield per plant and shattering percentage (−0.388 and −0.445) contributed negatively.

The PC-VI depicted proportion of variability as 6.35% and were positively associated with days to flowering 50% and 100% (0.138 and 0.095), days to maturity 50 and 100% (0.063 and 0.037), plant height (0.155), main raceme length (0.563), siliqua per main raceme (0.394), silique length/width ratio (0.090), 1000-seed weight (0.378), and seed yield per plant (0.125). Similarly the traits i.e. leaf length (−0.396), leaf width (−0.417), leaf length/width ratio (−0.021), leaves per plant (−0.196), stem thickness (−0.025), primary branches per plant (−0.160), silique length (−0.224), silique width (−0.259), seeds per silique (−0.238) and shattering percentage (−0.069) were negatively associated with PC-VI.

PC-VII accounted 5.15% of the total variation among agro-morphological traits. Days to maturity 50% and 100% (0.022 and 0.016), leaf length (0.004), leaf length/width ratio (0.192), leaves per plant (0.228), stem thickness (0.165), plant height (0.106), silique length/width ratio (0.022), seeds per silique (0.059) and seed yield per plant (0.248) were contributed positively to PC-VII. Whereas days to flowering 50% and 100% (−0.309 and −0.251), leaf width (−0.107), primary branches per plant (−0.116), main raceme length (−0.143), siliqua per main raceme length (−0.030), silique length (−0.509), silique width (−0.414), 1000-seed weight (−0.338) and shattering percentage (−0.295) depicted negative contribution.

4. Discussion

Maximum variation was calculated for shattering (427.2) followed by plant height (345.6), days to flowering 100% (336.2), main raceme length (210.0), days to flowering 50% (185.1), days to maturity 100% (122.9), days to maturity 50% (83.8). The results obtained from the present investigation were in closed agreement with that of [12] [13] [14] [15] .

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.

5. 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.

Acknowledgements

We acknowledge cooperation of Gene bank of PGRI, NARC, Islamabad, Pakistan for providing B. juncea germplasm.

Conflict of Interest

The author’s declares that there is no conflict of interests regarding the publication of this article.

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

References

  1. 1. Hopkins, C.J., Cogan, N.O.I., Hand, M., Jewell, E., Kaur, J., Li, X., Lim, G.A.C., Ling, A.E., Love, C., Mountford, H., Todorovic, M., Vardy, M., Spangenberg, G.C., Edwards, D. and Batley, J. (2007) Sixteen New Simple Sequence Repeat Markers from Brassica juncea Expressed in Sequences and Their Cross-Species Amplification. Molecular Ecology Notes, 7, 697-700.
    https://doi.org/10.1111/j.1471-8286.2007.01681.x

  2. 2. Rabbani, M.A., Iwabuchi, A., Murakami, Y., Suzuki, T. and Takayanagi, K. (1998) Phenotypic Variation and the Relationships among Mustard (Brassica juncea) Germplasm from Pakistan. Euphytica, 101, 357-366.
    https://doi.org/10.1023/A:1018305201279

  3. 3. Abbas, S.J., Farhatullah, M.K.B., Khan, I.A. and Munir, I. (2009) Molecular Analysis of Genetic Diversity in Brassica Species. Pakistan Journal of Botany, 41, 167-176.

  4. 4. Ahmad, M., Naeem, M. and Khan, I.A. (2013) Relative Abundance of Aphids Population on Different Brassica Genotypes. Sarhad Journal of Agriculture, 29, 133-138.

  5. 5. Shiva, V. (1994) Agriculture and Food Production. UNESCO/Environment Education Dossiers No. 9 (May), 2-3.

  6. 6. Jan, S.A., Shinwari, Z.K. and Rabbani, M.A. (2016) Morpho-Biochemical Evaluation of Brassica rapa Sub-Species for Salt Tolerance. Genetika, 8, 323-338.
    https://doi.org/10.2298/GENSR1601323J

  7. 7. Jan, S.A., Shinwari, Z.K. and Rabbani, M.A. (2016) Determining Genetic Divergence among Brassica rapa Ecotypes through Electrophoretic Mobility of Total Seed Proteins. Journal of Animal and Plant Sciences, 26, 1758-1764.

  8. 8. Khan, Q., Khan, A.S., Khurshid, H., Jan, S.A., Ahmad, N., Khan, S.A., Shah, S.H., Saleem, N., Ibrahim M.I., Ilyas, M. and Arif, M. (2016) Exploring Durable Genetic Resistance against Leaf Rust through Phenotypic Characterization and LR34 Linked STS Marker in Wheat Germplasm. Bioscience Journal, 32, 986-998.
    https://doi.org/10.14393/BJ-v32n4a2016-34046

  9. 9. Jan, S.A., Shinwari, Z.K., Rabbani, M.A., Shah, S.H., Ibrahim, M.I. and Ilyas, M. (2016) Optimization of an Efficient SDS-PAGE Protocol for Rapid Protein Analysis of Brassica rap. Journal of Biodiversity and Environmental Sciences, 9, 17-24.

  10. 10. Jan, S.A., Shinwari, Z.K. and Rabbani, M.A. (2016) Agro-Morphological and Physiological Responses of Brassica rapa Ecotypes to Salt Stress. Pakistan Journal of Botany, 48, 1379-1384.

  11. 11. González-Rubio, C., Rodríguez-Estrella, R., Lozano-Garza, O.A. and García-De-León, F.J. (2016) Genetic Diversity of the Endemic Xantus’ Hummingbird Using 16 Novel Polymorphic Microsatellite Loci, and Their Cross Amplification between Six Related Species. Open Journal of Genetics, 6, 19-27.
    https://doi.org/10.4236/ojgen.2016.61003

  12. 12. Khan, S.A., Iqbal, J., Khurshid, H., Saleem, N., Rabbani, M.A., Zia, M. and Shinwari, Z.K. (2014) The Extent of Intra-Specific Genetic Divergence in Brassica napus L. Population Estimated through Various Agro-Morphological Traits. European Academic Research, 2, 2254-2275.

  13. 13. Nasim, A., Ferhatullah, I.S., Shah, S. and Azam, S.M. (2013) Genetic Variability and Correlation Studies for Morpho-Physiological Traits in Brassica napus L. Pakistan Journal of Botany, 45, 1229-1234.

  14. 14. Vinu, V., Singh, N., Vasudev, S., Kumar, D., Kumar, Y.S., Naresh, S., Bhat, S.R. and Prabhu, K.V. (2013) Assessment of Genetic Diversity in Brassica juncea (Brassicaceae) Genotypes Using Phenotypic Differences and SSR Markers. Revista de Biología Tropical, 61, 1919-1934.

  15. 15. Abideen, S.N., Nadeem, F. and Abideen, S.A. (2013) Genetic Variability and Correlation Studies in Brassica napus L. Genotypes. International Journal of Innovation and Applied Studies, 2, 574-581.

  16. 16. Alemayehu, N. and Becker, H. (2002) Genotypic Diversity and Patterns of Variation in a Germplasm Material of Ethopian mustard (Brassica carinata A. Braun). Genetic Resources and Crop Evolution, 49, 573-582.
    https://doi.org/10.1023/A:1021204412404

  17. 17. Zare, M. (2011) Interrelationship between Grain Yield and Related Triats in Rapeseed (Brassica napus L.). African Journal of Agriculture Research, 6, 6684-6689.

  18. 18. Azadgoleh, E.M.A., Zamani, M. and Esmaeil, Y. (2009) Agronomical Important Traits Correlation in Rapeseed (Brassica napus L.) Genotypes. Research Journal of Agriculture & Biological Sciences, 5, 798-802.

  19. 19. Rabbani, M.A., Iwabuchi, A., Murakami, Y., Suzuki, T. and Takayanagi, K. (1999) Collection, Evaluation and Utilization of Oilseed Mustard (Brassica juncea L.) in Pakistan. Pakistan. Journal of Biological Sciences, 2, 88-94.
    https://doi.org/10.3923/pjbs.1999.88.94

  20. 20. Yasir, A., Farhatullah, R.H., Nasim, A., Azam, S.M. and Khan, A. (2013) Heritability and Correlation Analysis for Morphological and Biochemical Traits in Brassica carinata. Sarhad Journal o Agriculture, 29, 359-370.

  21. 21. Gupta, V.P., Sekhon, M.S. and Satija, D.R. (1991) Studies on Genetic Diversity, Heterosis and Combining Ability in Indian Mustard (Brassica juncea (L.) Czern. & Coss). Indian Journal of Genetics, 51, 448-453.

  22. 22. Balkaya, A., Yanmaz. R. and Kar, A.A. (2005) Morphological Characterization of White Head Cabbage (Brassica oleracea var. capitata subvar. alba) Genotypes in Turkey. New Zealand Journal of Crop and Horticultural Science, 33, 333-341.
    https://doi.org/10.1080/01140671.2005.9514367

  23. 23. Amurrio, J.M., De Ron, A.A. and Zeven, A.C. (1995) Numerical Taxonomy of Iberian Pea Landraces Based on Quantitative and Qualitative Characters. Euphytica, 82, 195-205.
    https://doi.org/10.1007/BF00029561

  24. 24. Dias, J.S., Monteiro, A.A. and Lima, M.B. (1993) Numerical Taxonomy of Portuguese tronchuda Cabbage and Galega Kale Landraces Using Morphological Characters. Euphytica, 69, 51-68.
    https://doi.org/10.1007/BF00021725