Abstract

Celery (Apium graveolens L.) is a plant belonging to the Apiaceae family, and a popular vegetable worldwide because of its abundant nutrients and various medical functions. Although extensive genetic and molecular biological studies have been conducted on celery, its genomic data remain unclear. Given the significance of celery and the growing demand for its genomic data, the whole genome of ‘Q2-JN11’ celery (a highly inbred line obtained by artificial selfing of ‘Jinnan Shiqin’) was sequenced using HiSeq 2000 sequencing technology. For the convenience of researchers to study celery, an online database of the whole-genome sequences of celery, CeleryDB, was constructed. The sequences of the whole genome, nucleotide sequences of the predicted genes and amino acid sequences of the predicted proteins are available online on CeleryDB. Home, BLAST, Genome Browser, Transcription Factor and Download interfaces composed of the organizational structure of CeleryDB. Users can search the celery genomic data by using two user-friendly query tools: basic local alignment search tool and Genome Browser. In the future, CeleryDB will be constantly updated to satisfy the needs of celery researchers worldwide.

Database URL: http://apiaceae.njau.edu.cn/celerydb

Introduction

Celery (Apium graveolens L.) is a plant belonging to the Apiaceae family originated from the Middle East and the Mediterranean, and is one of the most important vegetables worldwide (1). Celery is widely cultivated owing to its low calorie count and abundant celluloses, vitamins and carotenes. Previous studies have found that celery possesses numerous medicinal functions, such as inhibiting cancer cell growth and decreasing blood pressure (2, 3). Celery cultivated in China is mainly classified into two groups, namely, Chinese celery (also known as local celery) and Western celery (introduced from Western countries). The ‘Q2-JN11’, a local celery, is a highly inbred line obtained by artificial selfing of ‘Jinnan Shiqin’. Physiological and molecular investigations are necessary to address the increasing demand for celery.

