Journal of
Horticulture and Forestry

  • Abbreviation: J. Hortic. For.
  • Language: English
  • ISSN: 2006-9782
  • DOI: 10.5897/JHF
  • Start Year: 2009
  • Published Articles: 281

Full Length Research Paper

Morphological characterization and Simple Sequence Repeats (SSRs) based DNA fingerprinting of selected mango (Mangifera indica L.) genotypes in Bangladesh

Md. Rezwan Molla
  • Md. Rezwan Molla
  • Plant Genetic Resources Centre, Bangladesh Agricultural Research Institute, Gazipur, Bangladesh.
  • Google Scholar
Iftekhar Ahmed
  • Iftekhar Ahmed
  • Plant Genetic Resources Centre, Bangladesh Agricultural Research Institute, Gazipur, Bangladesh.
  • Google Scholar
Md. Amjad Hossain
  • Md. Amjad Hossain
  • Plant Genetic Resources Centre, Bangladesh Agricultural Research Institute, Gazipur, Bangladesh.
  • Google Scholar
Md. Shafiqul Islam
  • Md. Shafiqul Islam
  • Regional Horticulture Research Station, Bangladesh Agricultural Research Institute, Chapainawabganj, Bangladesh.
  • Google Scholar
Md. Aziz Zilani Chowdhury
  • Md. Aziz Zilani Chowdhury
  • Crops Division, Bangladesh Agricultural Research Council, Farmgate, Dhaka, Bangladesh.
  • Google Scholar
Dilruba Shabnam
  • Dilruba Shabnam
  • Department of Agricultural Extension, Plant Quarantine Wing, Narayanganj, Bangladesh.
  • Google Scholar
Md. Motiar Rohman
  • Md. Motiar Rohman
  • Plant Breeding Division, Bangladesh Agricultural Research Institute, Gazipur, Bangladesh.
  • Google Scholar

  •  Received: 18 July 2019
  •  Accepted: 09 September 2019
  •  Published: 31 October 2019


Nineteen genotypes of mango including nine released varieties viz. BARI Aam-1, BARI Aam-2 (Laxmanbhog), BARI Aam-3, BARI Aam- 4 (Hybrid), BARI Aam-5, BARI Aam-6, BARI Aam-7, BARI Aam-8, BARI Aam-9; one parental line viz. M- 3896 and nine Geographical Indication Crops (GIs) viz. Haribhanga, Surjapuri, Fazli, Gourmoti, Ashwina, Khirsapat, Gopalbhog, Langra and Ranipasand were characterized with a view to identifying the degree of morphological and molecular variation of mango within genotypes with their historical background their historical background, and to establish a permanent database for documentation of mango in Bangladesh. Wide variations were observed among GI crops and released varieties included in this study for plant, leaf, flower and fruit characters. Among 19 mango genotypes, eight were distinct by two traits and 11 by only single character. Molecular characterization was carried out with SSR markers. Using 21 primers across 19 genotypes a total of 80 alleles with an average number of 3.81 alleles per locus were found of which MIAC-6 and MIAC-11 showed the highest number of alleles (6) (size ranging from 244 to 312 and 133 to 167 bp, respectively). However, the lowest number of allele (2) with size ranging 237 to 366 and 118 to 125 bp was observed in the locus MiSHRS-39 and MIAC-11, respectively. The polymorphic information content (PIC values) ranged from 0.349 to 0.781, with a mean value of 0.602 for all loci. Of the 21 SSR primers, 13 were highly informative (PIC value ≥0.6). The distinct level of heterozygosity indicated higher level of diversity among the genotypes. Band patterns corresponding to individual genotype have been identified to discriminate the genotype. The genotypes presented genetic distances between 0.260 and 1.557. The dendrogram generated from UPGMA cluster analysis broadly placed 19 mango genotypes into two major groups, “A” and “B” in which only one poly-embryonic genotype namely BARI Aam-8 congregated in a distinct group “B” and other 18 mono-embryonic genotypes clustered in group “A”. The dendrogram revealed that Gourmoti and Ashwina were the most similar hybrids with 21% similarity. Contrary to this, hybrids BARI Aam-5 and BARI Aam-8 were the most divergent with a diversity value of 1.56.

Key words: Geographical Indication Crops (GIs), historical background, morphology, identity.



