Identification of SSRs and differentially expressed genes in two cultivars of celery (Apium graveolens L.) by deep transcriptome sequencing

Abstract

Celery (Apium graveolens L.) is one of the most important and widely grown vegetables in the Apiaceae family. Due to the lack of comprehensive genomic resources, research on celery has mainly utilized physiological and biochemical approaches, rather than molecular biology, to study this crop. Transcriptome sequencing has become an efficient and economic technology for obtaining information on gene expression that can greatly facilitate molecular and genomic studies of species for which a sequenced genome is not available. In the present study, 15 893 516 and 19 818 161 high-quality sequences were obtained by RNA-seq from two celery varieties ‘Ventura’ and ‘Jinnan Shiqin’, respectively. The obtained reads were assembled into 39 584 and 41 740 unigenes with mean lengths of 683 bp and 690 bp, respectively. A total of 1939 simple sequence repeat (SSR) markers were identified in ‘Ventura’ and 2004 SSRs in ‘Jinnan Shiqin’. Di-nucleotide repeats were the most common repeat motif, accounting for 55.49% and 54.84% in ‘Ventura’ and ‘Jinnan Shiqin’, respectively. A comparison of expressed genes between the two libraries, identified 338 differentially expressed genes (DEGs). Three hundred and three of the DEGs were annotated based on a sequence similarity search utilizing eight public databases. Additionally, the expression profile of eight annotated DEGs was characterized in response to abiotic stresses. The collective data generated in the present research represent a valuable resource for further genetic and molecular studies in celery.

References

[1] Krishnamurthy K . Celery. In: Chemistry of Spices. Cambridge, MA: CABI; 2008. p401.

[2] Craig WJ . Health-promoting properties of common herbs. Am J Clin Nutr 1999; 70(3 Suppl): 491S-499S.

[3] Raffo A, Sinesio F, Moneta E, Nardo N, Peparaio M, Paoletti F . Internal quality of fresh and cold stored celery petioles described by sensory profile, chemical and instrumental measurements. Eur Food Res Technol 2005; 222: 590-599.

[4] Hostetler GL, Riedl KM, Schwartz SJ . Effects of food formulation and thermal processing on flavones in celery and chamomile. Food Chem 2013; 141: 1406-1411.

[5] Toth G, Gaspari Z, Jurka J . Microsatellites in different eukaryotic genomes: survey and analysis. Genome Res 2000; 10: 967-981.

[6] Rongwen J, Akkaya MS, Bhagwat AA, Lavi U, Cregan PB . The use of microsatellite DNA markers for soybean genotype identification. Theor Appl Genet 1995; 90: 43-48.

[7] Gharghani A, Zamani Z, Talaie A et al. Genetic identity and relationships of Iranian apple (Malus × domestica Borkh.) cultivars and landraces, wild Malus species and representative old apple cultivars based on simple sequence repeat (SSR) marker analysis . Genet Resour Crop Evol 2009; 56: 829-842.

[8] Chen XM, Luo YH, Xia XC et al. Chromosomal location of powdery mildew resistance gene Pm16 in wheat using SSR marker analysis . Plant Breeding 2005; 124: 225-228.

[9] Marquez Y, Brown JW, Simpson C, Barta A, Kalyna M . Transcriptome survey reveals increased complexity of the alternative splicing landscape in Arabidopsis. Genome Res 2012; 22: 1184-1195.

[10] Garg R, Patel RK, Jhanwar S et al. Gene discovery and tissue-specific transcriptome analysis in chickpea with massively parallel pyrosequencing and web resource development. Plant Physiol 2011; 156: 1661-1678.

[11] Yuan Y, Song L, Li M et al. Genetic variation and metabolic pathway intricacy govern the active compound content and quality of the Chinese medicinal plant Lonicera japonica thunb . BMC Genomics 2012; 13: 195.

[12] Yang SS, Tu ZJ, Cheung F et al. Using RNA-Seq for gene identification, polymorphism detection and transcript profiling in two alfalfa genotypes with divergent cell wall composition in stems. BMC Genomics 2011; 12: 199.

[13] Zhang R, Marshall D, Bryan GJ, Hornyik C . Identification and characterization of miRNA transcriptome in potato by high-throughput sequencing. PloS ONE 2013; 8: e57233.

[14] Venturini L, Ferrarini A, Zenoni S et al. De novo transcriptome characterization of Vitis vinifera cv. Corvina unveils varietal diversity . BMC Genomics 2013; 14: 41.

[15] Chen X, Zhu W, Azam S et al. Deep sequencing analysis of the transcriptomes of peanut aerial and subterranean young pods identifies candidate genes related to early embryo abortion. Plant Biotechnol J 2013; 11: 115-127.

[16] Fu N, Wang Q, Shen HL . De novo assembly, gene annotation and marker development using illumina paired-end transcriptome sequences in celery (Apium graveolens L.) . PloS ONE 2013; 8: e57686.

[17] Li H. A celery variety-‘Jinnan Shiqin’. Sci Technol Tianjin Agric For 2004; (04): 21.

[18] Grabherr MG, Haas BJ, Yassour M et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol 2011; 29: 644-652.

[19] Thiel T, Michalek W, Varshney RK, Graner A . Exploiting EST databases for the development and characterization of gene-derived SSR-markers in barley (Hordeum vulgare L.) . TAG Theor Appl Genet 2003; 106: 411-422.

[20] Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B . Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods 2008; 5: 621-628.

[21] Audic S, Claverie JM . The significance of digital gene expression profiles. Genome Res 1997; 7: 986-995.