With the development of sequencing technology and molecular biology, many genetic researches on celery were reported. The transcriptome sequences of A. graveolens cv. ‘Ventura’ leaves at different stages were de novo assembled, and the results provided useful information on lignin accumulation in celery (4). Fu et al. (5) recognized several molecular markers of celery by transcriptome sequencing. Simple sequence repeat markers and differentially expressed genes were identified from two celery cultivars using deep transcriptome sequencing (6). Comparative proteomic analysis was conducted on celery to understand the defense system under temperature stresses (7). High-throughput sequencing of small RNAs of celery varieties identified the abiotic stress-related microRNAs (8, 9). The application of next-generation sequencing (NGS) technology to plants has provided considerable information on the genetic resources for researchers worldwide (10). However, to our knowledge, a public genomic database of celery is currently unavailable. On the basis of the whole-genome sequences of A. graveolens cv. ‘Q2-JN11’, we constructed CeleryDB (http://apiaceae.njau.edu.cn/celerydb), a genomic database for celery.

The de novo assembled whole-genome sequences of A. graveolens cv. ‘Q2-JN11’ are available on CeleryDB. The CeleryDB website consists of five interfaces, namely, Home, BLAST, Genome Browser, transcription factor (TF) and Download. The collective data presented in this database is expected to provide valuable resources for genetic and genomic studies on celery.

Data resources

Plant materials

The data presented in the current version of CeleryDB was derived from the genome sequence of Q2-JN11 celery (a highly inbred line obtained by artificial selfing of ‘Jinnan Shiqin’). Celery seeds were deposited and sowed at the State Key Laboratory of Crop Genetics and Germplasm Enhancement, Nanjing Agricultural University. The celery seedling was grown under the condition of 12-h light at 22°C and 12-h dark at 18°C, with a relative humidity of 60–70%. The whole genome sequences, nucleotide sequences of predicted genes and amino acid sequences of predicted proteins are available on CeleryDB.

Genome sequencing and assembly

The genomic DNA was extracted from the young leaves of ‘Q2-JN11’ celery using CTAB method with some modifications (11). The genomic DNA was sequenced using the HiSeq 2000 platform (BGI-Shengzhen). Raw data were filtered for adaptor contamination and low quality by using CutAdapt prior to assembly. Then, the cleaned data were assembled using SOAPdenovo2 software (http://soap.genomics.org.cn/soapdenovo.html) (12). Finally, we obtained an assembly of 3.18 Gb.

Gene prediction and annotation

A total of 34 277 putative genes were predicted using Augustus 3.2.2 software and SNAP program (13, 14). The average length of putative genes and numbers of exons per putative genes was 3267 bp and 5.27, respectively. The predicted genes were annotated based on the alignments to the NCBI non-redundant protein sequence, Swiss-Prot, TrEMBL and gene ontology (GO) databases with BLASTp at E values of 1 × 10−4 (15, 16). The GO IDs for the predicted genes were obtained using Blast2GO (17).

Identification and classification of TF

TFs bind to special sites of the target gene to regulate the gene transcription during plant biological processes (18, 19). TFs can be grouped into various families on the basis of the DNA-binding domains of the protein sequences (20). The sequences of other known TFs were downloaded from PlnTFDB (21). The conserved domains of various TF families were used as queries to search against the predicted celery proteins to obtain the celery TFs. Prediction of TFs in the CeleryDB database was accomplished using HMMER software (22). In the celery genome, a total of 1698 TFs that were classified under 40 families were identified and provided (Figure 1). Among the 40 families, the FAR1 family showed the largest number of TF members, followed by the ERF and MYB TF families. The numbers of TFs in the FAR1, ERF and MYB families were 185, 161 and 154, respectively.

Distribution of various TF families in celery.
Figure 1.

Distribution of various TF families in celery.

Database construction

Implementation of CeleryDB

To make the genomic data of celery available, the user-friendly website database CeleryDB was developed and constructed. We used the Linux (CentOS6.2) system as the server, Apache HTTP as the web server and PHP5 for Web development. Perl scripts, HTML and JavaScript were also used to build the website. The basic local alignment search tool (BLAST) and Genome Browser services were installed in this website database (23). CeleryDB allows users to access BLAST and browse the celery genome data, and to download the sequences of putative genes and proteins. The CeleryDB website consists of five interfaces, namely, Home, BLAST, Genome Browser, TF and Download. The Home interface provided an introduction and images of celery (Figure 2). Here, users can acquire the basic information on celery.

Homepage of CeleryDB.
Figure 2.

Homepage of CeleryDB.

Basic local alignment search tool

The BLAST program was embedded in the CeleryDB Web interface to allow users to perform sequence alignment (24). Users can obtain target genes from the celery database on the basis of sequence similarity by using BLAST program. The sequences of the whole genome, putative genes, and putative proteins of celery are available through the BLAST program. Prior to BLAST, users should enter the query sequence in FASTA format, select the algorithm (BLASTp or BLASTx), and set the associated parameters (Expect threshold and Matrix) (Figure 3). Then, clicking the Submit icon allows users to go to the query result interface. The IDs corresponding to the query sequence of celery are listed and ordered based on the alignment scores. Clicking the Bits icon behind the gene IDs provides the link to the alignment results interface.

Detailed procedures for BLAST search on CeleryDB.
Figure 3.

Detailed procedures for BLAST search on CeleryDB.

Genome browser

For exhibiting the annotation of the celery genome, a Genome Browser that well integrates the database and interactive web pages was embedded in CeleryDB (23). The Genome Browser allowed users to track the annotations of genes, mRNA, coding sequence (CDS) and transcripts of each scaffold. The various annotations of the scaffolds were marked with different icons in this browser. Users can acquire the detailed features of various annotations by clicking the corresponding icons (Figure 4).

Interface of Genome Browser on CeleryDB.
Figure 4.

Interface of Genome Browser on CeleryDB.

Transcription factor

TFs are vital regulators with highly conserved DNA-binding domains during plant growth development and stress response (18, 25). Previous studies demonstrated that numerous known TFs play significant roles in various biological processes. For example, the TFs of Dof and TCP families are involved in plant growth and development; some TFs of the NAC family are related to stress response; and numerous TFs of the MYB family function in the biosynthesis of secondary flavonoid metabolites (26–29). For convenience, the 1698 celery TFs belonging to 40 families are listed in a table on the TF interface (Figure 5). The numbers of different TF families were indicated in parenthesis. Users can download the TF data in the Download section.

Number of different TF families and the interface of TF.
Figure 5.

Number of different TF families and the interface of TF.

Download

The Download interface allows users to freely download the celery genomic data for further analysis. The available data included the scaffolds of the genome assembly, the CDS of predicted genes, the amino acid sequences of predicted proteins and the gene annotation (GO, InterPro, The best BLAST hit of NR and TF) of the celery genome (Figure 6).

Interface of Download on CeleryDB.
Figure 6.

Interface of Download on CeleryDB.

Allergenic protein genes in celery

Celery is recognized as a healthy vegetable because of its abundant nutrients over the world, whereas celery is also one of the common plant food sources that cause allergic reactions in central European human (30–32). The allergenic reaction to celery induces many symptoms such as human oral allergy, severe cases exhibited life-threatening anaphylactic reactions (33, 34). So far, several allergens were identified from celery, including Api g 1 (35), Api g 2 (36), Api g 4 (37), Api g 5 (38) and Api g 6 (39). Numerous studies indicated that Api g 1 is the most important allergen to human in celery (40, 41). To examine the quality of CeleryDB database, we used the sequence of known Api g 1 (GenBank accession No. Z75662.2) as query to search against the celery genome to obtain the allergenic protein gene. Search results indicated that Agr42308 in CeleryDB was the most similar to the sequence of Api g 1, with similarity of 100%. Agr42308 sequence contained a 480 bp open reading frame that encoded 159 amino acids. The above retrieved sequence in CeleryDB will be useful in future studies of celery allergens. These results indicated that the data in CeleryDB is efficient and accurate for celery genome research.

Discussion and future plans

Celery is a popular vegetable worldwide because of its abundant nutrients and various medicinal effects (1). Despite the numerous genetic and molecular biology studies on celery, no public database of the celery genome is currently available worldwide. In view of the significance of celery and the development of bioinformatics, an online database based on the genome of ‘Q2-JN11’ celery named CeleryDB was constructed. To our knowledge, CeleryDB is the first public genome database with functional annotations for celery. The sequences of the whole genome, putative genes and putative proteins of celery are available on CeleryDB. In addition, CeleryDB identifies and provides the putative TFs from the whole genome sequence of celery. To enable users to obtain the celery genomic data, we embedded two user-friendly query tools, namely, BLAST and Genome Browser, into CeleryDB.

CeleryDB was targeted to be a flexible computational platform for future genetic studies on celery. With the development of sequencing technology and the increasing studies on celery, various celery transcriptomes will be published in the coming years. Therefore, CeleryDB will be constantly updated with new information to satisfy the increasing research demands. We hope our efforts will make CeleryDB a helpful database for celery genome research.

Funding

The research was supported by the New Century Excellent Talents in University (NCET-11-0670); National Natural Science Foundation of China (31272175); Jiangsu Natural Science Foundation (BK20130027); Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).

Conflict of interest. None declared.

References

1

Li
M.Y.
,
Hou
X.L.
,
Wang
F.
et al.  (
2018
)
Advances in the research of celery, an important Apiaceae vegetable crop
.
Crit. Rev. Biotechnol
.,
38
,
172
183
.

2

Dianat
M.
,
Veisi
A.
,
Ahangarpour
A.
,
Fathi Moghaddam
H.
(
2015
)
The effect of hydro-alcoholic celery (Apium graveolens) leaf extract on cardiovascular parameters and lipid profile in animal model of hypertension induced by fructose
.
Avicenna J. Phytomed
.,
5
,
203
209
.

3

Shukla
S.
,
Gupta
S.
(
2007
)
Apigenin-mediated modulations of PI3K-Akt and MAPK signaling pathways causes growth Inhibition and cell cycle arrest in human prostate cancer cells
.
Cancer Res
.,
67
,
3350
.

4

Jia
X.L.
,
Wang
G.L.
,
Xiong
F.
et al.  (
2015
)
De novo assembly, transcriptome characterization, lignin accumulation, and anatomic characteristics: novel insights into lignin biosynthesis during celery leaf development
.
Sci. Rep
.,
5
,
8259
.

5

Fu
N.
,
Wang
Q.
,
Shen
H.L.
(
2013
)
De novo assembly, gene annotation and marker development using Illumina paired-end transcriptome sequences in celery (Apium graveolens L.)
.
PLoS One
,
8
,
e57686.

6

Li
M.Y.
,
Wang
F.
,
Jiang
Q.
et al.  (
2014
)
Identification of SSRs and differentially expressed genes in two cultivars of celery (Apium graveolens L.) by deep transcriptome sequencing
.
Hortic. Res
.,
1
,
10.

7

Huang
W.
,
Ma
H.Y.
,
Huang
Y.
et al.  (
2017
)
Comparative proteomic analysis provides novel insights into chlorophyll biosynthesis in celery under temperature stress
.
Physiol. Plant
,
161
,
468
485
.

8

Li
M.Y.
,
Wang
F.
,
Xu
Z.S.
et al.  (
2014
)
High throughput sequencing of two celery varieties small RNAs identifies microRNAs involved in temperature stress response
.
BMC Genomics
,
15
,
242
.

9

Jiang
Q.
,
Wang
F.
,
Li
M.Y.
et al.  (
2014
)
High-throughput analysis of small RNAs and characterization of novel microRNAs affected by abiotic stress in a local celery cultivar
.
Sci. Hortic
.,
169
,
36
43
.

10

McCouch
S.R.
,
McNally
K.L.
,
Wang
W.
,
Sackville Hamilton
R.
(
2012
)
Genomics of gene banks: a case study in rice
.
Am. J. Bot
.,
99
,
407
423
.

11

Rogers
S.O.
,
Bendich
A.J.
(
1985
)
Extraction of DNA from milligram amounts of fresh, herbarium and mummified plant tissues
.
Plant Mol. Biol
.,
5
,
69
76
.

12

Luo
R.
,
Liu
B.
,
Xie
Y.
et al.  (
2012
)
SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler
.
Gigascience
,
1
,
18.

13

Stanke
M.
,
Keller
O.
,
Gunduz
I.
et al.  (
2006
)
AUGUSTUS: ab initio prediction of alternative transcripts
.
Nucleic Acids Res
.,
34
,
W435
W439
.

14

Korf
I.
(
2004
)
Gene finding in novel genomes
.
BMC Bioinformatics
,
5
,
59
.

15

Ashburner
M.
,
Ball
C.A.
,
Blake
J.A.
et al.  (
2000
)
Gene ontology: tool for the unification of biology. The Gene Ontology Consortium
.
Nat. Genet
.,
25
,
25
29
.

16

Bairoch
A.
,
Apweiler
R.
(
2000
)
The SWISS-PROT protein sequence database and its supplement TrEMBL in 2000
.
Nucleic Acids Res
.,
28
,
45
48
.

17

Conesa
A.
,
Gotz
S.
,
Garcia-Gomez
J.M.
et al.  (
2005
)
Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research
.
Bioinformatics
,
21
,
3674
3676
.

18

Chen
K.
,
Rajewsky
N.
(
2007
)
The evolution of gene regulation by transcription factors and microRNAs
.
Nat. Rev. Genet
.,
8
,
93
103
.

19

Zhang
J.Z.
(
2003
)
Overexpression analysis of plant transcription factors
.
Curr. Opin. Plant Biol
.,
6
,
430
440
.

20

Riechmann
J.L.
,
Heard
J.
,
Martin
G.
et al.  (
2000
)
Arabidopsis transcription factors: genome-wide comparative analysis among eukaryotes
.
Science
,
290
,
2105
2110
.

21

Riano-Pachon
D.M.
,
Ruzicic
S.
,
Dreyer
I.
,
Mueller-Roeber
B.
(
2007
)
PlnTFDB: an integrative plant transcription factor database
.
BMC Bioinformatics
,
8
,
42.

22

Eddy
S.R.
(
2011
)
Accelerated Profile HMM Searches
.
PLoS Comput Biol
.,
7
,
e1002195.

23

Stein
L.D.
,
Mungall
C.
,
Shu
S.
et al.  (
2002
)
The generic genome browser: a building block for a model organism system database
.
Genome Res
.,
12
,
1599
1610
.

24

Johnson
M.
,
Zaretskaya
I.
,
Raytselis
Y.
et al.  (
2008
)
NCBI BLAST: a better web interface
.
Nucleic Acids Res
.,
36
,
W5
W9
.

25

Singh
K.
,
Foley
R.C.
,
Onate-Sanchez
L.
(
2002
)
Transcription factors in plant defense and stress responses
.
Curr. Opin. Plant Biol
,
5
,
430
436
.

26

Feng
K.
,
Xu
Z.S.
,
Que
F.
et al.  (
2018
)
An R2R3-MYB transcription factor, OjMYB1, functions in anthocyanin biosynthesis in Oenanthe javanica
.
Planta
,
247
,
301
315
.

27

Noguero
M.
,
Atif
R.M.
,
Ochatt
S.
,
Thompson
R.D.
(
2013
)
The role of the DNA-binding One Zinc Finger (DOF) transcription factor family in plants
.
Plant Sci
.,
209
,
32
45
.

28

Kosugi
S.
,
Ohashi
Y.
(
2002
)
DNA binding and dimerization specificity and potential targets for the TCP protein family
.
Plant J
,
30
,
337
348
.

29

Nuruzzaman
M.
,
Manimekalai
R.
,
Sharoni
A.M.
et al.  (
2010
)
Genome-wide analysis of NAC transcription factor family in rice
.
Gene
,
465
,
30
44
.

30

Han
D.
,
Row
K.H.
(
2011
)
Determination of luteolin and apigenin in celery using ultrasonic-assisted extraction based on aqueous solution of ionic liquid coupled with HPLC quantification
.
J. Sci. Food Agric
.,
91
,
2888
2892
.

31

Andre
F.
,
Andre
C.
,
Colin
L.
et al.  (
1994
)
Role of new allergens and of allergens consumption in the increased incidence of food sensitizations in France
.
Toxicology
,
93
,
77
83
.

32

Wuthrich
B.
(
2005
)
Frequency of food allergies over time - longitudinal statistics from 1978-1988
.
Allergologie
,
28
,
355
358
.

33

Pauli
G.
,
Bessot
J.C.
,
Braun
P.A.
et al.  (
1988
)
Celery allergy: clinical and biological study of 20 cases
.
Ann. Allergy
,
60
,
243
246
.

34

Ballmer-Weber
B.K.
,
Vieths
S.
,
Luttkopf
D.
et al.  (
2000
)
Celery allergy confirmed by double-blind, placebo-controlled food challenge: a clinical study in 32 subjects with a history of adverse reactions to celery root
.
J. Allergy Clin. Immunol
.,
106
,
373
378
.

35

Hoffmann-Sommergruber
K.
,
Ferris
R.
,
Pec
M.
et al.  (
2000
)
Characterization of api g 1.0201, a new member of the Api g 1 family of celery allergens
.
Int. Arch. Allergy Immunol
.,
122
,
115
123
.

36

Gadermaier
G.
,
Egger
M.
,
Girbl
T.
et al.  (
2011
)
Molecular characterization of Api g 2, a novel allergenic member of the lipid-transfer protein 1 family from celery stalks
.
Mol. Nutr. Food Res
.,
55
,
568
577
.

37

Scheurer
S.
,
Wangorsch
A.
,
Haustein
D.
,
Vieths
S.
(
2000
)
Cloning of the minor allergen Api g 4 profilin from celery (Apium graveolens) and its cross-reactivity with birch pollen profilin Bet v 2
.
Clin. Exp. Allergy
,
30
,
962
971
.

38

Bublin
M.
,
Radauer
C.
,
Wilson
I.B.
et al.  (
2003
)
Cross-reactive N-glycans of Api g 5, a high molecular weight glycoprotein allergen from celery, are required for immunoglobulin E binding and activation of effector cells from allergic patients
.
Faseb J
.,
17
,
1697
1699
.

39

Vejvar
E.
,
Himly
M.
,
Briza
P.
et al.  (
2013
)
Allergenic relevance of nonspecific lipid transfer proteins 2: identification and characterization of Api g 6 from celery tuber as representative of a novel IgE-binding protein family
.
Mol. Nutr. Food Res
.,
57
,
2061
2070
.

40

Breiteneder
H.
,
Hoffmann-Sommergruber
K.
,
O'Riordain
G.
et al.  (
1995
)
Molecular characterization of Api g 1, the major allergen of celery (Apium graveolens), and its immunological and structural relationships to a group of 17-kDa tree pollen allergens
.
Eur. J. Biochem
.,
233
,
484
489
.

41

Hoffmann-Sommergruber
K.
,
Demoly
P.
,
Crameri
R.
et al.  (
1999
)
IgE reactivity to Api g 1, a major celery allergen, in a Central European population is based on primary sensitization by Bet v 1
.
J. Allergy Clin. Immunol
.,
104
,
478
484
.

Author notes

Citation details: Feng,K., Hou,X.-L., Li,M.-Y. et al. CeleryDB: a genomic database for celery. Database (2018) Vol. 2018: article ID bay070; doi:10.1093/database/bay070

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.