Mango (Mangifera indica L. 2n=40)  is  a  member  of  the family  Anacardiaceae in the order Sapindales, a family of mainly tropical species with a few representatives in temperate regions (Viruel, 2005). The mango is considered as one of the oldest cultivated trees in the world. Historical records provide conflicting accounts for origin and distribution of mango. Although some authors have considered India as the centre of origin due to the high degree of mango diversity observed in that country (Ravishankar et al., 2000), taxonomic and molecular evidence also supports an evolution of mango within a larger area including northwestern Myanmar, Bangladesh and Northeastern India (Mukherjee 1997). Total world mango production is 26 million tons, and it is one of the most important fruits in the world, along with bananas, oranges, grapes and apples (Ukoskit, 2007). Its popularity and importance can easily be realized by the fact that it is often referred as the “the king of fruits” in the tropical world (Purseglove, 1972). The King among fruit is thriving very well in Bangladesh. Mangoes grow widely all over Bangladesh and there are innumerable varieties to charm the connoisseur. Each variety has its own admirers. Each has its distinctive flavour and arguments about the superiority of one over the other can get very serious. Though these are table varieties, meant to be relished as cut fruit, there are others that are used for making jam, jelly, squash, chutney, and pickle. The raw green mango is even added to dal or curry to enhance the flavour. According to the latest statistics provided by BBS (2017), indicated that, the production of mango in Bangladesh is1288000 metric tons. Mango contributes 21.77% to total fruit production in Bangladesh. Several local and exotic cultivars are grown in the country. The cultivars are mostly location specific. Bangladesh Agricultural Research Institute (BARI) has developed eleven improved varieties. Some cultivars have been originated in some localities and are being cultivated in those areas from more than several hundred years, which can be termed as geographical indication crops (GIs) for those localities. It is essential to characterize the GIs and released varieties of mangoes both in morphological and molecular level for establishment of Intellectual Property Right (IPR). Moreover, selection and correct identification of genotypes is essential for any breeding and improvement effort, is difficult, inefficient and inaccurate when based on morphological traits only. Even though a high number of descriptors are used (Thomas et al., 1994), this is due to some phenotypic traits are difficult to describe, and phenotypic data may be influenced by environmental factors and growing conditions, in addition to quantitative inheritance, or partial and complete dominance often confound the expression of genetic traits. Recently, as in other fruit tree species molecular identification of mango cultivars has been  carried   out  with  different  molecular  systems  as isozymes, minisatellites (Jintanawongse and Changatragoon, 2000), RFLPs (Chunwongse et al., 2000; Capote et al., 2003; Ravishankar et al., 2004), AFLPs (Kashkush et al., 2001; Yamanaka et al., 2006), ISSRs (Singh et al., 2007; Bajpai et al., 2008; Samant et al., 2010).  While DNA profiles based on polymorphic band patterns from Random Amplified polymorphic DNA (RAPD) analysis have been described for several fruit species including mango (Rahman et al., 2007; Abirami et al., 2008; Bajpai et al., 2008; Pruthvish and Chikkaswamy, 2016). But these markers have some limitation like RFLPs requires the use of radioactivity and is labour intensive. RAPDs and AFLPs identify only dominant alleles and are sensitive to PCR amplification. Different thermocyclers, Taq polymerases, DNA primer concentrations and even the skill of the experimenter can influence the results of RAPD marker (Sefc et al., 2001). Considering these aspects, a project was undertaken by Agricultural Research Council (BARC) to characterize GIs and released varieties of some crops grown in the country with the financial support of Sponsored Public Goods Research-National Agricultural Technology Project (SPGR-NATP) Phase-1. Bangladesh Agricultural Research Institute (BARI) has been assigned with 10 of its mandated crops and mango was one of them. Nineteen mango genotypes (one advanced line, nine released varieties and nine GIs) have been included in this programme. The present study was, therefore, undertaken to identify distinct morphological characteristics along with establish allelic patterns and estimate ge­netic distances based on microsatellite markers for 19 mango genotypes to generate a refer­ence database to support cultivar protection and settle possible commercial disputes as well as to guide breeding programmes and genetic resources of the species.



Selection of trees for GI released varieties

Nineteen genotypes of mango including nine released varieties viz. BARI Aam-1, BARI Aam-2 (Laxambogh), BARI Aam-3, BARI Aam-4 (Hybrid), BARI Aam-5, BARI Aam-6, BARI Aam-7, BARI Aam-8, BARI Aam-9; one parental line viz. M-3896 (Male parent of BARI Aam-4) and nine GIs viz. Haribhanga, Surjapuri, Fazli, Gourmoti, Ashwina (female parent of BARI Aam-4), Khirsapat, Gopalbhog, Langra and Ranipasand were characterized both at morphological and molecular level. Historical background was not recorded in two genotypes viz Gourmoti and Ranipasand. Morphological characterization of seventeen genotypes was done from standing trees. The centre of diversity of most concentrated with experienced fruit of the respective GI was identified through discussion with experienced fruit scientists and Department of Agriculture Extension (DAE) Officials at district and upazila level. A team of scientists, involved in this programme, visited particular geographical locations and located the targeted trees of selected GIs. Then three plants were labeled with laminated paper sheet as plant number 1, 2 and 3 for each GI. For released varieties and parental line, the team visited Regional Horticulture Research Station (RHRS), Chapainawabganj and Fruit and conserved original plant. Discussing with the station heads and working scientific personnel, the original mother tree (s) (OMT) was identified and selected for data collection. In cases where there was only one OMT, daughter mother trees (DMT) were also selected for data collection.  Each OMT and DMT was labeled as plant number 1, 2 and 3. The list of mango genotypes used in this study with its  data  collecting  site  is given in Table 1.



Recording historical background

Scientists discussed with aged people of the growing areas to find out the historical background of the respective GIs of mango. The team also located some very old trees (≥100 years) as indicated by the nearby people and symptoms on the tree like canopy coverage, trunk circumference, extra-rough trunk surface and galls on trunk etc. Typical plant photographs of Ashwina and Surjapuri were shown in Figures 1 and 2.




Management practices

Farmers normally sell their crop as total plantation either in orchard or in homesteads just after harvest of the previous year’s fruit. Then entire liability of the plantation goes to the traders. Management practices like pruning, weeding, irrigation and fertilizer and pesticide application etc. were done by the traders. Chemical fertilizers like urea, triple super phosphate (TSP), muriate of potash (MoP), gypsum etc. were applied at different rates. None of the traders used recommended doses of fertilizers. Pesticides as prescribed by the dealers or experienced traders were used indiscriminately even one or two days before harvesting especially for controlling fruit fly. In the research stations application of fertilizers and other cultural practices like ploughing, weeding, irrigation, pruning etc. were done as per recommendation of BARI.

Observation, data collection and record keeping

The selected trees were visited frequently at different stages of growth, flowering and fruiting. Passport information and morphological  data  in  respect  of  plant, leaf, flower/inflorescence,  fruit and stone characters were recorded following IPGRI descriptors for mango (IPGRI, 2006). The photographs of the specific trait considered to be helpful for identification of the variety/cultivar were taken from each genotype at appropriate time of traits to compare the distinctness among them. Data related to distinctness in morphological traits were photographed on each of the 19 mango genotypes.

Studies on molecular characterization

Extraction of genomic DNA