[22] Romualdi C, Bortoluzzi S, d’Alessi F, Danieli GA . IDEG6: a web tool for detection of differentially expressed genes in multiple tag sampling experiments. Physiol Genomics 2003; 12: 159-162.

[23] Ashburner M, Ball CA, Blake JA et al. Gene Ontology: tool for the unification of biology. Nat Genet 2000; 25: 25-29.

[24] Tatusov RL, Fedorova ND, Jackson JD et al. The COG database: an updated version includes eukaryotes. BMC Bioinformatics 2003; 4: 41.

[25] Kanehisa M, Goto S . KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res 2000; 28: 27-30.

[26] Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M . Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics (Oxford, England) 2005; 21: 3674-3676.

[27] Ye J, Fang L, Zheng H et al. WEGO: a web tool for plotting GO annotations. Nucleic Acids Res 2006; 34: W293-W297.

[28] Pfaffl MW . A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res 2001; 29: e45.

[29] Weber JL . Informativeness of human (dC-dA)n.(dG-dT)n polymorphisms. Genomics 1990; 7: 524-530.

[30] Wang Z, Gerstein M, Snyder M . RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet 2009; 10: 57-63.

[31] Nicolai M, Pisani C, Bouchet JP, Vuylsteke M, Palloix A . Discovery of a large set of SNP and SSR genetic markers by high-throughput sequencing of pepper (Capsicum annuum) . Genet Mol Res 2012; 11: 2295-2300.

[32] Strozycki PM, Szymanski M, Szczurek A, Barciszewski J, Figlerowicz M . A new family of ferritin genes from Lupinus luteus-comparative analysis of plant ferritins, their gene structure, and evolution . Mol Biol Evol 2010; 27: 91-101.

[33] Bachlava E, Taylor CA, Tang S et al. SNP discovery and development of a high-density genotyping array for sunflower. PloS ONE 2012; 7: e29814.

[34] Wei W, Qi X, Wang L et al. Characterization of the sesame (Sesamum indicum L.) global transcriptome using Illumina paired-end sequencing and development of EST-SSR markers . BMC Genomics 2011; 12: 1-13.

[35] Jung S, Abbott A, Jesudurai C, Tomkins J, Main D . Frequency, type, distribution and annotation of simple sequence repeats in Rosaceae ESTs. Funct Integr Genomics 2005; 5: 136-143.

[36] Fraser LG, Harvey CF, Crowhurst RN, Silva HN . EST-derived microsatellites from Actinidia species and their potential for mapping . Theor Appl Genet 2004; 108: 1010-1016.

[37] Gupta PK, Rustgi S, Sharma S, Singh R, Kumar N, Balyan HS . Transferable EST-SSR markers for the study of polymorphism and genetic diversity in bread wheat. Mol Genet Genomics 2003; 270: 315-323.

[38] Scott KD, Eggler P, Seaton G et al. Analysis of SSRs derived from grape ESTs. Theor Appl Genet 2000; 100: 723-726.

[39] Cordeiro GM, Casu R, McIntyre CL, Manners JM, Henry RJ . Microsatellite markers from sugarcane (Saccharum spp.) ESTs cross transferable to erianthus and sorghum . Plant Sci 2001; 160: 1115-1123.

[40] Harr B, Schlotterer C . Long microsatellite alleles in Drosophila melanogaster have a downward mutation bias and short persistence times, which cause their genome-wide underrepresentation . Genetics 2000; 155: 1213-1220.

[41] Schorderet DF, Gartler SM . Analysis of CpG suppression in methylated and nonmethylated species. Proc Natl Acad Sci USA 1992; 89: 957-961.

[42] Depeiges A, Goubely C, Lenoir A et al. Identification of the most represented repeated motifs in Arabidopsis thaliana microsatellite loci . Theor Appl Genet 1995; 91: 160-168.

[43] Xin D, Sun J, Wang J et al. Identification and characterization of SSRs from soybean (Glycine max) ESTs . Mol Biol Rep 2012; 39: 9047-9057.

[44] Morgante M, Hanafey M, Powell W . Microsatellites are preferentially associated with nonrepetitive DNA in plant genomes. Nat Genet 2002; 30: 194-200.

[45] Zhou M, Tong CF, Shi JS . A preliminary analysis of synonymous codon usage in poplar species. Zhi Wu Sheng Li Yu Fen Zi Sheng Wu Xue Xue Bao 2007; 33: 285-293.

[46] Morton BR, Wright SI . Selective constraints on codon usage of nuclear genes from Arabidopsis thaliana . Mol Biol Evol 2007; 24: 122-129.

[47] Tao X, Gu YH, Wang HY et al. Digital gene expression analysis based on integrated de novo transcriptome assembly of sweet potato [ Ipomoea batatas (L.) Lam] . PloS ONE 2012; 7: e36234.

[48] Ashraf M. Inducing drought tolerance in plants: recent advances. Biotechnol Adv 2010; 28: 169-183.

[49] Ashraf M, Foolad MR . Crop breeding for salt tolerance in the era of molecular markers and marker-assisted selection. Plant Breeding 2013; 132: 10-20.

[50] Yue B, Xue W, Xiong L et al. Genetic basis of drought resistance at reproductive stage in rice: separation of drought tolerance from drought avoidance. Genetics 2006; 172: 1213-1228.

[51] Yang Q, Chen ZZ, Zhou XF et al. Overexpression of SOS (Salt Overly Sensitive) genes increases salt tolerance in transgenic Arabidopsis. Mol Plant 2009; 2: 22-31.