Gummi-Gutta by Next-Generation Sequencing and Cross-Species
Gummi-Gutta by Next-Generation Sequencing and Cross-Species
Gummi-Gutta by Next-Generation Sequencing and Cross-Species
RESEARCH ARTICLE
1
ICAR-Indian Institute of Horticultural Research, Bengaluru 560 089, India
2
Department of Forest Biology and Tree Improvement, College of Forestry Sirsi, University of Agricultural Sciences,
Dharwad 581 401, India
3
Biodiversity International, 93.4, Dharhara Marg, Fulbari, Ward no.11, Pokhara, Nepal
4
Bioversity International, University of Horticultural Sciences, GKVK, Bengaluru 560 065, India
Abstract
Garcinia gummi-gutta (L.) Roxb. (Clusiaceae) is an endemic, semidomesticated, fruit-yielding tree species distributed in the
Western Ghats of India and Sri Lanka. Various bioactive phytochemicals, such as garcinol, benzophenones and xanthones
are isolated from G. gummi-gutta and have shown antibacterial, antiviral and antioxidant activities. We sequenced the total
genomic DNA using Illumina Hiseq 2000 platform and examined 241,141,804 bp high quality data, assembled into 773,889
contigs. In these contigs, 27,313 simple-sequence repeats (SSRs) were identied, among which mononucleotide repeats were
predominant (44.98%) followed by dinucleotide and trinucleotide repeats. Primers were designed for 9964 microsatellites
among which 32 randomly selected SSR primer pairs were standardized for amplication. Polymerase chain reaction (PCR)
amplication of genomic DNA in 30 G. gummi-gutta genotypes revealed polymorphic information content (PIC) across all
32 loci ranging from 0.867 to 0.951, with a mean value of 0.917. The observed and expected heterozygosity ranged from 0.00
to 0.63 and 0.896 to 0.974, respectively. Alleles per locus ranged from 12 to 27. This is the rst report on the development of
genomic SSR markers in G. gummi-gutta using next-generation sequencing technology. The genomic SSR markers developed
in this study will be useful in identication, mapping, diversity and breeding studies.
[Ravishankar K. V., Vasudeva R., Hemanth B., Sandya B. S., Sthapit B. R., Parthasarathy V. A. and Rao V. R. 2017 Isolation and charac-
terization of microsatellite markers in Garcinia gummi-gutta by next-generation sequencing and cross-species amplication. J. Genet. 96,
xxxx]
Introduction rind of the uppage fruit has been traditionally used in India
and Sri Lanka as a culinary additive and sh preserva-
Garcinia gummi-gutta (L.) Roxb., popularly known as Mal-
tive (Samarajeewa and Shanmugapirabu 1983; Bhagyavanth
abar Gamboge (in English) uppage (in Kannada) belongs
et al. 2010). Many studies showed that hydroxycitric acid
to the family Clusiaceae. It is a moderate-sized dioecious
(HCA), a secondary compound present in the rind of uppage
tree with round canopy, drooping branches and smooth dark
fruit is effective in weight loss (Jena et al. 2002). The
barks. It is common in lower Shola forests of the West-
ern Ghats, India, up to an altitude of 1800 mean sea level. extract obtained from this fruits has exhibited the property of
G. gummi-gutta is recognized in Ayurveda, the traditional antiobesity. Recently, a number of studies have shown that
Indian system of medicine, for better digestion and is pre- G. gummi-gutta fruit extracts are rich in HCA and are effec-
scribed against abdominal disorders and heart diseases. The tive in reducing body weight (Shara et al. 2004; Saito et al.
fruit rind is ground and used as a sour avouring spice in 2005).
the preparation of curries and to garnish sh preparations. The Germplasm utilization and conservation requires precise
information on the genetics, genetic relationships and diver-
For correspondence. E-mail: kv_ravishankar@yahoo.co.in; sity. Earlier, a few studies attempted to examine diversity
kvravi@iihr.res.in. using random amplied polymorphic DNA (RAPD) and
inter-simple sequence repeat (ISSR) markers (Mohan et al. PCR and genotyping
2012; Parthasarathy et al. 2013). However, development of
We selected 50 SSR primer sets randomly and synthe-
microsatellite markers have helped in ngerprinting unique
sized with M13 tail. These M13 tailed primers were
trees, assessing degree of diversity in the populations, and
rst screened for amplication using pooled total genomic
identifying marker tightly linked to the important agronomic
DNA from ve randomly selected genotypes. We used
traits like active ingredients, and resistance to biotic and
uorescence-based M13 tailing PCR method following
abiotic stresses. Therefore, the development of markers has
Schuelke (2000) to amplify the microsatellites in a quick,
become a prerequisite for genetic studies (Bohra et al. 2011;
accurate and efcient manner. The forward primer tailed
Dutta et al. 2011). Keeping this in view, here, we report the
with 5 -GTAAAACGACGGCCAGT-3 and reverse primer
development and standardization of microsatellite markers
tailed with 5 -GTTTCTT-3 . PCR was carried out in 20 L
of G. gummi-gutta using next-generation sequencing (NGS)
reaction volume containing 2 L of 10 reaction buffer,
and their cross species transferability.
2.0 L of 1 mM dNTPs, 0.9 L (5 pmol) of forward, 0.9 L
reverse primers (5 pmol), labelled M13 probes (HEX, NED,
Materials and methods VIC, TET) 1.2 L (5 pmol), 5.0 L (5075 ng) of tem-
plate genomic DNA, 0.8 L (2 U) of Taq DNA polymerase
Plant materials
and 7.2 L of nuclease free water. The PCR cycling prole
The total genomic DNA from G. gummi-gutta was used for was: initial denaturation at 94 C for 2 min, followed by 35
genome sequencing. We have also included Garcinia indica cycles of 94 C for 30 s, 55 C for 30 s, 72 C for 1 min and a
and Garcinia morella to examine cross species transferabil- nal extension at 72 C for 5 min. PCR reaction was carried
ity of isolated microsatellite markers. The leaf material was out using thermocycler (Eppendrof Master Cycler Gradient,
obtained from the germplasm collection of the College of Germany). Amplied products were initially separated on
Forestry, Sirsi (University of Agricultural Sciences, Dhar- 3% agarose gel to conrm the amplication. Finally, 32 SSR
wad), India (see details in table 1 in electronic supplementary primers were selected based on amplication of clear PCR
material at http://www.ias.ac.in/jgenet/). The leaf samples products. These primers were employed for amplication of
were obtained from Jaddi Gadde collection (14 44 13.2 N 30 genotypes of G. gummi-gutta and one genotype each of
latitude and 74 43 03.2 E longitude with altitude of 498 m). G. indica and G. morella. The PCR products were separated
The authenticated herbarium specimens were deposited to using automated DNA Sequencer (Applied Biosystems, ABI
the herbarium of College of Forestry, Department of Forest 3730 DNA Analyzer) through capillary electrophoresis, at
Biology and Tree improvement. M/S Eurons facility, Bengaluru.
Journal of Genetics
Table 3. Genetic analysis of microsatellite markers developed for G. gummi-gutta.
Journal of Genetics
GGGGA CACCT
GG_KVRj151 AACAACGAGGGCGTCGT CATATCACCATCACCACC (TG)6 18 119215 0.308 0.918 0.893 0.0151 NA A
GGAAC AACACGA
GG_KVRj152 AAATGGTGGCAAGACACAC TGCCACCTTCCAGGCAC (AC)6 15 146243 0.091 0.943 0.916 0.0204 A NA
ATACCTT CATGT
GG_KVRi741 ACCACTGTTCCAACCAT AGAAGGTGGTGGAGGTGTA (CT)6 20 391498 0.115 0.949 0.927 0.0089 NA NA
GGGCA GGTGT
Microsatellite markers in Garcinia gummi-gutta
Cross species
City, USA) for determining the exact allele size. Allele sizes
amplication
GM
NA
NA
NA
NA
NA
NA
NA
for each SSR loci were used for genetic analysis using
A
A
Cervus 3.0 software (Kalinowski et al. 2007). We have esti-
NA
NA
NA
NA
mated the number of alleles, observed heterozygosity (Ho ),
GI
A
A
A
A
expected heterozygosity (He ) and PIC. Probability of iden-
Expected Polymorphic Probability
of allele size range heterozygosity heterozygosity information of identity
tity (PI) was analysed for each SSR loci using Identity 1.0
0.0070
0.0163
0.0622
0.0144
0.0093
0.0165
0.0245
0.0089
0.0114
(PI)
0.932
0.897
0.888
0.916
0.867
0.923
0.89
0.91
0.953
0.919
0.923
0.941
0.943
0.916
0.932
0.896
0.945
(He )
0.233
0.208
0.593
0.179
0.208
0.259
0.118
0.115
0.63
175268
103192
105178
227302
226337
309411
110211
24
17
13
18
18
12
18
14
20
(TCA)5
(TCA)5
(TTG)6
Repeat
type
TGGACCTAGCCATGCCC
TGGTGCCATCCAGTTAG
TTGCCACCTAAGCATGC
AGGGTGAGTTTTTGGC
GGGTTT
AGCCT
GGCAT
ATCCT
CTTTT
TCCCCATCACCTTCCAC
TGCTGAGATGGCAGCG
GCAACATCCCAAGTGT
ATGCCGCACACCATGC
Forward sequence 5 3
AAACGA
GGGGGT
GTGCCT
TGCTCA
ATGCTT
GCTATT
GGCCA
CACCA
GG_KVRk715
GG_KVRk769
GG_KVRk831
GG_KVRj345
GG_KVRj447
GG_KVRj541
GG_KVRj603
GG_KVRj938
Locus name
Journal of Genetics
Microsatellite markers in Garcinia gummi-gutta
Journal of Genetics
K. V. Ravishankar et al.
(-)-hydroxycitric acid extract and its effect on body weight, Untergrasser A., Cutcutache I., Koressaar T., Ye J., Faircloth B. C.,
selected organ weights, hepatic lipid peroxidation and DNA frag- Remm M. and Rozen S. G. 2012 Primer3 - new capabilities and
mentation, hematology and clinical chemistry, and histopatholog- interfaces. Nucleic Acids Res. 40, e115.
ical changes over a period of 90 days. Mol. Cell. Biochem. 260, Varshney R. K., Graner A. and Sorrells M. E. 2005 Genic
171186. microsatellites markers in plants: features and application. Trend
Sonah H., Deshmukh R. K., Sharma A., Singh V. P., Gupta D. K., Biotechnol. 23, 4855.
Gacche R. N. et al. 2011 Genome-wide distribution and organi- Wagner H. W. and Sefc K. M. 1999 IDENTITY 1.0 Centre for
zation of microsatellites in plants: an insight into marker devel- Applied Genetics. University of Agricultural Sciences, Vienna,
opment in brachypodium. PLoS One 6, e21298. Austria (http://www.uni-graz.at/~sefck/).
Thiel T., Michalek W., Varshney R. and Graner A. 2003 Exploiting Zalapa J. E., Cuevas H., Zhu H., Steffan S., Senalik D., Zeldin E.
EST databases for the development and characterization of gene- et al. 2012 Using next-generation sequencing approaches to iso-
derived SSR-markers in barley (Hordeum vulgare L). Theor. late simple sequence repeat (SSR) loci in the plant sciences. Am.
Appl. Genet. 106, 411422. J. Bot. 99, 193208.
Received 10 December 2015, in nal revised form 8 July 2016; accepted 11 July 2016
Unedited version published online: 13 July 2016
Final version published online: 17 April 2017
Journal of Genetics