Abstract

Simple sequence repeats (SSRs) are regions in DNA sequence that contain repeating motifs of length 1–6 nucleotides. These repeats are ubiquitously present and are found in both coding and non-coding regions of genome. A total of 534 complete chloroplast genome sequences (as on 18 September 2014) of Viridiplantae are available at NCBI organelle genome resource. It provides opportunity to mine these genomes for the detection of SSRs and store them in the form of a database. In an attempt to properly manage and retrieve chloroplastic SSRs, we designed ChloroSSRdb which is a relational database developed using SQL server 2008 and accessed through ASP.NET. It provides information of all the three types (perfect, imperfect and compound) of SSRs. At present, ChloroSSRdb contains 124 430 mined SSRs, with majority lying in non-coding region. Out of these, PCR primers were designed for 118 249 SSRs. Tetranucleotide repeats (47 079) were found to be the most frequent repeat type, whereas hexanucleotide repeats (6414) being the least abundant. Additionally, in each species statistical analyses were performed to calculate relative frequency, correlation coefficient and chi-square statistics of perfect and imperfect SSRs. In accordance with the growing interest in SSR studies, ChloroSSRdb will prove to be a useful resource in developing genetic markers, phylogenetic analysis, genetic mapping, etc. Moreover, it will serve as a ready reference for mined SSRs in available chloroplast genomes of green plants.

Database URL:www.compubio.in/chlorossrdb/

Introduction

Chloroplasts are semiautonomous organelles having their own genome ( 1 ) and considered to be derived from cyanobacteria through endosymbiosis ( 2 ). Apart from their well-known function of photosynthesis, i.e. the conversion of light energy to chemical energy, chloroplasts are known to play a role in the synthesis of starch, fatty acids, pigments and amino acids ( 3 ). Moreover, chloroplast genome sequences have also been widely used in plant systematics ( 4–7 ) and simple sequence repeats (SSRs) mining ( 8 , 9 ).

SSRs also known as microsatellites are the specific portions of DNA sequence that contain clusters of tandem repeating motifs of length 1–6 nucleotides ( 10 ). These repeats are supposed to be generated by slippage during replication ( 11 ) and are present in both coding as well as non-coding regions of DNA. These repeats show less polymorphism in coding sequences as compared to non-coding sequences ( 12 ). The specificity, reproducibility, co-dominance and hypervariability of SSRs make them potential molecular markers ( 13 ). The conserved flanking sequences of SSRs help in the designing of PCR primers which can be further used for the amplification of repeat sequence ( 14 ). SSRs can be used for genotyping and population level evolutionary studies ( 15 ). Moreover, these repeats play an important role in gene regulation and the importance of SSRs in the evolution of coding and non-coding regions has been proved ( 16 , 17 ). Chloroplastic SSRs (cpSSRs) also play an important role in population genetics and evolutionary studies of plants ( 18 ).

With the increase in availability of expressed sequence tags (ESTs) and complete genome sequences in biological databases, in silico mining approaches proved to be useful in the identification of SSRs ( 8 , 9 , 19 , 20 ). Consequently, a large number of SSR-specific databases including MICdb ( 21 ), Cotton Marker Database ( 22 ), EuMicroSatdb ( 23 ), PIPEMicroDB ( 24 ), ChloroMitoSSRDB ( 25 ) and MitoSatPlant ( 26 ) have been developed.

This study is an attempt to develop a comprehensive, user friendly, specialized database of cpSSRs mined from complete chloroplast genome sequences of green plants (Viridiplantae). To the best of our knowledge, among all the SSR-specific databases available, this is the only database of SSRs which provides information about perfect, imperfect and compound SSRs along with statistical analyses of the repeats identified. The database includes pre-calculated density of SSRs, average length of SSRs, repeat type frequencies, chi-square statistics, relative values, their correlation coefficient, SSR-specific PCR primers, etc.

Materials and Methods

Data mining and statistical analysis

