African Journal of
Agricultural Research

  • Abbreviation: Afr. J. Agric. Res.
  • Language: English
  • ISSN: 1991-637X
  • DOI: 10.5897/AJAR
  • Start Year: 2006
  • Published Articles: 6865

Full Length Research Paper

Phenotypic traits detect genetic variability in Okra (Abelmoschus esculentus. L. Moench)

Aaron T. Asare
  • Aaron T. Asare
  • Department of Molecular Biology and Biotechnology, School of Biological Sciences, CANS, UCC, Cape Coast, Ghana.
  • Search for this author on:
  • Google Scholar
Elvis Asare-Bediako
  • Elvis Asare-Bediako
  • Department of Crop Science, School of Agriculture, CANS, UCC, Cape Coast, Ghana.
  • Search for this author on:
  • Google Scholar
Faustina Agyarko
  • Faustina Agyarko
  • Department of Crop Science, School of Agriculture, CANS, UCC, Cape Coast, Ghana.
  • Search for this author on:
  • Google Scholar
Kingsley Taah
  • Kingsley Taah
  • Department of Crop Science, School of Agriculture, CANS, UCC, Cape Coast, Ghana.
  • Search for this author on:
  • Google Scholar
Edmund O. Osei
  • Edmund O. Osei
  • Plant Genetic Resources Institute (PGRRI), Bunso, Ghana.
  • Search for this author on:
  • Google Scholar


  •  Received: 28 April 2016
  •  Accepted: 16 June 2016
  •  Published: 18 August 2016

 ABSTRACT

There is low production of okra in Ghana due to lack of improved varieties and biotic constraints. This study was conducted to characterize okra genotypes to predict genetic variation in the crop. Field trial was conducted to determine genetic variability in 21 okra genotypes. The experiment was based on the randomized complete block design (RCBD) involving planting distance of 0.6 × 0.6 m. Thirty-one quantitative and qualitative data were used to generate a dendrogram. Variations in leaf shape, leaf rib colour, petiole colour, petal colour, colour of the darkest ridges and stem colour were distinctive among the okra genotypes. The mean plant height, canopy diameter, leaf length and breadth, petiole length, internode length, number of branches, days to 50% flowering and fruit yield differed significantly (p ≤ 0.05) among the 21 okra genotypes. These were discriminated into three clusters in a dendrogram with GH3731 as the most diverse. UCCC1, UCCC2, UCCC3, UCCC4 and UCCC5 appeared genetically similar with low fruit yield but early maturity. However, GH5332 had a significantly (p ≤ 0.05) the highest fruit yield of 11.88 t ha-1 but late maturing. UCCC5 or similar genotypes with early maturity trait can be hybridized with GH5332 to improve the yield and earliness.

Key words: Breeding, germplasm, genotypes, genetic diversity and hybridization.


 INTRODUCTION

Okra (Abelmoschus esculentus L. Moench) of the family Malvaceae, is an important and widely cultivated annual crop in both the tropical and sub-tropical regions of the world (Eshiet and Brisibe, 2015; Ali et al., 2014). It is a vegetable  rich  in  organic  and  inorganic  nutrients   that sustain human health and as feed for animals (Chattopadhyay et al., 2011; Ofoefule, 2001; Rahman et al., 2012; Wamanda, 2007; Siesmonsma and Kouame, 2004; Saifullah and Rabbani, 2009). In Ghana, okra is often consumed in the diet by both children and adults in both rural and urban communities. However, the yield of this crop is low due to lack of improved varieties, biotic and abiotic stresses. Yield potential of 2000 to 3000 kg ha-1 has been reported for Okra (MoFA, 2007), depending on the cultivar, harvesting frequency and period for harvesting (Cudjoe et al., 2005). However, actual yields of okra are usually low and also decreased over the years in Ghana, in spite of its economic importance and health benefits (Asare-Bediako et al., 2014).
 
The development of new varieties with better adaptation and yield potential are crucial for sustainable production of okra. Genetic variation in okra is a necessary requirement to improve the crop. Omonhinmin and Osawaru (2005) reported that high degree of wide morphological variation was found among accessions of okra, especially in West African type. There are numerous cultivars of okra with varied plant height, degree of branching and pigmentation of the various parts, period of maturity, and pod shape and size (AdeOluwa and Kehinde, 2011). In addition, Bisht et al. (1995) observed that pigmentation and pubescence of stem, leaf, pods and seeds were important components of variability in okra germplasm.
 