Young and mature leaf samples were collected from the particular plants which were used for morphological characterization at the particular geographical locations and the research station from where the variety was released. The genomic DNA was isolated from leaf tissues following a standard protocol described by Uddin et al. (2014) with minor modifications. Essentially, the extraction buffer-1 composition was [0.4 M glucose, 20 mM ethylenediamine tetra acetic acid (EDTA; pH 8.0), 3% (w/v) polyvinyl pyrrolidone (PVP)-40 (molecular weight 40,000) and 0.2% (v/v) β-mercaptoethanol]. Preheated (65°C) solution-2 [2% cetyl trimethyl ammonium bromide (CTAB) (w/v), 100 mM Tris (pH 8.0), 20 mM EDTA (pH 8.0), 1.4 M NaCl, and 0.15% (v/v) β- Mercaptoethanol] was added as extraction buffer-2, 0.15% (v/v) β-mercaptoethanol was added and the mixture was mixed gently and incubated at 65°C in a water bath for 1 h with intermittent shaking. DNA was precipitated with ice-cold and extra pure isopropyl alcohol and purification with absolute ethanol (Plus sodium acetate, 3 M) and 70% ethanol chronologically. DNA sample of each mango germplasm was dissolved in 50 μl of TE buffer. When the DNA pellet was totally dissolved in TE buffer, 4 μl RNase (10 mg/ml) was added to isolated DNA and incubated at 37°C for 1.5 h. Finally, DNA sample was stored at -20°C.

Quantification and optimization of DNA concentration

The presence of genomic DNA was confirmed on 1% agarose gel qualitatively. The gels were visualized under UV light and photographed using photo documentation system (UV Transilluminator, Uvitec, UK). All of the DNA samples were found to be in good quality in this study. The amount of genomic DNA was quantified using UV a spectrophotometer (Spectronic® GENESYS™  10 Bio) at 260 nm. Using the absorbance reading obtained for DNA sample of each mango genotypes, the original DNA concentrations were determined.

Selection of microsatellite primers

Twenty five SSR primer pairs described previously in the literature (Schnell et al., 2005, Duval et al., 2005, Viruel et al., 2005, Kittipat, 2007; Wahdan et al., 2011) were used for microsatellite analysis in the present study. Among the 25 primers pairs, 21 except MiSHRS-18, mMiCIR022, mMiCIR029 and mMiCIR020 (Table 2) showed better responsiveness with clear and expected amplified product sizes.



PCR standardization and amplification

Microsatellites amplification was performed in 10-μL volume containing 5X Green GoTaq® Reaction Buffer (Promega, USA) 15 mM MgCl2, 1.25 U Taq DNA polymerase (Thermo Scientific, USA), 0.4 mM each of the dNTPs (NEB, USA), 10 μM forward and reverse primers and 50 ng template DNA. The mixtures were prepared at 0°C and transferred to the thermal cycler. Amplification reactions of SSR loci were carried out in a Mastercycler® nexus Gradient thermal cycler (Eppendorf, Germany). The PCR profile included initial denaturation for 4 min at 94°C, followed by 32 cycles of denaturation at 94°C for 45 s, annealing at 46 to 55°C for 1 min, extension at 72°C for 1 min, and final extension at 72°C for 8 min. After completion of cycling program, reactions were held at 10°C. For checking amplification, the amplified products were resolved using on 2% agarose gel containing Ethidium bromide in electrophoresis chamber. If the primer was shown good band resolution intensity, less smearing, amplifying the target genomic region of template DNA, the PCR protocol was considered to be correct.

Electrophoretic separation and visualization of PCR products

After standardization of PCR, amplified products of 19 mango genotypes against each primer were electrophoresed on a 5% denaturing polyacrylamide gel containing 19:1 acrylamide: bis-acrylamide, 10X TBE buffer, 10% APS and ultrapure Temed. Electrophoresis was done using the Triple Wide Mini-Vertical Electrophoresis System, MGV-202-33 (CBS Scientific, USA). The gel is run at 80 to 90V and 20°C temperature maintained by a cooling system (Julabo, Germany) upon loading of PCR products for a specified period of time depending on the size of amplified DNA fragment (usually 1 h for 100 bp). After completion of electrophoresis, the gel was stained with Ethidium bromide and the individual bands were scored for analysis.

Scoring and analysis of microsatellite data

SSR markers were scored as codominant, so homozygous and heterozygous genotypes could be distinguished in individual plants. The bands representing particular alleles at the microsatellite loci were scored manually and designated the bands as A, B, C, etc. from the top to the bottom of the gel. The genotypes of different individuals were hypothetically scored as AA, BB, CC, etc. for homozygous or as AB, AC, BC etc. for heterozygous. A single genotypic data matrix was constructed for all loci.  Statistics of genetic variation (number of observed and effective alleles, Nei’s gene diversity, Shannon’s information index, heterozygosity and polymorphic) were calculated using allelic frequency estimates obtained     from  genotypic  frequencies  of   SSR   loci   using   the computer program POPGENE (Version 1.31) (Yeh et al., 1999). In addition, Chi-square test (1:2:1) for Hardy-Weinberg equilibrium for each population was obtained for SSR alleles using this programme. The microsatellite data matrix was used to calculate Nei’s distance (Nei, 1972), and to generate the corresponding matrix of genetic distance estimates among accessions and cluster analyses were performed on the genetic distance matrix by using UPGMA  method   to   determine  the  relationships  among accessions (dendrograms) using POPGENE (Version 1.31) (Yeh et al., 1999). The polymorphism information content (PIC) of the SSR used or gene diversity value was calculated as PIC= 1- Σf2ij; where fij is the frequency of the ith allele for the jth SSR locus (Anderson et al., 1993). PIC values provided an estimate of the discriminatory power of any locus by considering the number of alleles per locus and the relative frequencies of those alleles in the population.  The  software  DNA  FRAG  version  3.03  was used to estimate allelic length (Nash, 1991).






Historical background of GIs

