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

A comprehensive review on genomic resources in medicinally and industrially important major spices for future breeding programs: Status, utility and challenges

In the global market, spices possess a high-value but low-volume commodities of commerce. The food industry depends largely on spices for taste, flavor, and therapeutic properties in replacement of cheap synthetic ones. The estimated growth rate for spices demand in the world is ∼3.19%. Since spices grow in limited geographical regions, India is one of the leading producer of spices, contributing 25-30 percent of total world trade. Hitherto, there has been no comprehensive review of the genomic resources of industrially important major medicinal spices to overcome major impediments in varietal improvement and management. This review focuses on currently available genomic resources of 24 commercially significant spices, namely, Ajwain, Allspice, Asafoetida, Black pepper, Cardamom large, Cardamom small, Celery, Chillies, Cinnamon, Clove, Coriander, Cumin, Curry leaf, Dill seed, Fennel, Fenugreek, Garlic, Ginger, Mint, Nutmeg, Saffron, Tamarind, Turmeric and Vanilla. The advent of low-cost sequencing machines has contributed immensely to the voluminous data generation of these spices, cracking the complex genomic architecture, marker discovery, and understanding comparative and functional genomics. This review of spice genomics resources concludes the perspective and way forward to provide footprints by uncovering genome assemblies, sequencing and re-sequencing projects, transcriptome-based studies, non-coding RNA-mediated regulation, organelles-based resources, developed molecular markers, web resources, databases and AI-directed resources in candidate spices for enhanced breeding potential in them. Further, their integration with molecular breeding could be of immense use in formulating a strategy to protect and expand the production of the spices due to increased global demand.

Comprehensive analysis of HSF genes from celery (Apium graveolens L.) and functional characterization of AgHSFa6-1 in response to heat stress

High temperature stress is regarded as one of the significant abiotic stresses affecting the composition and distribution of natural habitats and the productivity of agriculturally significant plants worldwide. The HSF family is one of the most important transcription factors (TFs) families in plants and capable of responding rapidly to heat and other abiotic stresses. In this study, 29 AgHSFs were identified in celery and classified into three classes (A, B, and C) and 14 subgroups. The gene structures of AgHSFs in same subgroups were conserved, whereas in different classes were varied. AgHSF proteins were predicted to be involved in multiple biological processes by interacting with other proteins. Expression analysis revealed that AgHSF genes play a significant role in response to heat stress. Subsequently, AgHSFa6-1, which was significantly induced by high temperature, was selected for functional validation. AgHSFa6-1 was identified as a nuclear protein, and can upregulate the expression of certain downstream genes (HSP98.7, HSP70-1, BOB1, CPN60B, ADH2, APX1, GOLS1) in response to high-temperature treatment. Overexpression of AgHSFa6-1 in yeast and Arabidopsis displayed higher thermotolerance, both morphologically and physiologically. In response to heat stress, the transgenic plants produced considerably more proline, solute protein, antioxidant enzymes, and less MDA than wild-type (WT) plants. Overall, this study revealed that AgHSF family members perform a key role in response to high temperature, and AgHSFa6-1 acts as a positive regulator by augmenting the ROS-scavenging system to maintain membrane integrity, reducing stomatal apertures to control water loss, and upregulating the expression level of heat-stress sensitive genes to improve celery thermotolerance.

AgMYB5, an MYB transcription factor from celery, enhanced β-carotene synthesis and promoted drought tolerance in transgenic Arabidopsis

BACKGROUND

Water shortage caused by global warming seriously affects the yield and quality of vegetable crops. β-carotene, the lipid-soluble natural product with important pharmacological value, is abundant in celery. Transcription factor MYB family extensively disperses in plants and plays regulatory roles in carotenoid metabolism and water scarcity response.

RESULTS

Here, the AgMYB5 gene encoding 196 amino acids was amplified from celery cv. 'Jinnanshiqin'. In celery, the expression of AgMYB5 exhibited transactivation activity, tissue specificity, and drought-condition responsiveness. Further analysis proved that ectopic expression of AgMYB5 increased β-carotene content and promoted drought tolerance in transgenic Arabidopsis thaliana. Moreover, AgMYB5 expression promoted β-carotene biosynthesis by triggering the expression of AtCRTISO and AtLCYB, which in turn increased antioxidant enzyme activities, and led to the decreased contents of H₂O₂ and MDA, and the inhibition of O₂^- generation. Meanwhile, β-carotene accumulation promoted endogenous ABA biosynthesis of transgenic Arabidopsis, which resulted in ABA-induced stomatal closing and delayed water loss. In addition, ectopic expression of AgMYB5 increased expression levels of AtERD1, AtP5CS1, AtRD22, and AtRD29.

CONCLUSIONS

The findings indicated that AgMYB5 up-regulated β-carotene biosynthesis and drought tolerance of Arabidopsis.

Transcriptome profiling of celery petiole tissues reveals peculiarities of the collenchyma cell wall formation

Transcriptome and biochemical analyses are applied to individual plant cell types to reveal potential players involved in the molecular machinery of cell wall formation in specialized cells such as collenchyma. Plant collenchyma is a mechanical tissue characterized by an irregular, thickened cell wall and the ability to support cell elongation. The composition of the collenchyma cell wall resembles that of the primary cell wall and includes cellulose, xyloglucan, and pectin; lignin is absent. Thus, the processes associated with the formation of the primary cell wall in the collenchyma can be more pronounced compared to other tissues due to its thickening. Primary cell walls intrinsic to different tissues may differ in structure and composition, which should be reflected at the transcriptomic level. For the first time, we conducted transcriptome profiling of collenchyma strands isolated from young celery petioles and compared them with other tissues, such as parenchyma and vascular bundles. Genes encoding proteins involved in the primary cell wall formation during cell elongation, such as xyloglucan endotransglucosylase/hydrolases, expansins, and leucine-rich repeat proteins, were significantly activated in the collenchyma. As the key players in the transcriptome orchestra of collenchyma, xyloglucan endotransglucosylase/hydrolase transcripts were characterized in more detail, including phylogeny and expression patterns. The comprehensive approach that included transcriptome and biochemical analyses allowed us to reveal peculiarities of collenchyma cell wall formation and modification, matching the abundance of upregulated transcripts and their potential substrates for revealed gene products. As a result, specific isoforms of multigene families were determined for further functional investigation.

Origin, evolution, breeding, and omics of Apiaceae: a family of vegetables and medicinal plants

Many of the world's most important vegetables and medicinal crops, including carrot, celery, coriander, fennel, and cumin, belong to the Apiaceae family. In this review, we summarize the complex origins of Apiaceae and the current state of research on the family, including traditional and molecular breeding practices, bioactive compounds, medicinal applications, nanotechnology, and omics research. Numerous molecular markers, regulatory factors, and functional genes have been discovered, studied, and applied to improve vegetable and medicinal crops in Apiaceae. In addition, current trends in Apiaceae application and research are also briefly described, including mining new functional genes and metabolites using omics research, identifying new genetic variants associated with important agronomic traits by population genetics analysis and GWAS, applying genetic transformation, the CRISPR-Cas9 gene editing system, and nanotechnology. This review provides a reference for basic and applied research on Apiaceae vegetable and medicinal plants.