Various types of okra in Ghana are cultivated in the savannah and forest agro-ecological zones that require assessment of their genotypes. Ahiakpa et al. (2013) has done some morphological characterization of okra in Ghana but lacked collections from the central region and inclusion of exotic genotypes. The current work considered collections from the central region and other regions of Ghana as well as Togo to enhance assessment of fully harness variability in the germplasm for breeding and conservation. Indeed, genetic variation may serve as recipe for controlled hybridization to improve the crop. Assessment of variable phenotypic traits of okra would be useful in predicting genetic diversity towards molecular characterization to establish the genetic structure of okra in Ghana. This will facilitate routine breeding and germplasm conservation of the crop. The objective of the current work was to explore phenotypic characteristics to predict genetic variability among 21 okra genotypes in Ghana and to establish basis for molecular characterization.


 MATERIALS AND METHODS

Experimental site
 
The experiment was conducted at the Teaching and Research Farm of the School of Agriculture, at University of Cape Coast during the 2015 major crop season (June to October). This site is located within the coastal savannah vegetation zone, with Acrisol soil type and is a highly endemic site for viral diseases and flea beetle infestation. The area has a bi-modal rainy season from May to June and August to October with an annual rainfall ranging between 750 and 1000 mm and temperatures ranging between 23.2 and 33.2°C with an annual mean of 27.6°C (Owusu-Sekyere et al., 2011).
 
Plant
 
Twenty-one genotypes of okra (both landraces and improved) were used for the study (Table 1). These comprised of fifteen accessions from Plant Genetic Resource Research Institute (PGRR1) at Bunso (wide collections from regions of Ghana and Togo), four farmer varieties, a land race and an improved variety (Asontem) from the Central Region of Ghana. The local names and sources of 21 okra genotypes are shown in Table 1.
 
Field experiment
 
The randomized complete block design (RCBD) with twenty-one genotypes of okra  was sown in four replications. A total land area of 1344 m2 (84 × 16 m) was ploughed and harrowed to render the soil loose. It was then divided into four blocks and each block was further divided into 21 plots, with each plot measuring 3 × 3 m. A distance of 1.0 m was left as walkway between the blocks and 1 m between the plots. Planting was done in June, 2015. The 21 okra genotypes were sown directly at three seeds per hole at a planting distance of 0.6 × 0.6 m and a planting depth of not more than 0.5 cm. The seedlings were later thinned out leaving two seedlings per hill. Weed control was done as necessary using herbicides and a hoe (manual weeding). NPK fertilizer (15:15:15) was applied at a rate of 250 kg ha-1. Watering was done when necessary using sprinklers.
 
Data collection and analysis
 
The quantitative and qualitative data were collected based on the International Plant Genetic Resource Institute (IPGRI, 1991) okra descriptor list and adopted the procedure by Ahiakpa et al. (2013) and Nwangburuka et al. (2011) with some modifications (Table 2). The quantitative data including plant height, canopy diameter, leaf length, breadth, and stem base diameter at first flowering stage and petiole length, fruit length and fruit girth were obtained using meter rule or tape measure. The number of branches and fruits per plant were counted and fresh weight of matured fruits were determined by electronic balance (Radwag, WPT 12C1, Poland). The qualitative parameters were determined by visual estimation and rated (Table 1).
 
Data on plant growth and yield parameters were subjected to one-way analysis of variance (ANOVA) to determine significant differences among the 21 okra genotypes. The means were separated by the least significant difference method, using GenStat Discovery version 4 (VSN International). The GenStat or Minitab 15 statistical software was used for all computations. The Pearson’s correlation coefficients were estimated for the growth and yield parameters. The Eigen-vectors and factor scores were used respectively to measure the relative discriminative power of the PC-axes and their associated characters. The data involving 31 parameters (Table 2) were analyzed with the PowerMarker version 3.5 and the dendrogram generated in Molecular Evolutionary Genetics Analysis version 4 (MEGA4) to determine genetic relatedness among the okra genotypes (Tamura et al. 2007).
 
 
 


 RESULTS AND DISCUSSION