During  visit  to  Shibganj upazila under the district of Chapainawabganj, Bangladesh, scientists located some very old trees of Fazli (>150 years) and took photographs as evidence. Mr. Md. Emajuddin is the owner of this tree. As the statement of the local people this cultivar was first collected from an old lady Fazli Bibi residing in this upazila. An English Collectorate, Ravenosh has given the name of the cultivar in honour of Fazli Bibi. Similarly, Md. Fariduddin, a man of 80 years old said that he has been seeing Ashwina variety in Bangladesh since his childhood. He also said that he has been informed about this variety from his grandfather late Md Abdul Motaleb. The variety is extremely late and has been commercially cultivating in Chapainawabganj and neighbouring districts of Bangladesh from more than 150 years. Khirsapat mango of Rajshahi region (Northern region of Bangladesh) bears some special quality in respect of taste and flavour, which it might gain from that geographical location. Since fruit pulp of this variety resembles Khirsa (Bengali terminology) hence it is called Khirsa variety of mango. Md. Fariduddin, a man of 80 years old said that he has been seeing this variety in Bangladesh since his childhood. He also said that he has informed about this variety from his grandfather late Md Abdul Motaleb. Scientists located a Langra mango tree having 32 m × 29 m canopy coverage at Namo Chakpara village under Shibganj upazila of Chapainawabganj district. The local people said that the tree is more than 150 years old. Langra mangoes of Chapainawabganj district especially of Shibganj have some special quality in respect of taste and flavour, which it might gain from that geographical location. The name of this variety was given according to a lame man (synonym Langra in Bangla). Md. Fariduddin, a man of 70 years old said that he has been informed about this variety from his grandfather and it has been cultivated in Chapainawabganj from generation to generation. Gopalbhog is a well-known mango variety to everybody and has been commercially cultivated in Chapainawabganj and neighbouring districts. Laxmanbhog (BARI Aam-2) is also a popular mango cultivar commercially cultivated in greater Rajshahi districts but in limited scale all over the country. As per opinion of local people including aged ones, they have been noticed the trees of the cultivar in this area since their childhood. Very old trees located in Chapainawabganj area also indicated its presence in this region from more than 100 years. The cultivar was released as BARI Aam-2 by BARI in 1996 from RHRS, Chapainawabganj, As per description of local people one earthenware vessel maker (Kumar) of Unchubalua village under Mithapukur upazilla of Rangpur district planted a seed of Maldia mango near the garbage of broken earthenware vessel (in Bangla which is termed as ‘haribhanga’) more than 100 years back. According to the statement it can be assumed that the cultivar was originated from chance seedling. Because of excellent taste, the variety got popularity, and people of nearby villages   started   collecting   grafts   of   the  variety,  and spontaneously it was named as ‘Haribhanga’ according to its location. After the division of Indian sub-continent, that Kumar migrated to India, and one slaughter named Suhrab Kosai purchased his house. Late Nofol Uddin Paiker planted one graft of the variety in his homestead about 65 years ago, which is the oldest tree of the cultivar. The variety spread in other areas of the upazilla from this plant.  Md. Abdul Mannan (71 years old villager) described the history of ‘Haribhanga mango’, which was supported by other villagers. It is evident from a very old tree that the Surjapuri cultivar evolved more than 200 years ago.  Md. Nurul Islam son of late Sharif Uddin, Village- Harinmari Nayapara, Union- Harinmari, Upazila-Baliadangi, District- Thakurgaon is the owner of the tree. His fore father Kontu Mohammad (grandfather of his father) planted it. As per opinion of the local people and DAE officials, the variety was evolved from natural cross pollination, and is being cultivated in Baliadangi and Ranishankail Upazilas since several centuries back.

Historical background of released variety and parental line

BARI Aam-1

The origin of this variety is in Chapainawabganj Sadar. It is clear that this variety evolved from chance seedling more than 100 years ago. Its old name is Satiar Kara. Bangladesh Agricultural Research Institute registered this variety as BARI Aam-1 in 1996 from National Seed Board of Bangladesh after a long-time evaluation at Regional Horticulture Research Station, Chapainawabganj. The variety possesses countrywide adaptability, and high export potentiality.

BARI Aam-3

The origin of this variety is in India. It is an Indian hybrid which is known as Amrapali. Bangladesh Agricultural Research Institute registered the variety for cultivation as BARI Aam-3 in 1996 from National Seed Board of Bangladesh after a long time evaluation at Regional Horticulture Research Station, Chapainawabganj and Fruit Research Station, Binodpur, Rajshahi. The variety is commercially cultivated all over the country.

BARI Aam-4

This is the only hybrid variety of Bangladesh. Scientist of RHRS, BARI, Chapainawabganj developed this variety in 1993 from a crossing between Ashwina (a commercial cultivar) and M-3896 (a Florida line). After evaluation, the hybrid was registered as BARI Aam-4 in 2003 from National Seed Board of Bangladesh. It is an outstanding variety. The variety possesses country wide adaptability, and  at  present   commercially   cultivated   all   over   the country. It is a late variety.

BARI Aam-5

The variety is developed from a chance seedling found in the office premises of Regional Agricultural Research Station, Jessore. After evaluation the hybrid was registered as BARI Aam-5 in 2010 from National Seed Board of Bangladesh. It is an early variety. The original mother tree of the variety is about 40 years old, and is still alive. It is an early coloured variety with less juice in the pulp and less sweetness having high export potentiality.

BARI Aam-6

The origin of this variety is in Chapainawabganj, Bangladesh. It was evolved from chance seedling about 100 years ago. Locally the cultivar is known as ‘Bou Bholani’. Bangladesh Agricultural Research Institute organized a National Mango Show at RHRS, Chapainawabganj in 1993. ‘Bou Bholani’ won first prize in ‘Other’ group. RHRS scientists collected the cultivar, and evaluate for more than 15 years. Then it is registered as BARI Aam-6 in 2010 from National Seed Board of Bangladesh. Cultivation of this variety is limited at Chapainawabganj.

BARI Aam-7

The origin of this variety is in Chapainawabganj, Bangladesh. It was evolved from chance seedling 20 years ago. Bangladesh Agricultural Research Institute registered this variety as BARI Aam-7 in 2010 from National Seed Board of Bangladesh. Cultivation of this variety is limited at Chapainawabganj. It is a coloured variety.

BARI Aam-8