Combined Analysis of the Metabolome and Transcriptome to Explore Heat Stress Responses and Adaptation Mechanisms in Celery (Apium graveolens L.)

Celery is an important leafy vegetable that can grow during the cool season and does not tolerate high temperatures. Heat stress is widely acknowledged as one of the main abiotic stresses affecting the growth and yield of celery. The morphological and physiological indices of celery were investigated in the present study to explore the physiological mechanisms in response to high temperatures. Results showed that the antioxidant enzyme activity, proline, relative conductivity, and malondialdehyde were increased, while chlorophyll and the water content of leaves decreased under high-temperature conditions. Short-term heat treatment increased the stomatal conductance to cool off the leaves by transpiration; however, long-term heat treatment led to stomatal closure to prevent leaf dehydration. In addition, high temperature caused a disordered arrangement of palisade tissue and a loose arrangement of spongy tissue in celery leaves. Combined metabolomic and transcriptomic analyses were further used to reveal the regulatory mechanisms in response to heat stress at the molecular level in celery. A total of 1003 differential metabolites were identified and significantly enriched in amino acid metabolism and the tricarboxilic acid (TCA) cycle. Transcriptome sequencing detected 24,264 different genes, including multiple transcription factor families such as HSF, WRKY, MYB, AP2, bZIP, and bHLH family members that were significantly upregulated in response to heat stress, suggesting that these genes were involved in the response to heat stress. In addition, transcriptional and metabolic pathway analyses showed that heat stress inhibited the glycolysis pathway and delayed the TCA cycle but increased the expression of most amino acid synthesis pathways such as proline, arginine, and serine, consistent with the results of physiological indicators. qRT-PCR further showed that the expression pattern was similar to the expression abundance in the transcriptome. The important metabolites and genes in celery that significantly contributed to the response to high temperatures were identified in the present study, which provided the theoretical basis for breeding heat-resistant celery.

Integrative genome, transcriptome, microRNA, and degradome analysis of water dropwort (Oenanthe javanica) in response to water stress