Qualitative and quantitative characteristic variations exist among the 21 okra genotypes. Variations in matured leaf colour, leaf shape, leaf rib colour, petiole colour, petal colour, colour of the darkest ridges and stem colour were distinctive differentiation characters. In addition, differences in flowering span, fruit colour, fruit shape and number of ridges per fruit, pubescence, position of fruits on main stem and branching position at main stem were evident among the okra genotypes (Table 3).
 
The mean plant height, canopy diameter, leaf length and breadth, petiole length, internode length and number of branches as well as days to first flowering, 50% flowering and fruit yield differed significantly (p < 0.05) among the 21 okra genotypes. However, the stem diameter did not show significant (p > 0.05) variation among the okra genotypes. The highest average plant height of 68.95 cm was noted for UCCC3 and the least plant height of 28.72 cm was for GH4374 (Table 4). Plant height at flowering and fruiting are of particular interest for breeding programmes, because the presence of plants with tall and thin stems will increase the rate of lodging near harvesting and this could lead to loss of dry matter and subsequent decrease in fruit yield (Esthiet and Brisibe, 2015). In fact, Verma (1993), Ariyo et al. (1987) and Perdosa (1983) intimated that plant height is controlled by genetic factors and is closely associated with number of flowering node, average fruits per plant and number of internodes.
 
 
 

The leaves of okra serve as the main sites for photosynthesis, an increase or a decrease in their size could affect production of assimilates in the crop. Larger size leaves in any okra genotype may have higher  ability to intercept solar radiation to assume higher photosynthetic capacity, which may enhance growth and crop yield. GH3734 had the widest canopy diameter of 100.86 cm compared to the least of UCCC5 (59.14 cm). The highest leaf length of 21.82 cm was produced by GH3734 and the lowest of 14.84 cm was associated with UCCC5. The mean leaf breadth of 29.55 cm was the highest observed for GH3760 compared to  the least average leaf breadth of 19.89 cm for UCCC5  (Table 4). According to Ahiakpa (2013), an increased leaf area index and a resultant higher fraction of intercepted radiation and its utilization efficiency may increase crop yield. Significant (p < 0.05) correlations were observed between leaf length and canopy diameter (r = 0.72), breadth and canopy diameter (r = 0.53), leaf length and breadth (r = 0.65), petiole length and canopy diameter (r = 0.49) and stem diameter and canopy diameter (r = 0.70) in the okra germplasm could be determinants for plant vigour and yield indicators. GH2026 had the longest internodes of 21.53 cm per plant compared with the shortest for UCCC4 (14.69 cm per plant).

The mean petiole length of 20.64 cm was the highest for GH2052 compared with that of the lowest average petiole length of 14.29 cm for UCCC5. GH2026 had highest mean internodes (21.53 cm) and number of branches (5), respectively. However, GH2057, GH2063, GH4372 and GH5793 had the least number of 2 branches per plant. Variations in petiole length, leaf  size, canopy diameter, number of branches and stem diameter may have implications for crop yield and stability to control lodging. 

According to Ariyo and Odulaja (1991), variability in okra germplasm is more prominent in days to flowering, plant height and various fruit characteristics and these traits could be important in differentiating varieties of A. esculentus. Similarly, in the current study, fruit length, girth and weight as well as the days to first flowering and days to 50% flowering differ significantly (p < 0.05) among the 21 genotypes of okra (Table 5). However, UCCC1, UCCC2, UCCC3, UCCC4, and UCCC5 were very similar in early flowering and days to 50% maturity as well as average fruit number, weight and size. Generally, all okra genotypes with high vegetative growth delayed flowering and maturity. UCCC5 was first to flower at 47 days and attained 50% flowering at 53 days, respectively after sowing seeds, which were significantly early compared to others. On the contrary  UCCC6 was very late to first flower at 139 days and 50% flowering at 141 days, which were significantly (p < 0.05) different from all the other okra genotypes.

In this study, the significantly (p < 0.05) high yielding okra genotype, GH5332, produced 20 fruits per plant, with the highest fruit weight of 11.88 t ha-1, which compared well with the size of the fruits (mean fruit length of 14.2 cm, girth 20 mm, and weight of fruit per plant of 21.6 g) among the okra genotypes. However, GH5332 is late maturing with 50% flowering at 101 days. Indeed, Esthiet and Brisibe (2015) reported that  fruit  length,  pod number and pod weight are the most important determinants of production or yield in okra. It has been suggested that the number of days and plant height at flowering are controlled by the same genetic variables (Choudhary et al., 2006; Hussain et al., 2006).