The origin of this variety is in Myanmar. It is a Burmese poly-embryonic variety, which is known as Ranguaichi. Bangladesh Agricultural Research Institute registered the variety for cultivation as BARI Aam-8 in 2010 from National Seed Board of Bangladesh after a long time evaluation at different Regions of the country (Chapainawabganj, Jessore, Chittagong, Rangamati Hill district, Khagrachari Hill district). The variety is commercially cultivated in hilly areas of the country. As it has wider adaptability, it can be grown commercially all over the country.

BARI Aam-9

The    origin   of    this   variety   is   in  Chapainawabganj, Bangladesh. It was evolved from chance seedling more 100 years ago. Locally the cultivar is called Tikkafaras. Bangladesh Agricultural Research Institute registered this variety as BARI Aam-9 in 2011 from National Seed Board of Bangladesh after a long time evaluation at the Regional Horticulture Research Station, Chapainawabganj. Cultivation of this variety is limited at Chapainawabganj. It is sweet at green stage (Kanchamitha).


The origin of this line is in Florida, USA. It was introduced to Bangladesh during the 1980s with the financial and technical support of FAO Mango Improvement Project. It is a coloured variety and showed very good performance in respect of yield and quality under Bangladesh condition. Because of high incidence of diseases, the line was not registered for cultivation in the country, and exploited as gene donor. It is the male parent of BARI Aam-4.

Studies on morphological traits

Nineteen varieties/cultivars were used for DUS (Distinctness, Uniformity and Stability) experiment to study the morphological traits and attempt to distinguish one from the other individual on the basis of traits as per IBPGR descriptor. According to the descriptor a total of 102 descriptor traits (15 plants, 20 leaf, 22 flower, 33 fruit, 10 stone and 5 seed descriptor) was used to characterize 19 mango genotypes. Nineteen are important among 102 traits for distinctness of the studied genotypes (Table 3). Among 19 mango genotypes, eight were distinct by two traits and 11 by only single character. The genotypes with respective distinct traits are given in Table 4. Fazli had the maximum fruit weight (780 g) and weak stalk attachment and these two traits distinguish this genotype from the others. Fruiting duration is an important trait which indicates the availability of this genotype. Considering this trait Ashwina was extremely late and its fruiting duration was recorded end of July to end of August which was an identifying character of this cultivar. Similarly, slightly prominent fruit beak discriminated Laxambough from all the genotypes. Pentamerous type flower was found in 16 genotypes while Haribhanga produced both pentamerous and tetramerous type flower which was considered as a unique trait for this cultivar identification. Fruit skin colour of ripe fruit is an important character and three genotypes viz. Surjapuri, BARI Aam-7 and M-3896 showed green with red blush, yellow with red blush and reddish yellow, respectively in ripening stage and these types of skin colour were distinct among the other 14 genotypes. Slightly prominent fruit sinus was a distinguishable character for M-3896. Considering colour   of   young  leaf  all  the  studied  genotypes  were classified into five groups and among these three were comprised of more than one genotype. As for instance, five genotypes (Fazli, Gopalbough, BARI Aam-3, 8 and 9) exhibited reddish brown, eight genotypes (Ashwina, Khirshapat, Laxmanbhog, Haribhanga, BARI Aam-4, 5, 6 and 7) light brick red, two genotypes (Langra and Surjapuri) light green with brownish tinge while the rest two groups contained only one genotype where green with red patches and tinge colour young leaf was observed in BARI Aam-1 and M-3896, respectively which was dissimilar with another one. Again, leaf attitude in relation to branch separated another two genotypes. All the genotypes except BARI Aam-5 and 8 exposed semi-erect leaf attitudes in relation to branch whereas two unique traits, that is, horizontal and semi-drooping attitude were recorded in respect of these two genotypes. In addition, poly-embryonic type has given to BARI Aam-8 for its discrepancy in comparison to another one. BARI Aam-9 showed green pulp colour in ripe fruit, which was divergent from the other genotypes. However, sometimes single character is not sufficient to differentiate one genotype. Therefore, it is effective to combine more than one character to distinguish completely one genotype from others. Such as, light green with brownish tinge colour was found in young leaf of Khirshapat and Surjapuri. When another trait viz. fruit shape or pulp colour of ripe fruit was considered with young leaf colour then this combination of two traits could easily distinguish Khirshapat from Surjapuri. Similarly, fruit shape and pulp colour of ripe fruit showed distinctness for BARI Aam-6. On the other hand, perceptible fruit stalk was observed in BARI Aam-3 and Laxmanbhog, while oblique type was common in BARI Aam-4 and Haribhanga. Hence, BARI Aam-3 and 4 could be identified from other genotypes when another trait like fruit beak type combined with fruit stalk. Again, colour of young leaf and fruit shape showed  distinctness for Langra whereas Gopalbhog was unique for fruit neck prominence and fruit beak type. Gourmoti showed distinctness from the other genotypes based on fruit shape and pulp colour of ripe fruit whereas skin and pulp colour of ripe fruit showed uniqueness for Ranipasand (Figure 3).





Studies on molecular traits

Twenty five (25) SSR primers were used for generating banding profile. Out of which 21 primers (Table 2) were selected in the analysis for their reproducible and polymorphic DNA amplification patterns among genotypes. Two typical SSR profiles are shown in Figure 4. Analysis of the variability parameters for the 21 SSRs in the 19 mango genotypes are shown in Table 5. A total of 80 alleles with an average number of 3.81 alleles per locus were found in the present  study.  The  number  of   alleles  detected  varied form 2 (‘MiSHRS-39’ and ‘MIAC-11’) to 6 (‘MIAC-6’ and ‘MIAC-11’). The allele size ranged from 99 (MIAC-4) to 366 bp (MiSHRS-39). Earlier, Schnell et al. (2005) and Wahdan et al. (2011) reported similar values of SSR polymorphism, number of alleles and allele size in mango cultivars. In the present study, most of the SSR primers detected multiple loci, which can be attributed to the allopolyploid nature of mango (Mukherjee, 1950).