Water dropwort (Liyang Baiqin, Oenanthe javanica (BI.) DC.) is an aquatic perennial plant from the Apiaceae family with abundant protein, dietary fiber, vitamins, and minerals. It usually grows in wet soils and can even grow in water. Here, whole-genome sequencing of O. javanica via HiSeq 2000 sequencing technology was reported for the first time. The genome size was 1.28 Gb, including 42,270 genes, of which 93.92% could be functionally annotated. An online database of the whole-genome sequences of water dropwort, Water dropwortDB, was established to share the results and facilitate further research on O. javanica (database homepage: http://apiaceae.njau.edu.cn/waterdropwortdb ). Water dropwortDB offers whole-genome and transcriptome sequences and a Basic Local Alignment Search Tool. Comparative analysis with other species showed that the evolutionary relationship between O. javanica and Daucus carota was the closest. Twenty-five gene families of O. javanica were found to be expanded, and some genetic factors (such as genes and miRNAs) related to phenotypic and anatomic differentiation in O. javanica under different water conditions were further investigated. Two miRNA and target gene pairs (miR408 and Oja15472, miR171 and Oja47040) were remarkably regulated by water stress. The obtained reference genome of O. javanica provides important information for future work, thus making in-depth genetic breeding and gene editing possible. The present study also provides a foundation for the understanding of the O. javanica response to water stress, including morphological, anatomical, and genetic differentiation.

Impact of Genomic and Transcriptomic Resources on Apiaceae Crop Breeding Strategies

The Apiaceae taxon is one of the most important families of flowering plants and includes thousands of species used for food, flavoring, fragrance, medical and industrial purposes. This study had the specific intent of reviewing the main genomics and transcriptomic data available for this family and their use for the constitution of new varieties. This was achieved starting from the description of the main reproductive systems and barriers, with particular reference to cytoplasmic (CMS) and nuclear (NMS) male sterility. We found that CMS and NMS systems have been discovered and successfully exploited for the development of varieties only in Foeniculum vulgare, Daucus carota, Apium graveolens and Pastinaca sativa; whereas, strategies to limit self-pollination have been poorly considered. Since the constitution of new varieties benefits from the synergistic use of marker-assisted breeding in combination with conventional breeding schemes, we also analyzed and discussed the available SNP and SSR marker datasets (20 species) and genomes (8 species). Furthermore, the RNA-seq studies aimed at elucidating key pathways in stress tolerance or biosynthesis of the metabolites of interest were limited and proportional to the economic weight of each species. Finally, by aligning 53 plastid genomes from as many species as possible, we demonstrated the precision offered by the super barcoding approach to reconstruct the phylogenetic relationships of Apiaceae species. Overall, despite the impressive size of this family, we documented an evident lack of molecular data, especially because genomic and transcriptomic resources are circumscribed to a small number of species. We believe that our contribution can help future studies aimed at developing molecular tools for boosting breeding programs in crop plants of the Apiaceae family.

The Accumulation of Lutein and β-Carotene and Transcript Profiling of Genes Related to Carotenoids Biosynthesis in Yellow Celery

Carotenoids are the general term of natural pigments. The formation of plant color is probably related to the components of carotenoids. As the yellow variety of celery, it is rich in the composition and content of carotenoids. However, the transcript profiling and roles of the genes related to carotenoids biosynthesis in yellow celery remain unclear. In this study, three yellow celery cultivars at different growth stages were used to analyze the content and composition of carotenoids and transcriptional changes of carotenoid biosynthesis-related genes. The lutein and β-carotene were detected in yellow celery cultivar, while α-carotene and lycopene were not detected. The contents of lutein and β-carotene were higher in leaf blades than in petioles. During the growth and development, the contents of lutein and β-carotene gradually decreased in celery. Compared with the other two cultivars, the contents of lutein and β-carotene were the highest in 'Huangtaiji' of 65 days after sowing (DAS) and 85 DAS and 'Liuhehuangxinqin' of 105 DAS, respectively. The expression levels of AgLCYB and AgPSY2 genes were significantly correlated with lutein and β-carotene contents. This work provided a reference for the further study on carotenoid metabolisms in yellow celery and also made sense on the way of cultivating yellow celery with high carotenoids content.

The celery genome sequence reveals sequential paleo-polyploidizations, karyotype evolution and resistance gene reduction in apiales

Celery (Apium graveolens L. 2n = 2x = 22), a member of the Apiaceae family, is among the most important and globally grown vegetables. Here, we report a high-quality genome sequence assembly, anchored to 11 chromosomes, with total length of 3.33 Gb and N50 scaffold length of 289.78 Mb. Most (92.91%) of the genome is composed of repetitive sequences, with 62.12% of 31 326 annotated genes confined to the terminal 20% of chromosomes. Simultaneous bursts of shared long-terminal repeats (LTRs) in different Apiaceae plants suggest inter-specific exchanges. Two ancestral polyploidizations were inferred, one shared by Apiales taxa and the other confined to Apiaceae. We reconstructed 8 Apiales proto-chromosomes, inferring their evolutionary trajectories from the eudicot common ancestor to extant plants. Transcriptome sequencing in three tissues (roots, leaves and petioles), and varieties with different-coloured petioles, revealed 4 and 2 key genes in pathways regulating anthocyanin and coumarin biosynthesis, respectively. A remarkable paucity of NBS disease-resistant genes in celery (62) and other Apiales was explained by extensive loss and limited production of these genes during the last ~10 million years, raising questions about their biotic defence mechanisms and motivating research into effects of chemicals, for example coumarins, that give off distinctive odours. Celery genome sequencing and annotation facilitates further research into important gene functions and breeding, and comparative genomic analyses in Apiales.

The gene encoding lycopene epsilon cyclase of celery enhanced lutein and β-carotene contents and confers increased salt tolerance in Arabidopsis

Celery (Apium graveolens L.) is a leafy vegetable of Apiaceae, which is greatly popular because of its rich nutrients. Lutein and β-carotene are two important carotenoids. Lycopene epsilon cyclase (LCY-ε) is a key branch point enzyme in the carotenoid biosynthetic pathway. In this study, we cloned the AgLCY-ε gene from celery and overexpressed it in Arabidopsis. The results showed that both lutein and β-carotene accumulation increased significantly in transgenic Arabidopsis hosting AgLCY-ε gene, compared with wild type (WT) plants. The transcription levels of AtPSY and AtCRTISO genes involved in carotenoids biosynthesis also increased in transgenic lines. One-month-old transgenic Arabidopsis seedlings were treated with 200 mM NaCl. The malondialdehyde (MDA) content in transgenic Arabidopsis plants after salt treatment was significantly lower, and the activities of the two antioxidant enzymes, superoxide dismutase (SOD) and peroxidase (POD), were significantly increased than that of WT plants. Overexpression of AgLCY-ε gene showed increased lutein and β-carotene accumulations, and enhanced salt tolerance in transgenic plants.

AgNAC1, a celery transcription factor, related to regulation on lignin biosynthesis and salt tolerance

The NAC transcription factor participates in various biotic and abiotic stress responses and plays a critical role in plant development. Lignin is a water-insoluble dietary fiber, but it is second only to cellulose in abundance. Celery is the main source of dietary fiber, but its quality and production are limited by various abiotic stresses. Here, AgNAC1 containing the NAM domain was identified from celery. AgNAC1 was found to be a nuclear protein. Transgenic Arabidopsis thaliana plants hosting AgNAC1 have longer root lengths and stomatal axis lengths than the wide type (WT). The evidence from lignin determination and expression levels of lignin-related genes indicated that AgNAC1 plays a vital role in lignin biosynthesis. Furthermore, the results of the physiological characterization and the drought and salt treatments indicate that AgNAC1-overexpressing plants are significantly resistive to salt stress. Under drought and salt treatments, the AgNAC1 transgenic Arabidopsis thaliana plants presented increased superoxide dismutase (SOD) and peroxidase (POD) activities and decreased malondialdehyde (MDA) content and size of stomatal apertures relatively to the WT plants. The AgNAC1 served as a positive regulator in inducing the expression of stress-responsive genes. Overall, the overexpressing AgNAC1 enhanced the plants' resistance to salt stress and played a regulatory role in lignin accumulation.

Advances in AP2/ERF super-family transcription factors in plant

In the whole life process, many factors including external and internal factors affect plant growth and development. The morphogenesis, growth, and development of plants are controlled by genetic elements and are influenced by environmental stress. Transcription factors contain one or more specific DNA-binding domains, which are essential in the whole life cycle of higher plants. The AP2/ERF (APETALA2/ethylene-responsive element binding factors) transcription factors are a large group of factors that are mainly found in plants. The transcription factors of this family serve as important regulators in many biological and physiological processes, such as plant morphogenesis, responsive mechanisms to various stresses, hormone signal transduction, and metabolite regulation. In this review, we summarized the advances in identification, classification, function, regulatory mechanisms, and the evolution of AP2/ERF transcription factors in plants. AP2/ERF family factors are mainly classified into four major subfamilies: DREB (Dehydration Responsive Element-Binding), ERF (Ethylene-Responsive-Element-Binding protein), AP2 (APETALA2) and RAV (Related to ABI3/VP), and Soloists (few unclassified factors). The review summarized the reports about multiple regulatory functions of AP2/ERF transcription factors in plants. In addition to growth regulation and stress responses, the regulatory functions of AP2/ERF in plant metabolite biosynthesis have been described. We also discussed the roles of AP2/ERF transcription factors in different phytohormone-mediated signaling pathways in plants. Genomic-wide analysis indicated that AP2/ERF transcription factors were highly conserved during plant evolution. Some public databases containing the information of AP2/ERF have been introduced. The studies of AP2/ERF factors will provide important bases for plant regulatory mechanisms and molecular breeding.

Transcriptome analysis reveals immune-related gene expression changes with age in giant panda (Ailuropoda melanoleuca) blood

The giant panda (Ailuropoda melanoleuca), an endangered species endemic to western China, has long been threatened with extinction that is exacerbated by highly contagious and fatal diseases. Aging is the most well-defined risk factor for diseases and is associated with a decline in immune function leading to increased susceptibility to infection and reduced response to vaccination. Therefore, this study aimed to determine which genes and pathways show differential expression with age in blood tissues. We obtained 210 differentially expressed genes by RNA-seq, including 146 up-regulated and 64 down-regulated genes in old pandas (18-21yrs) compared to young pandas (2-6yrs). We identified ISG15, STAT1, IRF7 and DDX58 as the hub genes in the protein-protein interaction network. All of these genes were up-regulated with age and played important roles in response to pathogen invasion. Functional enrichment analysis indicated that up-regulated genes were mainly involved in innate immune response, while the down-regulated genes were mainly related to B cell activation. These may suggest that the innate immunity is relatively well preserved to compensate for the decline in the adaptive immune function. In conclusion, our findings will provide a foundation for future studies on the molecular mechanisms underlying immune changes associated with ageing.

Identification and Selection of Reference Genes for Quantitative Transcript Analysis in Corydalis yanhusuo

Corydalis yanhusuo is a medicinal plant frequently used in traditional Chinese medicine, which has effective medical effects in many aspects. Real-time polymerase chain reaction (RT-PCR) has been one of the most widely used methods in biosynthesis research due to its high sensitivity and quantitative properties in gene expression analysis. To obtain accurate normalization, reference genes are often selected in advance; however, no reference genes are available in C. yanhusuo. Herein, 12 reference gene candidates, named cyclophilin 2 (CYP2), elongation factor 1-α (EF1-α), protein phosphatase 2 (PP2A), SAND protein family (SAND), polypyrimidine tract-binding protein (PTBP), TIP41-like protein (TIP41), lyceraldehyde-3-phosphate hydrogenase (GAPDH), ubiquitin-conjugating enzyme 9 (UBC9), cyclophilin 1 (CYP1), tubulin beta (TUBA), thioredoxin (YLS8), and polyubiquitin 10 (UBQ10), were selected for stability analysis. After being treated with hormone, UV, salt, metal, oxidative, drought, cold (4 °C), and hot stresses (40 °C), the qRT-PCR data of the selected genes was analyzed with NormFinder, geNorm, and BestKeeper. The result indicated that GAPDH, SNAD, and PP2A were the top three most stable reference genes under most treatments. This study selected and validated reliable reference genes in C. yanhusuo under various environmental conditions, which can provide great help for future research on gene expression normalization in C. yanhusuo.

Celery and Celeriac: A Critical View on Present and Future Breeding

Cultivated for the crispy petioles and round, fleshy, and flavored hypocotyl celery and celeriac have over two centuries of breeding history in Europe. In this review paper we summarized the most recent advances touching when necessary the historical context of celery and celeriac breeding. In the post genomic era of research, the genome sequence of celery is only partially available. We comprised however in this paper the most important aspects of celery genetics that are available today and have applicability in celery modern cultivars development. We discussed the problems and traits that drive the main celery and celeriac breeding goals, like hybrid seed production, disease resistance, and interesting enlarged hypocotyl and petiole characteristics. Besides the classical breeding traits we covered the potential of integration of existing cultivars as sources for consumer oriented traits like nutraceuticals and health promoting substances. Sustainability is a subject that is continuously growing in popularity and we looked at the genetic base of celery and celeriac that makes them sources for abiotic stress resistance and candidates for phytoremediation. We explored the fundamental concepts gained in various fields of celery and related species research, as resources for future improvement of celery and celeriac germplasm. We forecast what the next years will bring to Apium breeding.

The genome sequence of celery (Apium graveolens L.), an important leaf vegetable crop rich in apigenin in the Apiaceae family

Celery (Apium graveolens L.) is a vegetable crop in the Apiaceae family that is widely cultivated and consumed because it contains necessary nutrients and multiple biologically active ingredients, such as apigenin and terpenoids. Here, we report the genome sequence of celery based on the use of HiSeq 2000 sequencing technology to obtain 600.8 Gb of data, achieving ~189-fold genome coverage, from 68 sequencing libraries with different insert sizes ranging from 180 bp to 10 kb in length. The assembled genome has a total sequence length of 2.21 Gb and consists of 34,277 predicted genes. Repetitive DNA sequences represent 68.88% of the genome sequences, and LTR retrotransposons are the main components of the repetitive sequences. Evolutionary analysis showed that a recent whole-genome duplication event may have occurred in celery, which could have contributed to its large genome size. The genome sequence of celery allowed us to identify agronomically important genes involved in disease resistance, flavonoid biosynthesis, terpenoid metabolism, and other important cellular processes. The comparative analysis of apigenin biosynthesis genes among species might explain the high apigenin content of celery. The whole-genome sequences of celery have been deposited at CeleryDB (http://apiaceae.njau.edu.cn/celerydb). The availability of the celery genome data advances our knowledge of the genetic evolution of celery and will contribute to further biological research and breeding in celery as well as other Apiaceae plants.

Constructing a Core Collection of the Medicinal Plant Angelica biserrata Using Genetic and Metabolic Data

Angelica biserrata is an important medicinal plant in Chinese traditional medicine. Its roots, which are known as Duhuo in Chinese, are broadly applied to treat inflammation, arthritis, and headache. With increasing market demand, the wild resources of A. biserrata have been overexploited, and conservation, assessment of genetic resources and breeding for this species is needed. Here, we sequenced the transcriptome of A. biserrata and developed simple sequence repeat (SSR) markers from it to construct a core collection based on 208 samples collected from Changyang-related regions. A total of 132 alleles were obtained for 17 SSR loci used with the polymorphic information content (PIC) ranging from 0.44 to 0.83. Abundant genetic diversity was inferred by Shannon's information index (1.51), observed (0.57) and expected heterozygosity (0.72). The clustering analysis resulted into two sample groups and analysis of molecular variance (AMOVA) showed only 6% genetic variation existed among populations. A further metabolic analysis of these samples revealed the main coumarin contents, such as osthole and columbianadin. According to the genetic and metabolic data, we adopted the least distance stepwise sampling strategy to construct seven preliminary core collections, of which the 20CC collection, which possessed 42 A. biserrata individuals accounting for 90.20% of the genetic diversity of the original germplasm, represented the best core collection. This study will contribute to the conservation and management of A. biserrata wild germplasm resources and provide a material basis for future selection and breeding of this medicinal plant.

Comparative de novo transcriptomics and untargeted metabolomic analyses elucidate complicated mechanisms regulating celery (Apium graveolens L.) responses to selenium stimuli

Presently, concern regarding the effects of selenium (Se) on the environment and organisms worldwide is increasing. Too much Se in the soil is harmful to plants. In this study, Illumina RNA sequencing and the untargeted metabolome of control and Se-treated celery seedlings were analyzed. In total, 297,911,046 clean reads were obtained and assembled into 150,218 transcripts (50,876 unigenes). A total of 36,287 unigenes were annotated using different databases. Additionally, 8,907 differentially expressed genes, including 5,319 up- and 3,588 downregulated genes, were identified between mock and Se-treated plants. “Phenylpropanoid biosynthesis” was the most enriched KEGG pathway. A total of 24 sulfur and selenocompound metabolic unigenes were differentially expressed. Furthermore, 1,774 metabolites and 237 significant differentially accumulated metabolites were identified using the untargeted metabolomic approach. We conducted correlation analyses of enriched KEGG pathways of differentially expressed genes and accumulated metabolites. Our findings suggested that candidate genes and metabolites involved in important biological pathways may regulate Se tolerance in celery. The results increase our understanding of the molecular mechanism responsible for celery’s adaptation to Se stress.

Isolation, purification and characterization of an ascorbate peroxidase from celery and overexpression of the AgAPX1 gene enhanced ascorbate content and drought tolerance in Arabidopsis

BACKGROUND

Celery is a widely cultivated vegetable abundant in ascorbate (AsA), a natural plant antioxidant capable of scavenging free radicals generated by abiotic stress in plants. Ascorbate peroxidase (APX) is a plant antioxidant enzyme that is important in the synthesis of AsA and scavenging of excess hydrogen peroxide. However, the characteristics and functions of APX in celery remain unclear to date.

RESULTS

In this study, a gene encoding APX was cloned from celery and named AgAPX1. The transcription level of the AgAPX1 gene was significantly upregulated under drought stress. AgAPX1 was expressed in Escherichia coli BL21 (DE3) and purified. The predicted molecular mass of rAgAPX1 was 33.16 kDa, which was verified by SDS-PAGE assay. The optimum pH and temperature for rAgAPX1 were 7.0 and 55 °C, respectively. Transgenic Arabidopsis hosting the AgAPX1 gene showed elevated AsA content, antioxidant capacity and drought resistance. Less decrease in net photosynthetic rate, chlorophyll content, and relative water content contributed to the high survival rate of transgenic Arabidopsis lines after drought.

CONCLUSIONS

The characteristics of APX in celery were different from that in other species. The enhanced drought resistance of overexpressing AgAPX1 in Arabidopsis may be achieved by increasing the accumulation of AsA, enhancing the activities of various antioxidant enzymes, and promoting stomatal closure. Our work provides new evidence to understand APX and its response mechanisms to drought stress in celery.

Elevated gibberellin altered morphology, anatomical structure, and transcriptional regulatory networks of hormones in celery leaves

Gibberellins (GAs), as one of the important hormones in regulating the growth and development of higher plants, can significantly promote cell elongation and expansion. Celery is a widely grown leafy vegetable crop with rich nutritional value. However, the effect of gibberellins on celery leaves is unclear. In this paper, the celery variety "Jinnan Shiqin" plants were treated with gibberellic acid (GA₃) and paclobutrazol (PBZ, a gibberellin inhibitor). Our results showed that GA₃ treatment promoted the growth of celery leaves and caused lignification of celery leaf tissue. In addition, the transcript levels of genes associated with gibberellins, auxin, cytokinins, ethylene, jasmonic acid, abscisic acid, and brassinolide were altered in response to increased or decreased exogenous gibberellins or inhibitor. GA₃ may regulate celery growth by interacting with other hormones through crosstalk mechanisms. This study provided a reference for further study of the regulation mechanism of gibberellins metabolism, and exerted effects on understanding the role of gibberellins in the growth and development of celery.

Selection of appropriate reference genes for RT-qPCR analysis under abiotic stress and hormone treatment in celery

Celery is one of the most important vegetable crop and its yield and quality is influenced by many environmental factors. Researches on gene expression not only help to unravel the molecular regulatory mechanism but also identify the key genes in the biological response. RT-qPCR is a commonly used technology to quantify the gene expression. Selecting an appropriate reference gene is an effective approach to improve the accuracy of RT-qPCR assay. To our knowledge, the evaluation of reference genes under different treatments in celery has not been reported yet. In this study, the expression stabilities of eight candidate reference genes (ACTIN, eIF-4α , GAPDH, TBP, TUB-A, UBC, TUB-B, and EF-1α ) under abiotic stresses (heat, cold, drought, and salt) and hormone treatments (SA, MeJA, GA, and ABA) were detected. The expression stabilities of candidate genes were compared and ranked by geNorm, NormFinder, BestKeeper, ΔCt, and RefFinder programs. The results calculated by different programs were not completely consistent. Considering the comprehensive analysis results, ACTIN was the most stable reference gene and TUB-B showed the worst expression stabilities under the selected abiotic stress and hormone treatments in celery. The reliability of reference genes was further confirmed by the normalization of CAT1 gene under drought stress. This study presented evidences and basis to select the appropriate reference genes under different treatments in celery.

Genomic identification of AP2/ERF transcription factors and functional characterization of two cold resistance-related AP2/ERF genes in celery (Apium graveolens L.)

This study analyzed the AP2/ERF transcription factors in celery and showed that two dehydration-responsive-element-binding (DREB) transcription factors, AgDREB1 and AgDREB2, contribute to the enhanced resistance to abiotic stress in transgenic Arabidopsis. The AP2/ERF family is a large family of transcription factors (TFs) in higher plants that plays a central role in plant growth, development, and response to environmental stress. Here, 209 AP2/ERF family members were identified in celery based on genomic and transcriptomic data. The TFs were classified into four subfamilies (i.e., DREB, ERF, RAV, and AP2) and Soloist. Evolution analysis indicated that the AP2/ERF TFs are ancient molecules and have expanded in the long-term evolution process of plants and whole-genome duplication events. AgAP2/ERF proteins may be associated with multiple biological processes as predicted by the interaction network. The expression profiles and sequence alignment analysis of the TFs in the DREB-A1 group showed that eight genes could be divided into four branches. Two genes, AgDREB1 and AgDREB2, from the DREB-A1 group were selected for further analysis. Subcellular localization assay suggested that the two proteins are nuclear proteins. Yeast one hybrid assay demonstrated that the two proteins could bind to the dehydration-responsive element (DRE). The overexpression of AgDREB1 and AgDREB2 in Arabidopsis induced the increased tolerance to cold treatment and the up-regulation of the COR genes expression. AgDREB1 and AgDREB2 might function as transcriptional activators in regulating the downstream genes by binding to corresponding DRE to enhance stress tolerance in celery.

Advances in research on the carrot, an important root vegetable in the Apiaceae family

Carrots (Daucus carota L.), among the most important root vegetables in the Apiaceae family, are cultivated worldwide. The storage root is widely utilized due to its richness in carotenoids, anthocyanins, dietary fiber, vitamins and other nutrients. Carrot extracts, which serve as sources of antioxidants, have important functions in preventing many diseases. The biosynthesis, metabolism, and medicinal properties of carotenoids in carrots have been widely studied. Research on hormone regulation in the growth and development of carrots has also been widely performed. Recently, with the development of high-throughput sequencing technology, many efficient tools have been adopted in carrot research. A large amount of sequence data has been produced and applied to improve carrot breeding. A genome editing system based on CRISPR/Cas9 was also constructed for carrot research. In this review, we will briefly summarize the origins, genetic breeding, resistance breeding, genome editing, omics research, hormone regulation, and nutritional composition of carrots. Perspectives about future research work on carrots are also briefly provided.

Transcriptome profiling reveals the association of multiple genes and pathways contributing to hormonal control in celery leaves

Celery is a vital vegetable belonging to the Apiaceae family. The leaves of celery are its main edible part with high nutritional value. Hormone signaling plays diverse and critical roles in controlling plant growth and development. However, the molecular mechanism of hormone regulating growth and development in celery leaves has not been investigated. Here, we aimed to understand the molecular functions of genes related to hormone metabolism in celery leaf growth and development. A total of 77 hormone-related differentially expressed genes (DEGs) were identified from the transcriptome of celery leaves at three development stages. The roles and interactions of DEGs in the growth and development of celery leaves were discussed. The contents of multiple hormones (IAA, ZR, ABA, BR, GA3, and MeJA) in celery leaf development were also detected. The changes of endogenous hormone level during the development of celery leaves could be regulated by the expressions of hormone-related genes. Our results indicated that the plant hormones had a complex regulatory mechanism for the growth of celery leaves. Our current findings will provide potential valuable references for the future research on celery leaf development.

Genome-wide mining seed-specific candidate genes from peanut for promoter cloning

Peanut seeds are ideal bioreactors for the production of foreign recombinant proteins and/or nutrient metabolites. Seed-Specific Promoters (SSPs) are important molecular tools for bioreactor research. However, few SSPs have been characterized in peanut seeds. The mining of Seed-Specific Candidate Genes (SSCGs) is a prerequisite for promoter cloning. Here, we described an approach for the genome-wide mining of SSCGs via comparative gene expression between seed and nonseed tissues. Three hundred thirty-seven SSCGs were ultimately identified, and the top 108 SSCGs were characterized. Gene Ontology (GO) analysis revealed that some SSCGs were involved in seed development, allergens, seed storage and fatty acid metabolism. RY REPEAT and GCN4 motifs, which are commonly found in SSPs, were dispersed throughout most of the promoters of SSCGs. Expression pattern analysis revealed that all 108 SSCGs were expressed specifically or preferentially in the seed. These results indicated that the promoters of the 108 SSCGs may perform functions in a seed-specific and/or seed-preferential manner. Moreover, a novel SSP was cloned and characterized from a paralogous gene of SSCG29 from cultivated peanut. Together with the previously characterized SSP of the SSCG5 paralogous gene in cultivated peanut, these results implied that the method for SSCG identification in this study was feasible and accurate. The SSCGs identified in this work could be widely applied to SSP cloning by other researchers. Additionally, this study identified a low-cost, high-throughput approach for exploring tissue-specific genes in other crop species.

The Population Genetics of Alternaria tenuissima in Four Regions of China as Determined by Microsatellite Markers Obtained by Transcriptome Sequencing

A total of 32,284 unigenes were obtained from the transcriptome of Alternaria tenuissima, a pathogenic fungus causing foliar disease in tomato, using next-generation sequencing (NGS) technology. In total, 24,670 unigenes were annotated using five databases, including NCBI non-redundant protein, Swiss-Prot, euKaryotic Orthologous Groups, Kyoto Encyclopedia of Genes and Genomes, and the Gene Ontology. A total of 1,140 simple sequence repeats were also identified for use as molecular markers. Sixteen of the simple sequence repeat loci were selected to study the population structure of A. tenuissima. A population genetic analysis of 191 A. tenuissima isolates, sampled from four geographic regions in China, indicated that A. tenuissima had a high level of genetic diversity, and that the selected simple sequence repeat markers could reliably capture the genetic variation. The null hypothesis of random mating was rejected for all four geographic regions in China. Isolation by distance was observed for the entire data set, but not within clusters, which is indicative of barriers to gene flow among geographic regions. The analyses of Bayesian and principal coordinates, however, did not separate four geographic regions into four separate genetic clusters. The different levels of historical migration rates suggest that isolation by distance did not represent a major biological obstacle to the spread of A. tenuissima. The potential epidemic spread of A. tenuissima in China may occur through the transport of plant products or other factors. The presented results provide a basis for a comprehensive understanding of the population genetics of A. tenuissima in China.

Exogenous brassinosteroids altered cell length, gibberellin content, and cellulose deposition in promoting carrot petiole elongation

Brassinosteroid (BR) is a predominant plant hormone in regulating cell elongation and cell size. BR-deficient mutants display reduced plant growth and dwarfism in Arabidopsis and rice. In carrot, BRs promote petiole elongation, but its underlying mechanism involving exogenous BR remains unknown. Here, weighted gene co-expression network analysis and promoter region analysis were adopted to identify the potential genes that interacted with DcBZR1/BES1. Bioactive gibberellin (GA) level and cellulose deposition were also determined in the control and treated plants. Quantitative real-time PCR was performed to detect the expression profiles of GA biosynthesis-related genes, GA signaling genes, and cellulose synthase genes. Bioactive GA level and cellulose deposition were upregulated after the petioles were treated with 24-epibrassinolide (24-EBL). The most putative DcBZR1/BES1 genes were clustered in yellow module. The expression level of DCAR_009411 (a GA5-like gene) was significantly induced after 3 h of treatment. The expression levels of DCAR_019754 and DCAR_013973 (CESA-like genes) were also significantly induced after 3 h of 24-EBL treatment. Our results suggested that the effect of BR on carrot petiole growth was quick. These results also provided potential insights into the mechanism by which BRs modulate GA and cellulose synthesis to promote cell elongation in carrot petioles.

AgMYB2 transcription factor is involved in the regulation of anthocyanin biosynthesis in purple celery (Apium graveolens L.)

This study showed that an R2R3-MYB transcription factor, AgMYB2, functions in anthocyanin biosynthesis and accumulation in purple celery. Anthocyanins are involved in tissue coloration and stress response in plants. Foods containing high anthocyanin content are also beneficial to human health. Purple celery accumulated amounts of anthocyanins in the petioles. The biosynthesis of anthocyanin in plants is mainly regulated by the R2R3-MYB transcription factor (TF). However, the R2R3-MYB TF that controls anthocyanin accumulation in purple celery remains unclear. In this study, an R2R3-MYB TF gene, AgMYB2, was cloned from purple celery and characterized as anthocyanin biosynthetic regulator. Sequence analysis indicated that AgMYB2 contained highly conserved R2R3 domain and two anthocyanin characteristic motifs, ANDV motif and KPRPR[S/T]F motif. The relative expression level of AgMYB2 in purple celery was significantly higher than that in non-purple celery at three developmental stages. Heterologous expression of AgMYB2 in Arabidopsis generated more anthocyanins and resulted in dark-purple leaves and flowers. The expression levels of anthocyanin biosynthetic genes and the antioxidant activity of transgenic Arabidopsis carrying AgMYB2 were up-regulated. The determination of anthocyanin glycosylation activity of Arabidopsis crude enzyme verified the anthocyanin biosynthesis regulatory function of AgMYB2 at the protein level. The interaction between AgMYB2 and bHLH proteins was shown by yeast two-hybrid assay. The results will help to elucidate the molecular mechanism of anthocyanin biosynthesis in purple celery and provide an approach for cultivating plants with high anthocyanin content.

Development of novel microsatellites for population genetic analysis of Phenacoccus solenopsis Tinsley (Hemipeta: Pseudoccoccidae) based on genomic analysis

The cotton mealybug, Phenacoccus solenopsis Tinsley (Hemipeta: Pseudoccoccidae), is an aggressively invasive pest causing huge economic losses of crops around the world. In this study, we developed genome-wide microsatellites for population genetic analysis of P. solenopsis. We obtained a random genome of P. solenopsis with a size of 267.07 Mb and scaffold N50 of 14.12 Kb. In total 115,639 microsatellites were isolated from the genome, of which those with trinucleotide motifs were the most abundant. Forty-two polymorphic loci were selected for primer validation based on three populations. Allele numbers varied from 2 to 5 with an average value of 2.5 per locus, and allelic richness ranged from 1.00 to 4.48. The observed heterozygosity (H₀) and expected heterozygosity (H_E) ranged from 0.00 to 0.92 and 0.00 to 0.73, respectively. Population genetic structure analysis based on the developed markers revealed strong differentiation between three populations of P. solenopsis collected from its invasive range in China. The microsatellites developed in our study should provide efficient genetic markers for population level studies of P. solenopsis to reveal invasion history and patterns of dispersal.

The leaf transcriptome of fennel (Foeniculum vulgare Mill.) enables characterization of the t-anethole pathway and the discovery of microsatellites and single-nucleotide variants

Fennel is a plant species of both agronomic and pharmaceutical interest that is characterized by a shortage of genetic and molecular data. Taking advantage of NGS technology, we sequenced and annotated the first fennel leaf transcriptome using material from four different lines and two different bioinformatic approaches: de novo and genome-guided transcriptome assembly. A reference transcriptome for assembly was produced by combining these two approaches. Among the 79,263 transcripts obtained, 47,775 were annotated using BLASTX analysis performed against the NR protein database subset with 11,853 transcripts representing putative full-length CDS. Bioinformatic analyses revealed 1,011 transcripts encoding transcription factors, mainly from the BHLH, MYB-related, C2H2, MYB, and ERF families, and 6,411 EST-SSR regions. Single-nucleotide variants of SNPs and indels were identified among the 8 samples at a frequency of 0.5 and 0.04 variants per Kb, respectively. Finally, the assembled transcripts were screened to identify genes related to the biosynthesis of t-anethole, a compound well-known for its nutraceutical and medical properties. For each of the 11 genes encoding structural enzymes in the t-anethole biosynthetic pathway, we identified at least one transcript showing a significant match. Overall, our work represents a treasure trove of information exploitable both for marker-assisted breeding and for in-depth studies on thousands of genes, including those involved in t-anethole biosynthesis.

Elevated CO2 induces alteration in lignin accumulation in celery (Apium graveolens L.)

Carbon dioxide (CO₂) is an important regulator of plant growth and development, and its proportion in the atmosphere continues to rise now. Lignin is one of the major secondary products in plants with vital biological functions. However, the relationship between CO₂ level and xylogenesis in celery is still unknown. In order to investigate the effects of increasing CO₂ concentration on lignin accumulation in celery, 'Jinnanshiqin' were exposed to two CO₂ applications, 400 (e₀) and 1000 μmol mol^-1 (e₁), respectively. Plant morphology and lignin distribution in celery plants treated with elevated CO₂ did not change significantly. There was an upward trend on lignin content in celery leaves, and the transcript abundance of 12 genes involved in lignin metabolism has altered in response to elevated CO₂. The effects of high level of CO₂ on different tissues were different. Our works confirmed that CO₂ may play an important role in lignin accumulation in celery leaves. The current study will offer new evidence to understand the regulation mechanism of lignin biosynthesis under elevated CO₂ and provide a reference to improve celery quality by adjusting the growth environment.

Isolation, purification, and characterization of AgUCGalT1, a galactosyltransferase involved in anthocyanin galactosylation in purple celery (Apium graveolens L.)

This study showed that a galactosyltransferase, AgUCGalT1, is involved in anthocyanin galactosylation in purple celery. Celery is a well-known vegetable because of its rich nutrients, low calories, and medicinal values. Its petioles and leaf blades are the main organs acting as nutrient sources. UDP-galactose: cyanidin 3-O-galactosyltransferase can transfer the galactosyl moiety from UDP-galactose to the 3-O-position of cyanidin through glycosylation. This process enhances the stability and water solubility of anthocyanins. In the present study, LC-MS data indicated that abundant cyanidin-based anthocyanins accumulated in the petioles of purple celery ('Nanxuan liuhe purple celery'). A gene encoding UDP-galactose: cyanidin 3-O-galactosyltransferase, namely AgUCGalT1, was isolated from purple celery and expressed in Escherichia coli BL21 (DE3). Sequence alignments revealed that the AgUCGalT1 protein contained a highly conserved putative secondary plant glycosyltransferase (PSPG) motif. The glycosylation product catalyzed by AgUCGalT1 was detected using UPLC equipment. The recombinant AgUCGalT1 had an optimal enzyme activity at 35 °C and pH 8.0, and showed highest enzyme activity toward cyanidin among the enzyme activities involving other substances, namely, peonidin, quercetin, and kaempferol. The expression levels of AgUCGalT1 were positively correlated with the total anthocyanin contents in purple and non-purple celery varieties. Crude enzymes extracted from purple celery exhibited glycosylation ability, whereas crude enzymes obtained from non-purple celery did not have this ability. This work provided evidence as a basis for investigations on the function of AgUCGalT1 in anthocyanin glycosylation in purple celery.

Transcriptome-wide identification and expression profile analysis of the bHLH family genes in Camellia sinensis

The tea plant is an important commercial horticulture crop cultivated worldwide. Yield and quality of this plant are influenced by abiotic stress. The bHLH family transcription factors play a pivotal role in the growth and development, including abiotic stress response, of plants. A growing number of bHLH proteins have been functionally characterized in plants. However, few studies have focused on the bHLH proteins in tea plants. In this study, 120 CsbHLH TFs were identified from tea plants using computational prediction method. Structural analysis detected 23 conservative residues, with over 50% identities in the bHLH domain. Moreover, 103 CsbHLH proteins were assumed to bind DNA and encompassed 98 E-Box binders and 85 G-Box binders. The CsbHLH proteins were grouped into 20 subfamilies based on phylogenetic analysis and a previous classification system. A survey of transcriptome profiling screened 22 and 39 CsbHLH genes that were upregulated under heat and drought stress. Nine CsbHLH genes were validated using qRT-PCR. Results were approximately in accordance with transcriptome data. These genes could be induced by one or more abiotic stresses.

Exploration and Exploitation of Novel SSR Markers for Candidate Transcription Factor Genes in Lilium Species

Lilies (Lilium sp.) are commercially important horticultural crops widely cultivated for their flowers and bulbs. Here, we conducted large-scale data mining of the lily transcriptome to develop transcription factor (TF)-associated microsatellite markers (TFSSRs). Among 216,768 unigenes extracted from our sequence data, 6966 unigenes harbored simple sequence repeats (SSRs). Seventy-one SSRs were associated with TF genes, and these were used to design primers and validate their potential as markers. These 71 SSRs were accomplished with 31 transcription factor families; including bHLH, MYB, C2H2, ERF, C3H, NAC, bZIP, and so on. Fourteen highly polymorphic SSRs were selected based on Polymorphic Information Content (PIC) values and used to study genetic diversity and population structure in lily accessions. Higher genetic diversity was observed in Longiflorum compared to Oriental and Asiatic populations. Lily accessions were divided into three sub-populations based in our structure analysis, and an un-rooted neighbor-joining tree effectively separated the accessions according to Asiatic, Oriental, and Longiflorum subgroups. Finally, we showed that 46 of the SSR-associated genes were differentially expressed in response to Botrytis elliptica infection. Thus, our newly developed TFSSR markers represent a powerful tool for large-scale genotyping, high-density and comparative mapping, marker-aided backcrossing, and molecular diversity analysis of Lilium sp.

Identification of Putative Precursor Genes for the Biosynthesis of Cannabinoid-Like Compound in Radula marginata

The liverwort Radula marginata belongs to the bryophyte division of land plants and is a prospective alternate source of cannabinoid-like compounds. However, mechanistic insights into the molecular pathways directing the synthesis of these cannabinoid-like compounds have been hindered due to the lack of genetic information. This prompted us to do deep sequencing, de novo assembly and annotation of R. marginata transcriptome, which resulted in the identification and validation of the genes for cannabinoid biosynthetic pathway. In total, we have identified 11,421 putative genes encoding 1,554 enzymes from 145 biosynthetic pathways. Interestingly, we have identified all the upstream genes of the central precursor of cannabinoid biosynthesis, cannabigerolic acid (CBGA), including its two first intermediates, stilbene acid (SA) and geranyl diphosphate (GPP). Expression of all these genes was validated using quantitative real-time PCR. We have characterized the protein structure of stilbene synthase (STS), which is considered as a homolog of olivetolic acid in R. marginata. Moreover, the metabolomics approach enabled us to identify CBGA-analogous compounds using electrospray ionization mass spectrometry (ESI-MS/MS) and gas chromatography mass spectrometry (GC-MS). Transcriptomic analysis revealed 1085 transcription factors (TF) from 39 families. Comparative analysis showed that six TF families have been uniquely predicted in R. marginata. In addition, the bioinformatics analysis predicted a large number of simple sequence repeats (SSRs) and non-coding RNAs (ncRNAs). Our results collectively provide mechanistic insights into the putative precursor genes for the biosynthesis of cannabinoid-like compounds and a novel transcriptomic resource for R. marginata. The large-scale transcriptomic resource generated in this study would further serve as a reference transcriptome to explore the Radulaceae family.

CeleryDB: a genomic database for celery

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.

Comparative proteomic analysis provides novel insights into chlorophyll biosynthesis in celery under temperature stress

Chlorophyll (Chl) is essential for light harvesting and energy transduction in photosynthesis. A proper amount of Chl within plant cells is important to celery (Apium graveolens) yield and quality. Temperature stress is an influential abiotic stress affecting Chl biosynthesis and plant growth. There are limited proteomic studies regarding Chl accumulation under temperature stress in celery leaves. Here, the proteins from celery leaves under different temperature treatments (4, 25 and 38°C) were analyzed using a proteomic approach. There were 71 proteins identified through MALDI-TOF-TOF analysis. The relative abundance of proteins involved in carbohydrate and energy metabolism, protein metabolism, amino acid metabolism, antioxidant and polyamine biosynthesis were enhanced under cold stress. These temperature stress-responsive proteins may establish a new homeostasis to enhance temperature tolerance. Magnesium chelatase (Mg-chelatase) and glutamate-1-semialdehyde aminotransferase (GSAT), related to Chl biosynthesis, showed increased abundances under cold stress. Meanwhile, the Chl contents were decreased in heat- and cold-stressed celery leaves. The inhibition of Chl biosynthesis may be due to the downregulated mRNA levels of 15 genes involved in Chl biosynthesis. The study will expand our knowledge on Chl biosynthesis and the temperature tolerance mechanisms in celery leaves.

Advances in the research of celery, an important Apiaceae vegetable crop

Celery (Apium graveolens L.), one of the most important vegetables in Apiaceae family, is cultivated worldwide and utilized in food and cosmetic industries because it is an excellent source of vitamins, phenolic compounds, volatile oils and other nutrients. Celery extracts possess various medicinal properties, such as antibacterial, anti-inflammatory and lowering blood glucose and serum lipid levels. With the rapid advancements in molecular biology and sequencing technology, studies on celery have been performed. Numerous molecular markers and regulatory genes have been discovered and applied to improve celery. Research advances, including genetic breeding, genomics research, function genes and chemical composition, regarding celery are reviewed in this paper. Further exploration and application trends are briefly described. This review provides a reference for basic and applied research on celery, an important Apiaceae vegetable crop.

Transcriptional profiling of genes involved in ascorbic acid biosynthesis, recycling, and degradation during three leaf developmental stages in celery

Ascorbic acid (AsA) is an important nutrient in the human body and performs various healthy functions. With considerable medicinal properties, celery (Apium graveolens L.) could be a good source of AsA for human health. However, the biosynthetic, recycling, and degradation pathways of AsA in celery have yet to be characterized. To study the metabolic pathways involved in AsA, the genes involved in AsA biosynthesis, recycling, and degradation were isolated from celery, and their expression profiles and AsA levels were analyzed in the leaf blades and petioles of two celery varieties at three different growth stages. AsA levels were higher in 'Ventura' compared with 'Liuhehuangxinqin' in both tissues possibly because of different transcription levels of genes, such as L-galactose dehydrogenase (GalDH), L-galactono-1,4-lactone dehydrogenase (GalLDH), and glutathione reductase (GR). Results revealed that the D-mannose/L-galactose pathway may be the predominant pathway in celery, and the D-galacturonic acid pathway appeared to contribute largely to AsA accumulation in petioles than in leaf blades in 'Liuhehuangxinqin.' AsA contents are regulated by complex regulatory mechanisms and vary at different growth stages, tissues, and varieties in celery. The results provide novel insights into AsA metabolic pathways in leaf during celery growth and development.

Temperature expression patterns of genes and their coexpression with LncRNAs revealed by RNA-Seq in non-heading Chinese cabbage

BACKGROUND

Non-heading Chinese cabbage (NHCC, Brassica rapa ssp. chinensis) is an important leaf vegetable grown worldwide. However, little is known about the molecular mechanisms underlying tolerance for extreme temperature in NHCC. The limited availability of NHCC genomic information has greatly hindered functional analysis and molecular breeding.

RESULTS

Here, we conduct comprehensive analyses of cold and heat treatments in NHCC using RNA-seq. Approximately 790 million paired-end reads representing 136,189 unigenes with N50 length of 1705 bp were obtained. Totally, 14,329 differentially expressed genes (DEGs) were detected. Among which, 10 DEGs were detected in all treatments, including 7 up-regulated and 3 down-regulated. The enrichment analyses showed 25 and 33 genes were enriched under cold and heat treatments, respectively. Additionally, 10,001 LncRNAs were identified, and 9,687 belonged to novel LncRNAs. The expression of miRNAs were more than that of pri-miRNAs and LncRNAs. Furthermore, we constructed a coexpression network for LncRNAs and miRNAs. It showed 67 and 192 genes were regulated by LncRNAs under cold and heat treatments, respectively. We constructed the flowchart for identifying LncRNAs of NHCC using transcriptome. Except conducting the de novo transcriptome analyses, we also compared these unigenes with the Chinese cabbage proteins. We identified several most important genes, and discussed their regulatory networks and crosstalk in cold and heat stresses.

CONCLUSIONS

We presented the first comprehensive characterization for NHCC crops and constructed the flowchart for identifying LncRNAs using transcriptome. Therefore, this study represents a fully characterized NHCC transcriptome, and provides a valuable resource for genetic and genomic studies under abiotic stress.

Transcriptome Analysis of Differentially Expressed Genes Provides Insight into Stolon Formation in Tulipa edulis

Tulipa edulis (Miq.) Baker is an important medicinal plant with a variety of anti-cancer properties. The stolon is one of the main asexual reproductive organs of T. edulis and possesses a unique morphology. To explore the molecular mechanism of stolon formation, we performed an RNA-seq analysis of the transcriptomes of stolons at three developmental stages. In the present study, 15.49 Gb of raw data were generated and assembled into 74,006 unigenes, and a total of 2,811 simple sequence repeats were detected in T. edulis. Among the three libraries of stolons at different developmental stages, there were 5,119 differentially expressed genes (DEGs). A functional annotation analysis based on sequence similarity queries of the GO, COG, KEGG databases showed that these DEGs were mainly involved in many physiological and biochemical processes, such as material and energy metabolism, hormone signaling, cell growth, and transcription regulation. In addition, quantitative real-time PCR analysis revealed that the expression patterns of the DEGs were consistent with the transcriptome data, which further supported a role for the DEGs in stolon formation. This study provides novel resources for future genetic and molecular studies in T. edulis.