It is critical to consider early maturity in the phase of erratic rainfall as essential trait to complement  yield for hybridization to produce climate-smart okra genotypes. Therefore, UCCC5 with very early flowering and fruiting traits can be hybridized with the high yielding but late maturing genotypes of GH5332 to improve the crop. A successful cross between unrelated varieties may result into an array of elite genotypes from which advantageous agronomic line may be selected (Ali et al., 2014).

The variation in the quantitative  characteristics which accounted for the total variance includes number of fruits per plant, mean plant height, canopy diameter, leaf length, breadth, stem diameter, petiole length, internode length, and number of branches. The proportion contributed by each quantitative variable to determine the total variation within each Principal Component (PC) axis is shown in Table 6. The variations in the quantitative characters contributed significantly (Eigen vector ≥ 0.2) to the variation within each of the four PC-axes as 38.00, 14.90, 12.40 and 11.70% for PC1, PC2, PC3 and PC4, respectively. The cumulative proportion of variation explained by the first four PC-axes, 77.00% (Table 6)  compared well with observations made by Campos et al. (2005) and Ogunbayo et al. (2005) that the PC-axes contributed 76.62 and 64.5% variations, respectively. Similarly, Ahiakpa et al. (2013) reported that the first four PC-axis contributed 82.97% of the variations in okra. The remaining six axes in the current study accounted for only 23.00% of the total variation. Indeed, canopy diameter, leaf length, breadth, stem diameter and petiole length contributed to the variation in PC1. Plant height, length of internodes and number of branches accounted for the variations observed in PC2 and fruit per plant as well as plant height contributed to the variations in PC4. These variations may suggest the existence of genetic diversity in okra that can be harnessed to improve the crop. Similar observation was made by Yonas et al. (2014). 

 

The 21 okra genotypes were distinguished into 3 main clusters (I, II and III) in the dendrogram at 43% genetic dissimilarity based on 31 quantitative and qualitative morphological characters (Figure 1). All the clusters were made up of varied sub-clusters with the exception of Cluster II, which had a single okra genotype (4.8%) involving GH3731, but the most diverse of all. However, cluster III made of 33.3% of the 21 okra genotypes appeared more closely related, including GH3734, GH6211, UCCC1, UCCC2, UCCC3, UCCC4 and UCCC5, which may suggest genetic similarity. The remaining 61.9% of the okra genotypes were distinguished in cluster I, which is the largest and made up of all the okra genotypes collected from the national gene bank, the Plant Genetic Resources Research Institute at Bunso with the exception of UCCC6. At 25% genetic dis-similarity, all the 21 okra genotypes were fully distinguished in the dendrogram. The dendrogram generated from genetic distance matrices gave an overall pattern of variations and relatedness among the okra genotypes, which agreed with the observation made by Nwangburuka et al. (2011).
 
Indeed, the dendrogram offered distinctive synopsis of the genetic relatedness in the okra germplasm, which is in agreement with the observation made by Aliyu and Fawole (2001) as well as Aremu et al. (2007). According to Ahiakpa et al. (2013), there is a direct relation between the eco-geographical origins of okra collections and  their clustering patterns in Ghana. Similarly, in the current work, the okra collections from the central region of Ghana clustered together and differentiated from those obtained from the PGRRI of Bunso. In addition, the germplasm  collections  from  Togo  especially   GH2057, GH2052 and GH2063 being stored in PGRRI may be genetically similar for their close relatedness in cluster I. The distribution of okra genotypes might be influenced by farmer-consumer preferences, as well as okra trade and germplasm collection activities.

 
 

 


 CONCLUSIONS