The PIC value provides an estimate of the discriminatory power of a marker by taking into account not only the number of alleles at a locus but also the relative frequencies of these alleles (Huda et al., 2019). All studied SSRs were polymorphic among mango genotypes and informative for describing their genotypic variation (that is, PIC values different from zero). PIC values ranged from 0.349 to 0.781 (Table 5), with a mean PIC of 0.602. Thirteen of these SSRs were very informative (PIC>0.6), with the highest PIC value recorded for mMiCIR014 (0.781) and followed by mMiCIR032 (0.776), MIAC-5 (0.731) and LMMA9 (0.691) which were higher than the average PIC value reported by Wahdan et al. (2011) and Kumar et al. (2013). High PIC values were observed might be due to use of di-nucleotide repeats and also due to genotypic differences (Molla el al., 2010). Indeed, the very informative markers are extremely useful for genetic studies and determination of the level of polymorphism on a specific marker locus (Sundaram et al., 2007).

According to the banding patterns obtained with 21 selected primer pairs, one or two bands  were  present  in each genotype; the amplification pattern seems to indicate the detection of a single locus. Mango has been described as allopolyploid (Mukherjee, 1997) and these results suggest a complete depolarization in this species. The genotypes studied were considered homozygous and heterozygous when one or two fragments were present per locus, respectively (Callen et al., 1993). Consequently, the average observed heterozygosity in each SSR locus considering all studied genotypes were 0.584. Similarly, higher level of heterozygosity (0.587) also found in 19 genotypes considering all SSR locus under study (Table 5) and the 17 out of 19 genotypes showed heterozygosity higher than 0.50 (Table 6). The great heterozygosity can be attributed to the mating system of this species that is normally out cross pollination with some self-pollination. The higher level of heterozygosity observed in the present study has also been reported by Shiran et al. (2007) and Wahdan et al. (2011), probably due to the greater diversity of genotypes used in the present study. For the 19 mango genotypes, a total of 80 alleles were detected using the 21 loci, ranging from 25 for genotype ‘BARI Aam-5’ to 37 for ‘Haribhanga’ (Table 6). The average values for observed number of alleles per genotype were 1.540. The analysis of allelic pattern showed that the number of polymorphic loci within accessions ranged from 4 (19.05) in BARI Aam-5 to 16 (76.19) in Haribhanga and Gopalbhog. The mean Shannon’s information index (I) was 0.409, and ranged from 0.132 in BARI Aam-5 to 0.555 in Gopalbhog (Table 6).