The information included in the ChloroSSRdb was mined from 534 completely sequenced chloroplast genomes of Viridiplantae. These genomes belong to 39 algae, 9 bryophytes, 17 pteridophytes, 41 gymnosperms and 428 angiosperms. Genbank (*.gbk) and fasta (*.fna) formatted files of these genomes were retrieved from NCBI ( ftp://ftp.ncbi.nih.gov/genomes/Chloroplasts/plastids/ ). For SSRs scanning, a minimum length criterion of repeating motif ≥12 mono, ≥6 di, ≥4 tri and ≥3 for tetra, penta and hexa nucleotide repeats was used. Perfect and compound SSRs were identified using microsatellite identification tool (MISA; http://pgrc.ipk-gatersleben.de/misa/down load/misa. pl ) which fetches non-redundant perfect and compound microsatellites from a given DNA sequence. The number of intervening nucleotides between motifs of compound microsatellites was set to 0. Imperfect Microsatellite Extractor (IMEx 2.0) ( 27 ) with the imperfection percentage of 10% was employed to get imperfect SSRs. In-house developed Perl scripts were used to parse results of MISA and IMEx along with additional information from Genbank and fasta formatted files. Primer3 ( http://bioinfo.ut.ee/primer3-0.4.0/ ) ( 28 ) with default parameters was used to generate PCR primers considering 200 bases of SSR flanking regions. In addition to this, chi-square (χ 2 ) values of identified SSRs, their relative values and correlation coefficient ( r ) between the relative values of perfect and imperfect SSRs were calculated as follows:
where O i and E i are the observed and expected values of i th observation.
where X and Y are the relative frequencies of perfect and imperfect SSRs, respectively. The workflow of ChloroSSRdb is presented in Figure 1 .
The workflow of ChloroSSRdb.
Figure 1.

The workflow of ChloroSSRdb.

Database design and web interface

ChloroSSRdb is based on relational database management system and was developed using SQL server 2008. It follows client–server architecture in which the communication is one-to-one and takes place between client and server without any intermediate. The database contains a total of 19 tables. Each table uses the accession number of chloroplast genome sequence as a unique identifier (primary key). The database can be accessed through an interactive, easy to use interface developed in ASP.NET.

Results and Discussion

The front end of ChloroSSRdb provides a user-friendly browsing facility to look for the SSRs information in respective organism. Navigation to different pages are provided to link mined data with other information and every page contains a hyperlink to download the displayed information in MS Excel file.

The home page ( Figure 2 A) displays a broader classification of chloroplast genomes as algae, bryophytes, pteridophytes, gymnosperms and angiosperms where angiosperms are further classified as monocots and eudicots. The mined information of mono–hexa repeat types can be accessed for both perfect and imperfect SSRs ( Figure 2 A). The organism name is directly linked to the taxonomic page and its accession id to GenBank page at NCBI which enable user to fetch taxonomic and genomic DNA information, respectively. The repeat counts are hotlinked to the information page where frequency of mono–hexa repeats, perfect, imperfect, compound, overlapping compound [overlap of few bases of previous SSR with next SSR, e.g. (ACC) n (CT) n ] and perfect compound SSRs are displayed. Additional information such as coding (CDS, tRNA, rRNA), non-coding (intergenic and intragenic), coding–non-coding regions (occurrence of few bases in coding as well as in non-coding regions or vice versa), gene id, protein id, total density and average length of SSRs are provided ( Figure 2 B). The genomic location of SSRs provided with additional information will facilitate in the determination of their functional roles. Furthermore, links to designed PCR primers and alignment (to report substitutions or indels in imperfect SSRs) are available ( Figure 2 C and D). SSR flanking regions of 200 nucleotides are available with primer sequences. These primers can be used to develop SSR-based markers, for transferability studies across species, within a genus, across genera, and for the experimental validation of polymorphism. In addition to this, the designed primers can be used to check length polymorphism of SSRs in different species which can be helpful in species identification.

 Browsing activity of ChloroSSRdb. ( A ) Home page showing name of organisms with SSR (mono–hexa) frequency. ( B ) Information of selected organism. ( C ) Primer sequences of SSRs. ( D ) Alignment of imperfect SSR with expected perfect SSR. ( E ) Advanced search page. ( F ) Results of advanced search.
Figure 2.

Browsing activity of ChloroSSRdb. ( A ) Home page showing name of organisms with SSR (mono–hexa) frequency. ( B ) Information of selected organism. ( C ) Primer sequences of SSRs. ( D ) Alignment of imperfect SSR with expected perfect SSR. ( E ) Advanced search page. ( F ) Results of advanced search.

The advanced search option caters information on SSRs by motif specific search and filters the results on the basis of region, repeat type and the organism classification ( Figure 2 E and F). The descriptive statistics link ( Figure 2 B–D) opens a graphical representation ( Figure 3 ) which depicts the relative frequencies (mono–hexa) of perfect and imperfect SSRs, their correlation coefficient value and the chi-square statistics. The range of correlation coefficients are (−0.33, 0.98) where 50% of the species have significantly positive value. Moreover, some of the cases are having significantly negative correlation coefficient values. The calculated chi-square values in almost all cases are much greater than the tabulated value at 5% level of significance (5 df). This clearly rejects the null hypothesis, i.e. the data follow a natural specified distribution, which concludes that the observed values are significant and the differences are not only due to chance.

Chi-square statistic, relative values of perfect and imperfect SSRs along with their correlation coefficient.
Figure 3.

Chi-square statistic, relative values of perfect and imperfect SSRs along with their correlation coefficient.

For user’s ease, a tutorial is provided to easily browse ChloroSSRdb. The statistics tab displays overall statistical information of ChloroSSRdb which includes largest and smallest genomes mined, total number of SSRs mined, number of primers designed, most and least abundant repeats and organism with maximum SSR density. Additionally, a bar graph is provided that depicts the year-wise submission of chloroplast genomes of green plants in NCBI.

Conclusion

An easy to use, comprehensive database of cpSSRs mined from 534 chloroplast genomes has been developed. We hope that ChloroSSRdb will appear to be a useful resource for researchers interested in the study of cpSSRs. The statistical analyses of mined data will aid the scientific community to understand the distribution and pattern of cpSSRs evolution among plant lineages. ChloroSSRdb will be regularly updated in accordance with the sequence entries in NCBI and we will expand it by providing as much information as possible.

Acknowledgements

We are grateful to Professor Aditya Shastri, Vice Chancellor, Banasthali University, and Dr C. K. Jha, HOD, Department of Computer Science for encouragement. We thank three anonymous reviewers for constructive suggestions to improve the article.

Funding

This study was generously supported by University Grants Commission major research project grant (F.No. 42-138/2013) awarded to A.S.

Conflict of interest. None declared.

References

1

Kabeya
Y.
Miyagishima
S.Y.
(
2013
)
Chloroplast DNA replication is regulated by the redox state independently of chloroplast division in Chlamydomonas reinhardtii
.
Plant Physiol.
,
161
,
2102
2112
.

2

Wicke
S.
Schneeweiss
G.M.
dePamphilis
C.W.
et al.  . (
2011
)
The evolution of the plastid chromosome in land plants: gene content, gene order, gene function
.
Plant Mol. Biol.
,
76
,
273
297
.

3

Neuhaus
H.E.
Emes
M.J.
(
2000
)
Nonphotosynthetic metabolism in plastids
.
Annu. Rev. Plant Physiol. Plant Mol. Biol.
,
51
,
111
140
.

4

Olmstead
R.G.
Palmer
J.D.
(
1994
)
Chloroplast DNA systematics: a review of methods and data analysis
.
Am. J. Bot.
,
81
,
1205
1224
.

5

Shanker
A.
Sharma
V.
Daniell
H.
(
2011
)
Phylogenomic evidence of bryophytes’ monophyly using complete and incomplete datasets from chloroplast proteomes
.
J. Plant Biochem. Biotech.
,
20
,
288
292
.

6

Shanker
A.
(
2013a
)
Paraphyly of bryophytes inferred using chloroplast sequences
.
Arch. Bryol.
,
163
,
1
5
.

7

Shanker
A.
(
2013b
)
Combined data from chloroplast and mitochondrial genome sequences showed paraphyly of bryophytes
.
Arch. Bryol.
,
171
,
1
9
.

8

Shanker
A.
(
2013c
)
Mining of simple sequence repeats in chloroplast genome of a parasitic liverwort: Aneura mirabilis
.
Arch. Bryol.
,
196
,
1
4
.

9

Shanker
A.
(
2014
)
Computationally mined microsatellites in chloroplast genome of Pellia endiviifolia
.
Arch. Bryol.
,
199
,
1
5
.

10

Toth
G.
Gaspari
Z.
Zurka
J.
(
2000
)
Microsatellites in different eukaryotic genome, survey and analysis
.
Genome Res.
,
10
,
1967
1981
.

11

Tautz
D.
Schlotterer
C.
(
1994
)
Simple sequences
.
Curr. Opin. Genet. Dev.
,
4
,
832
837
.

12

Hancock
J.M.
(
1995
)
The contribution of slippage-like processes to genome evolution
.
J. Mol. Evol.
,
41
,
1038
1047
.

13

Squirrell
J.
Hollingsworth
P.M.
Woodhead
M.
et al.  . (
2003
)
How much effort is required to isolate nuclear microsatellites from plants? Mol
.
Ecol.
,
12
,
1339
1348
.

14

Botstein
D.
White
R.L.
Skolnick
M.
et al.  . (
1980
)
Construction of a genetic linkage map in man using restriction fragment length polymorphisms
.
Am. J. Hum. Genet.
,
32
,
314
331
.

15

Sung
W.
Tucker
A.
Bergeron
R.D.
et al.  . (
2010
)
Simple sequence repeat variation in the Daphnia pulex genome
.
BMC Genomics
,
11
,
691
.

16

Kashi
Y.
King
D.G.
(
2006
)
Simple sequence repeats as advantageous mutators in evolution
.
Trends Genet.
,
22
,
253
259
.

17

Farazi
T.A.
Juranek
S.A.
Tuschl
T.
(
2008
)
The growing catalog of small RNAs and their association with distinct Argonaute/Piwi family members
.
Development
,
135
,
1201
1214
.

18

Provan
J.
Powell
W.
Hollingsworth
P.M.
(
2001
)
Chloroplast microsatellite: new tool for studies in plant ecology and evolution
.
Trends Eco. Evol.
,
16
,
142
147
.

19

Shanker
A.
Bhargava
A.
Bajpai
R.
et al.  . (
2007
)
Bioinformatically mined simple sequence repeats in expressed sequences of Citrus sinensis
.
Sci. Hort.
,
113
,
353
361
.

20

Shanker
A.
Singh
A.
Sharma
V.
(
2007a
)
In silico mining in expressed sequences of Neurospora crassa for identification and abundance of microsatellites
.
Microbiol. Res.
,
162
,
250
256
.

21

Sreenu
V.B.
Alevoor
V.
Nagaraju
J.
et al.  . (
2003
)
MICdb: database of prokaryotic microsatellites
.
Nucleic Acids Res.
,
31
,
106
108
.

22

Blenda
A.
Scheffler
J.
Scheffler
B.
et al.  . (
2006
)
CMD: a cotton microsatellite database resource for Gossypium genomics
.
BMC Genomics
,
7
,
132
.

23

Aishwarya
V.
Grover
A.
Sharma
P.C.
(
2007
)
EuMicroSatdb: a database for microsatellites in the sequenced genomes of eukaryotes
.
BMC Genomics
,
8
,
225
.

24

Sarika
Arora
V.
Iquebal
M.A.
et al.  . (
2013
)
PIPEMicroDB: microsatellite database and primer generation tool for pigeonpea genome
.
Database
,
2013: article ID bas 054; doi:10.1093/database/bas054
.

25

Sablok
G.
Mudunuri
S.B.
Patnana
S.
et al.  . (
2013
)
ChloroMitoSSRDB: open source repository of perfect and imperfect repeats in organelle genomes for evolutionary genomics
.
DNA Res.
,
20
,
127
133
.

26

Kumar
M.
Kapil
A.
Shanker
A.
(
2014
)
MitoSatPlant: mitochondrial microsatellites database of viridiplantae
.
Mitochondrion
,

27

Mudunuri
S.B.
Nagarajaram
H.A.
(
2007
)
IMEx: imperfect microsatellite extractor
.
Bioinformatics
,
23
,
1181
1187
.

28

Untergrasser
A.
Cutcutache
I.
Koressaar
T.
et al.  . (
2012
)
Primer3: new capabilities and interfaces
.
Nucleic Acids Res.
,
40
,
e115
.

Author notes

Citation details: Kapil,A., Rai,P.K. and Shanker,A. ChloroSSRdb: a repository of perfect and imperfect chloroplastic simple sequence repeats (cpSSRs) of green plants. Database (2014) Vol. 2014: article ID bau107; doi:10.1093/database/bau107

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.