The 31 quantitative and qualitative characters distinguished all the 21 okra genotypes without identifying clones. The discriminatory ability of the 31 characters was evident in clustering of the 21 okra genotypes in the dendrogram. UCCC1, UCCC2, UCCC3, UCCC4, UCCC5, GH6211 and GH3734 appeared more closely related. On the whole, the most diverse okra genotype was GH3731. UCCC5 had 50% flowering at 53 days which suggests very early maturity, followed by UCCC3 (68 days), UCCC2 (70 days). Though, GH5332 is late maturing with 50% flowering at 101 days, it produced a significantly (p < 0.05) highest fruit yield of 11.88 t/ha at a rate of 20 fruit per plant which compared well with the size of the fruits among the okra genotypes. UCCC5, UCCC3 and UCCC2 with early maturing but low yield can be hybridized with the high fruit producing, but late maturing okra genotypes of GH5332 and GH6105 to improve earliness in fruiting and adapt the crop to escape terminal drought. The most diverse genotype GH3731 could also be incorporated into breeding to broaden the genetic base of the crop. The almost 50% of the okra genotypes that produced fruits and also had large leaves suitable for use as leafy vegetables and to feed cattle could serve a dual purpose.


 RECOMMENDATION

Though the phenotypic characters were  useful  to  detect genetic variations in okra germplasm collections, they are not absolutely reliable since the traits can be influenced by the environment. Hence, there is a need to employ molecular markers to characterize the okra germplasm including exotic genotypes to establish the genetic structure of the crop and establish baseline information for breeding and conservation of the crop.


 CONFLICT OF INTERESTS

The authors have not declared any conflict of interests.


 ACKNOWLEDGEMENT

The authors wish to acknowledge the University of Cape Coast in Ghana for sponsorship to conduct this initial genetic diversity study in okra.



 REFERENCES

AdeOluwa OO, Kehinde OB (2011). Genetic variability studies in West African okra (Abelmoschus caillei). Agric. Biol. J. N. Am. 2(10):1326-1335.
Crossref
 

Ahiakpa JK, Kaledzi PD, Adi EB, Peprah S, Dapaah HK (2013). Genetic diversity, correlation and path analyses of okra (Abelmoschus spp. (L.) Moench) germplasm collected in Ghana. Int. J. Dev. Sust. 2(2):1396-1415.

  
 

Ali AS, Shah H, Gul R, Ahmad H, Nangyal H, Sherwani KS (2014). Morpho-Agronomic Characterization of Okra (Abelmuscus esculentus L.). World App. Sc. J. 31(3):336-340.

  
 

Aliyu B, Fawole I (2001). Inheritance of pubescence in crosses between V. unguiculata and V. rhomboidea. Nig. J. Genet. 15:9-14.

  
 

Aremu CO, Adebayo MA, Ariyo OJ, Adewale BB (2007). Classification of Genetic diversity and choice of parents for hybridization in cowpea Vigna unguiculata (L.) Walp for humid savanna ecology. Afr. J. Biotechnol. 6(20): 2333-2339.
Crossref

  
 

Ariyo OJ, Akinova ME, Fatokun PE (1987). Plant character correlation and path analysis of pod yield in okra. Euphytica 36:677-686.
Crossref

  
 

Ariyo OJ, Odulaja A (1991). Numerical analysis of variation among accessions of okra. (A. esculentus [L.] Moench). Malvaceae. Ann. Bot. 67:527-531.

  
 

Asare-Bediako E, Addo-Quaye AA, Bi-Kusi A (2014a). Comparative efficacy of phytopesticides in the management of Podagrica spp and mosaic disease on okra (Abelmoschus esculentus L.). Am. J. Expt. Agric. 4(8):879-889.
Crossref

  
 

Bisht IS, Mahajan RK, Rana RS (1995). Genetic diversity in South Asian okra (Abelmoschus esculentus) germplasm collection. Ann. Appl. Biol. 126:239-550.
Crossref

  
 

Campos ET, Espinosa MAG, Warburton ML, Monter AV (2005). Characterisation of mandarin (Citrus spp) using morphological and AFLP markers", Interciencia 30(11):1-14.

  
 

Chattopadhyay A, Dutta S, Chatterjee S (2011). Seed yield and quality of okra as influenced by sowing dates. Afr. J. Biotechnol. 10:5461-5467.

  
 

Choudhary UN, Khanvilkar MH, Desai SD, Prabhudesai SS, Choudhary PR (2006). Performance of different okra hybrids under North Konkan coastal zone of Maharashtra. J. Soils Crops 16:375-378.

  
 