Allele sizing technologies are well established and can be readily used to size microsatellite alleles from any organism (Song et al., 1999). SSR genotypic data from a number of loci have the potential to provide unique allelic profiles or DNA fingerprints for precisely establishing genotypic identity. Distinguishing SSR marker band with their positions are shown in Table 7. The band patterns corresponding to individual genotype may help to recognize the genotype in question. When one primer would not distinguish individual variety from others, another primer should be considered and sometimes combination of more than one primer should be taken into account. Thus, additional primer or set of primers might be needed to test to identify all expected varieties. All the 19 mango genotypes were discriminated successfully by the 21 SSR markers. Among 80 alleles detected, 48 were specific to studied mango genotypes. From a total of 80 scorable alleles, 59 were polymorphic bands and 21 monomorphic. These results indicate that 48 out of the 80 (50%) fragments are considered putative genotypes specific markers for these genotypes. Most of the unique band patterns used for genotype identification was found at locus MIAC-11. Locus MIAC-11 alone discriminated seven genotypes (BARI Aam-3, M-3896, Surjapuri, Gourmoti, Khirsapat, Gopalbhog and Ranipasand). One specific allele was detected in the genotypes BARI Aam-1 (180 bp/mMiCIR032, 134 bp/MiSHRS-4, 125bp/MIAC-11), BARI Aam-2 (159 bp/mMiCIR014), BARI Aam-3 (204 bp/ mMiCIR032, 133 bp/MIAC-5),  BARI   Aam-4  (154  bp/mMiCIR014),  BARI Aam-5 (230 bp/LMMA7), BARI Aam-8 (214 bp/MiSHRS-48, 201 bp/LMMA1, 184 bp/LMMA9), M-3896 (195 bp/mMCIR030), Gourmoti (138 bp/MiSHRS-37, 167 bp/MIAC-5), Ashwina (231 bp/mMiCIR025 and  Khirsapat (207 bp/MiSHRS-48) (Table 7). Besides, BARI Aam-7, BARI Aam-9, Haribhanga, Surjapuri showed unique band with the primer MiSHRS-1 (228/211 bp), mMiCIR032 (198/180 bp), mMiCIR018 (246/215 bp), MIAC-5 (167/133 bp), respectively. Only two genotypes BARI Aam-6 and Fazli could be easily identified in combination of the primer mMiCIR014 and mMiCIR032, in which all primers showed heterozygous condition. Results of the present study represent one of the first attempts to find out a small set of microsatellite makers to discriminate mango genotypes of Bangladesh providing meaningful data that can be enlarged by additional mango genotypes and new microsatellite markers.



The analysis of molecular data showed different levels of genetic diversity among ten mango genotypes determined based on the Nei’s (1972) genetic distance. The genotypes presented genetic distances between 0.260 and 1.557, which reflects the high genetic variability of the collection of mango genotypes studied (Figure 5). Based upon more recent genetic analysis involving microsatellite marker, it is now estimated that monoembryonic cultivar descended from polyembryonic cultivars (Schnell et al., 2005, Viruel et al., 2005). Mango accessions showed genetic differences based on geographical origins and their known history (Viruel et al., 2005, Pandit et al., 2007, Duval et al., 2009, Hirano et al., 2010). On the whole, earlier studies clearly the differentiation of mango accessions regardless the marker system use to fingerprint based on type of embryo (mono- or poly-embryonic), geographical origins, or genetic status (cultivars, landraces, species). In this study, highest genetic distance value (1.557, 1.033, 1.261 and 1.149) were observed between the genotypes BARI Aam-8 vs BARI Aam-5, BARI Aam-2 and BARI Aam-1 (Figure 5). As in type of embryo BARI Aam-8 possess poly-embryonic and others are mono-embryonic. On the other hand, the minimal genetic distance (0.260) was observed between Gourmoti vs Ashwina, both two originated from similar location and popular as late variety.



The dendrogram generated from the unweighted pair group arithmetic average (UPGMA) cluster analysis broadly placed 19 mango genotypes into two major groups “A” and “B” in which only one poly-embryonic genotype namely BARI Aam-8 congregated in a distinct group “B” and other 18 mono-embryonic genotypes clustered in group “A” (Figure 5). Nevertheless, cluster “A” formed two sub-clusters “A1” and “A2”. Sub-cluster “A1”  subsequently   separated   into   another   two   sub-clusters “A1.a”, “A1.b” respectively. Similarly, based on their genetic distances different genotypes grouped into different sub-cluster.  Upon subsequent separation, the highest genetic dissimilarity coefficient was observed between BARI Aam-5 vs BARI Aam-8 (GD=1.557) which is assembled in sub-cluster “A2.a1” and also formed a distinct group “B” respectively. Subgroup “A2.a2” gathered 5 genotypes in which genotypes Gourmoti and Ashwina comprised sharp similarity (GD = 0.260). The genotypes had a distinct status in the dendrogram, because there might have effect of morphological traits and geographical sources. As for instance, BARI Aam-8, BARI Aam-5, BARI Aam-2 and BARI Aam-1 are different based on their morphological features like, fruit shape, shape of fruit apex, fruit cavity, embryo type and so on (Hossain et al., 2014). Moreover, these varieties have been released from different location and also different year.










Historical   background   of   geographical   indications  of mango as described by aged people of their most concentrate areas of cultivation indicated that the cultivars have been originated naturally in those areas. The cultivars possess high commercial value and are being cultivated widely around their areas of origin. GIs possess some special character, which they might gain from their habitat. Characterization of mango genotypes on the basis of DNA fingerprinting data in combination with morphological traits has become an efficient tool to link phenotypic and genotypic variation. Morphological traits and SSRs have been shown to be highly efficient for genotype identification and provide a positive assessment to the ability of SSR marker to produce distinctive DNA profiles of 19 mango genotypes. The results of the present study could be applied as baseline information to maintain the appropriate identity and the construction of a database of all mango cultivars grown in Bangladesh and in broad sense, to protect the mango germplasm of Bangladesh.



The authors have not declared any conflict of interests.



The research results presented here was supported by the institution, Bangladesh Agricultural Research Council, Farmgate, Dhaka for financial support of SPGR-NATP Phase 1 through the “Coordinated Sub-Project on Characterization of Important Plant Genetic Resources: BARI Component”. The technical suggestion provided by Dr. Lutfur Rahman, Retired Professor, Department of Genetics and Plant Breeding, Bangladesh Agricultural University, Mymensingh, is acknowledged with appreciation.



Abirami K, Singh SK, Singh R, Mohapatra T, Kumar AR (2008). Genetic diversity studies on polyembryonic and monoembryonic mango genotypes using molecular markers. Indian Journal of Horticulture 65(3):258-262.


Anderson JA, Churchill GA, Autrique JE, Tanksley SD, Sorrels ME (1993). Optimizing parental selection for genetic linkage maps. Genome 36:181-186.


Bajpai A, Srivastava N, Rajan S, Chandra R (2008). Genetic diversity and discrimination of mango accessions using RAPD and ISSR markers. Indian Journal of Horticulture 65(4):377-382.


Bangladesh Bureau of Statistics (BBS) (2017). Yearbook of Agricultural Statistics-2015, Bangladesh Bureau of Statistics. Ministry of Planning, Government of the People's Republic of Bangladesh, Dhaka pp. 157-158.


Callen DF, Thompson AD, Shen Y, Phillips HA, Richards RI, Mulley JC, Sutherland GR (1993). Incidence and origin of 'null' alleles in the (AC)n microsatellite markers. American Journal of Human Genetics 52:922-927.


Capote M, Cueto J, Becker D, Rohde W (2003). Development and application of various DNA marker types for the characterization of


genetic diversity within commercial mango varieties. Cuba Journal of Genetics and Breeding 57:175-184.


Chunwongse J, Phumichai C, Barbrasert C, Chunwongse C, Sukonsawan S, Boonreungrawd R (2000). Molecular mapping of mango cultivars 'Alphonso' and 'Palmar'. Acta Horticulturae 509:193-206.


Duval ME, Risterucci AM, Calabre C, Le Bellec F, Bunel J, Sitbon C (2009). Genetic diversity of Caribben Mangoes (Mangifera indica L.) using microsatellite markers. Acta Horticulturae 820:183-188.


Duval, MF, Bunel J, Sitbon C, Risterucci AM (2005). Development of microsatellite markers for mango (Mangifera indica L.), Molecular Ecology Notes 5:823-826.


Hirano R, Htun-Oo T, Watanabe KN (2010). Myanmar mango landraces reveal genetic uniqueness over common cultivars from Florida, India, South East Asia. Genome 53:221-330.


Hossain MA, Ahmed I, Molla MR (2014). Characterization of important Plant Genetic Resources: BARI component. Bangladesh Agricultural Research Institute, Joydebpur, Gazipur.


Huda MN, Hasan M, Abdullah HM, Sarker U (2019). Spatial distribution and genetic diversity of wild date palm (Phoenix sylvestris) growing in coastal Bangladesh. Tree Genetics and Genomes 15(1):3.


International Plant Genetic Resources Institute (IPGRI) (2006). Descriptors for Mango (Mangifera indica L.). International Plant Genetic Resources Institute, Rome, Italy.


Jintanawongse S, Changtragoon S (2000). Identification of cultivars and certification of hybrids in mango (Mangifera indica L.) by isozyme markers. Acta Horticulturae 509:177-184.


Kashkush K, Jinggui F, Tomer E, Hillel J, Lavi U (2001). Cultivar identification and genetic map of mango (Mangifera indica ). Euphytica 122:129-136.


Kittipat U (2007). Development of microsatellite markers in mango (Mangifera indica L.) using 5' anchored PCR. Thammasat International Journal of Science and Technology 12:1-7.


Kumar M, Ponnuswami V, Nagarajan P, Jeyakumar P, Senthil N (2013). Molecular characterization of ten mango cultivars using simple sequences repeat (SSR) markers. African Journal of Biotechnology 12(47):6568-6573.


Molla MR, Islam MN, Rohman MM, Rahman L (2010). Microsatellite allele size profiling to determine varietal identity and genetic diversity among groundnut varieties in Bangladesh. Nature and Science 8(12):123-129.


Mukherjee SK (1950). Mango: its allopolyploid nature. Nature 166:196-197.


Mukherjee SK (1997). Introduction: botany and importance. In: Litz RE (ed.). The Mango, Botany, Production and Uses. CAB International, Wallingford, UK.


Nash JHE (1991). DNAfrag, Version 3.03. Institute for biological sciences, National Research Council of Canada, Ottawa, Ontario, Canada.


Nei M (1972). Genetic distance between populations. American Naturalist 106:283 292.


Nei M (1973). Analysis of gene diversity in subdivided populations. Proceedings of the National Academy of Sciences of the United States of America 70:3321-3323.


Pandit SS, Mitra S, Giri AP, Pujari KH, Patil BP, Jambhale ND, Gupta VS (2007). Genetic diversity analysis of mango cultivars using inter simple sequence repeat markers. Current Science 93:1135-1141.


Pruthvish R, Chikkaswamy BK (2016). Genetic diversity and relationships among mango varieties using RAPD molecular markers. International Journal of Current Microbiology and Applied Science 5(1):778-787.


Purseglove JW (1972). Some aspects of mango culture in the Western Tropics. Proceedings of the International Symposium on Mango and Mango Culture, IARI. New Delhi.


Rahman ML, Rabbani MG, Siddique MNA, Rahman MA, Garvey EJ, Rahman EHMS (2007). Molecular Characterization of 28 Mango Germplasm Using RAPD. Plant Tissue Culture and Biotechnology 17(1):71-77.


Ravishankar KV, Anand L, Dinesh MR (2000). Assessment of genetic relatedness among mango cultivars of India using RAPD markers. Journal of Horticultural Science and Biotechnology 75:198-201.


Ravishankar KV, Chandrashekara P, Sreedhara SA, Dinesh MR, Anand


L, Saiprasad GVS (2004). Diverse genetic bases of Indian polyembryonic and monoembryonic mango (Mangifera indica L.) cultivars. Current Science 87:870-871.


Samant D, Singh AK, Srivastav M, Singh NK (2010). Assessment of genetic diversity in mango using Inter-simple Sequence Repeat markers. Indian Journal of Horticulture 67(4):1-8.


Schnell RJ, Olano CT, Quintanilla WE, Meerow AW (2005). Isolation and characterization of 15 microsatellite loci from mango (Mangifera indica L.) and cross-species amplification in closely related taxa. Molecular Ecology Notes 5:625-627.


Sefc KM, Lefort F, Grando MS, Scott KD, Steinkellner H, Thomas MR (2001). Microsatellite markers for grapevine: A state of the art. In: Molecular Biology & Biotecnology of Grapevine. (Roubelakis- Angelakis, KA., Ed.) Kluwer Academic Publishers, Dordrecht, The Netherlands pp. 1-30.


Shiran B, Amirbakhtiar N, Kiani S, Mohammadi SH (2007). Molecular characterization and genetic relationships among cultivars assessed by RAPD and SSR markers. Scientia Horticulturae 111:280-292.


Singh K, Devanshi Sharma P, Singh R, Singh B, Koundal KR, Singh NK (2007). Assessment of genetic diversity in Ziziphus mauritiana using inter-simple sequence repeat markers. Journal of Plant Biochemistry and Biotechnology 16 (1):35-40.


Song QJ, Quiley CV, Nelson RL, Carter TE, Boema HR, Strachan JL, Cregan PB (1999). A selected set of trinucleotide simple sequence repeat markers for soybean cultivar identification. Plant Varieties and Seeds 12:207-220.


Sundaram GR, Venkatesan T, Kunnummal KV (2007). Genetic diversity among cultivars, landraces and wild relatives of rice as revealed by microsatellite markers. Journal of Applied Genetics 48:337-345.


Thomas MR, Cain P, Scott NS (1994). DNA typing of grapevines: a universal methodology and database for describing cultivars and evaluating geneting relatedness. Plant Molecular Biology 25:939-949.


Uddin MS, Sun W, He X, Jaime A, Da Silva T, Cheng Q (2014). An improved method to extract DNA from mango Mangifera indica Biologia 69(2):133-138.


Ukoskit K (2007). Development of microsatellite markers in mango (Mungifera indica L.) using 5' anchore PCR. International Journal of Science and Technology 12(3):1-7.


Viruel MA, Escribano P, Barbieri M, Ferri M, Hormaza JI (2005). Fingerprinting, embryo type and geographic differentiation in mango (Mangifera indica L., Anacardiaceae) with microsatellites. Molecular Breeding 14:383-383.


Wahdan MT, Abdelsalam AZ, El-Naggar AA, Hussein MA (2011). Preliminary horticultural studies to describe and identify of two new Egyptian mango strains using DNA fingerprint. Journal of American Science 7(2):641-650.


Yamanaka N, Hasran M, Xu DH, Tsunematsu H, Idris S, Ban T (2006). Genetic relationship and diversity of four Mangifera species revealed through AFLP analysis. Genetic Resources and Crop Evolution 53:949-954.


Yeh FC, Yang RC, Boyle T (1999). POPGENE VERSION 1.31: Microsoft Window-based free Software for Population Genetic Analysis, 1999,