Cudjoe AR, Kyofa-Boamah M, Nkansah GO, Braun M, Owusu S, Adams E, Monney E, Attasi R, Owusu P, Sarpong S (2005). Commercial Okra Production in Ghana - Good Agricultural Practices/Code of Practice and IPM Strategies: In: Kyofa-Boamah M, Blay E, Braun M, Kuehn A (eds). Handbook of Crop Protection. Recommendations in Ghana, Ministry of Food and Agric. Accra pp. 75-92.

  
 

Esthiet JA, Brisibe AE (2015). Morphological Characterization and Yield Traits Analysis in Some Selected Varieties of Okra (Abelmoschus Esculentus L. Moench). Adv. Crop. Sci. Technol. 3(5):1-5.

  
 

Hussain S, Muhammad NS, Shah A, Iqbal Z (2006). Response of okra (Abelmoschus esculentus) cultivars to different sowing times. J. Agric. Biol. Sci. 01:55-59.

  
 

IPGRI (1991). Okra descriptor, Diversity for Development. Int. Plant Gen. Res Inst. Rome.

  
 

Nwangburuka CC, Kehinde OB, Ojo DK, Denton OA, Popoola AR (2011). Morphological classification of genetic diversity in cultivated okra, Abelmoschus esculentus (L) Moench using principal component analysis (PCA) and single linkage cluster analysis (SLCA). Afr. J. Biotechnol. 10(54):11165-11172.
Crossref

  
 

Ogunbayo SA, Ojo DK, Guei R, Oyelakin OO, Sanni KA (2005). Phylogenetic diversity and relationships among 40 rice accessions using morphological and RAPDs techniques. Afr. J. Biotechnol. 4(11):1234-1244.

  
 

Ofoefule SI, Chukwu AN, Anayakoha A, Ebebe IM (2001). Application of Abelmoschus esculentus in solid dosage forms 1: use as binder for poorly water soluble drug. Indian J. Pharm. Sci. 63:234-238.

  
 

Omonhinmin CA, Osawaru ME (2005). Morphological characterization of two species of Abelmoschus: Abelmoschus esculentus and Abelmoschus caillei. Genet. Resour. Newsl. 144:51-55.

  
 

Owusu-Sekyere JD, Alhassan M, Nyarko BK (2011). Assessment of climate shift and crop yields in the Cape Coast area in the Central Region of Ghana. ARPN. J. Agric. Biol. Sci. 6(2):49-54.

  
 

Perdosa J (1983). The morphological characters of okra introductions. Horticult. Bras. 1(1):14-23.

  
 

Rahman K, Waseem M, Kashif MS, Jilani M, Kiran G (2012) Performance of different okra (Abelmoschus esculentus L.) cultivars under the agro-climatic conditions of Defra Ismail Khan. Pak. J. Sci. 64:316-319.

  
 

Saifullah M, Rabbani MG (2009). Evaluation and Characterization of Okra (Abelmoschus esculentus L. Moench.) Genotypes. SAARC J. Agric. 7(1):92-99.

  
 

Siesmonsma JS, Kouame C (2004). Vegetables. In: Plant Resources of Tropical Africa 2 (Grubben, G.J.H. & Denton, O.A., Eds.) PROTA Foundation, Wageningen, Netherlands and Backhuys Publishers, Leinden, Netherlands and CTA, Wageningen, Netherlands pp. 20-29.

  
 

Tamura K, Dudley J, Nei M, Kumar S (2007). MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. (Accessed on 28th March, 2016). Mol. Biol. Evol. 24:1596-1599. 
Crossref

  
 

Verma VD (1993). Collecting eggplant and okra in Madhya Pradesh and Maharashtra. IBPGR. Newsletter for Asia Pacific and Oceania 13:14-15.

  
 

Wamanda DT (2007). Inheritance studies in collected local okra (Abelmoschus esculentus L. Moench) cultivars. In: Combining ability analysis and heterosis on diallel cross of okra. Afr. J. Agric. Res. 5:2108-2155.

  
 

Yonas M, Garedew W, Debela A (2014). Maltivariate analysis among Okra (Abelmoschus esculentus L. Moench) Collections in South Western Ethiopia. J. Plant Sci. 9(2):43-50.
Crossref

  

 




Copyright © 2024 Author(s) retain the copyright of this article.

This article is published under the terms of the Creative Commons Attribution License 4.0

Back to Vol. 11 No. 33
Back to articles
SUBSCRIBE TO RSS
SUBMIT MANUSCRIPT
CONFERENCES
          */?>