De novo assembly, gene annotation and marker development using Illumina paired-end transcriptome sequences in celery (Apium graveolens L.)

Cell-Type Specific miRNA Regulatory Network Responses to ABA Stress Revealed by Time Series Transcriptional Atlases in Arabidopsis

In plants, microRNAs (miRNAs) participate in complex gene regulatory networks together with the transcription factors (TFs) in response to biotic and abiotic stresses. To date, analyses of miRNAs-induced transcriptome remodeling are at the whole plant or tissue levels. Here, Arabidopsis's ABA-induced single-cell RNA-seq (scRNA-seq) is performed at different stages of time points-early, middle, and late. Single-cell level primary miRNAs (pri-miRNAs) atlas supported the rapid, dynamic, and cell-type specific miRNA responses under ABA treatment. MiRNAs respond rapidly and prior to target gene expression dynamics, and these rapid response miRNAs are highly cell-type specific, especially in mesophyll and vascular cells. MiRNA-TF-mRNA regulation modules are identified by identifying miRNA-contained feed-forward loops (M-FFLs) in the regulatory network, and regulatory networks with M-FFLs have higher co-expression and clustering coefficient (CC) values than those without M-FFLs, suggesting the hub role of miRNAs in regulatory networks. The cell-type-specific M-FFLs are regulated by these hub miRNAs rather than TFs through sc-RNA-seq network analysis. MiR858a-FBH3-MYB module inhibited the expression of MYB63 and MYB20, which related to the formation of plant secondary wall and the production of lignin, through M-FFL specifically in vascular. These results can provide prominent insights into miRNAs' dynamic and cell-type-specific roles in plant development and stress responses.

Genetic variation and structure of endemic and endangered wild celery (Kelussia odoratissima Mozaff.) quantified using novel microsatellite markers developed by next-generation sequencing

Kelussia odoratissima Mozaff. (Apiaceae) is a native plant that has been traditionally consumed in Iran's food and pharmaceutical industries. Overharvesting of the taxon, especially at the beginning of the growing season, due to its considerable medicinal and economic value, is believed to be the main reason for the extirpating of this plant. The consequences of the severe anthropogenic impacts on the genetic diversity of populations are poorly known. In order to investigate the level of genetic variation and patterns of the genetic structure of K. odoratissima, we developed novel microsatellite markers using the 454 Roche next-generation sequencing (NGS) platform for the first time. Out of 1,165 microsatellite markers bioinformatically confirmed, twenty-five were tested, of which 23 were used to screen genetic variation across 12 natural populations. Our results showed that the average number of alleles per locus and the polymorphic information content (PIC) were 10.87 (range 7 to 27), and 0.81 (range 0.67 to 0.94), respectively. The mean observed and expected heterozygosities (± SD) across all populations were 0.80 ± 0.31 and 0.72 ± 0.14, respectively. The average pairwise FST among the populations was 0.37 (range 0.04 to 0.81). Bayesian and distance-based clustering, and principal coordinate analyses revealed at least four major genetic clusters. Although high level of structure can be explained by landscape topography and geographic distance, presence of admixed populations can be associated to seed or pollen dispersal. Contrary to expectations, the high level of genetic variation and lack of inbreeding suggest that overexploitation has not yet significantly purged the allelic variability within the natural populations in protected areas.

Genomics empowering conservation action and improvement of celery in the face of climate change

MAIN CONCLUSION

Integration of genomic approaches like whole genome sequencing, functional genomics, evolutionary genomics, and CRISPR/Cas9-based genome editing has accelerated the improvement of crop plants including leafy vegetables like celery in the face of climate change. The anthropogenic climate change is a real peril to the existence of life forms on our planet, including human and plant life. Climate change is predicted to be a significant threat to biodiversity and food security in the coming decades and is rapidly transforming global farming systems. To avoid the ghastly future in the face of climate change, the elucidation of shifts in the geographical range of plant species, species adaptation, and evolution is necessary for plant scientists to develop climate-resilient strategies. In the post-genomics era, the increasing availability of genomic resources and integration of multifaceted genomics elements is empowering biodiversity conservation action, restoration efforts, and identification of genomic regions adaptive to climate change. Genomics has accelerated the true characterization of crop wild relatives, genomic variations, and the development of climate-resilient varieties to ensure food security for 10 billion people by 2050. In this review, we have summarized the applications of multifaceted genomic tools, like conservation genomics, whole genome sequencing, functional genomics, genome editing, pangenomics, in the conservation and adaptation of plant species with a focus on celery, an aromatic and medicinal Apiaceae vegetable. We focus on how conservation scientists can utilize genomics and genomic data in conservation and improvement.

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.

Transcriptomics and Genomics Analysis Uncover the Differentially Expressed Chlorophyll and Carotenoid-Related Genes in Celery

Celery (Apium graveolens L.), a plant from Apiaceae, is one of the most important vegetables and is grown worldwide. Carotenoids can capture light energy and transfer it to chlorophyll, which plays a central role in photosynthesis. Here, by performing transcriptomics and genomics analysis, we identified and conducted a comprehensive analysis of chlorophyll and carotenoid-related genes in celery and six representative species. Significantly, different contents and gene expression patterns were found among three celery varieties. In total, 237 and 290 chlorophyll and carotenoid-related genes were identified in seven species. No notable gene expansion of chlorophyll biosynthesis was detected in examined species. However, the gene encoding ζ-carotene desaturase (ZDS) enzyme in carotenoid was expanded in celery. Comparative genomics and RNA-seq analyses revealed 16 and 5 key genes, respectively, regulating chlorophyll and carotenoid. An intriguing finding is that chlorophyll and carotenoid-related genes were coordinately regulated by transcriptional factors, which could be distinctively classified into positive- and negative-regulation groups. Six CONSTANS (CO)-like transcription factors co-regulated chlorophyll and carotenoid-related genes were identified in celery. In conclusion, this study provides new insights into the regulation of chlorophyll and carotenoid by transcription factors.

Single-Cell RNA Sequencing for Plant Research: Insights and Possible Benefits

In recent years, advances in single-cell RNA sequencing (scRNA-seq) technologies have continued to change our views on biological systems by increasing the spatiotemporal resolution of our analysis to single-cell resolution. Application of scRNA-seq to plants enables the comprehensive characterization of both common and rare cell types and cell states, uncovering new cell types and revealing how cell types relate to each other spatially and developmentally. This review provides an overview of scRNA-seq methodologies, highlights the application of scRNA-seq in plant science, justifies why scRNA-seq is a master player of sequencing, and explains the role of single-cell transcriptomics technologies in environmental stress adaptation, alongside the challenges and prospects of single-cell transcriptomics. Collectively, we put forward a central role of single-cell sequencing in plant research.

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.

Transcriptomic analysis to reveal the differentially expressed miRNA targets and their miRNAs in response to Ralstonia solanacearum in ginger species

BACKGROUND

Bacterial wilt is the most devastating disease in ginger caused by Ralstonia solanacearum. Even though ginger (Zingiber officinale) and mango ginger (Curcuma amada) are from the same family Zingiberaceae, the latter is resistant to R. solanacearum infection. MicroRNAs have been identified in many crops which regulates plant-pathogen interaction, either through silencing genes or by blocking mRNA translation. However, miRNA's vital role and its targets in mango ginger in protecting bacterial wilt is not yet studied extensively. In the present study, using the "psRNATarget" server, we analyzed available ginger (susceptible) and mango ginger (resistant) transcriptome to delineate and compare the microRNAs (miRNA) and their target genes (miRTGs).

RESULTS

A total of 4736 and 4485 differential expressed miRTGs (DEmiRTGs) were identified in ginger and mango ginger, respectively, in response to R. solanacearum. Functional annotation results showed that mango ginger had higher enrichment than ginger in top enriched GO terms. Among the DEmiRTGs, 2105 were common in ginger and mango ginger. However, 2337 miRTGs were expressed only in mango ginger which includes 62 defence related and upregulated miRTGs. We also identified 213 miRTGs upregulated in mango ginger but downregulated in ginger, out of which 23 DEmiRTGS were defence response related. We selected nine miRNA/miRTGs pairs from the data set of common miRTGs of ginger and mango ginger and validated using qPCR.

CONCLUSIONS

Our data covered the expression information of 9221 miRTGs. We identified nine miRNA/miRTGs key candidate pairs in response to R. solanacearum infection in ginger. This is the first report of the integrated analysis of miRTGs and miRNAs in response to R. solanacearum infection among ginger species. This study is expected to deliver several insights in understanding the miRNA regulatory network in ginger and mango ginger response to bacterial wilt.

Development, characterization, functional annotation and validation of genomic and genic-SSR markers using de novo next generation sequencing in Melia dubia Cav

Melia dubia Cav. (Meliaceae), a fast-growing tropical tree finds use in plywood, pulp and high-value solid wood products. To increase its productivity, we must essentially capture genetic diversity and identify genotypes with superior wood properties. This study aimed to develop novel microsatellite markers from genomic data and validate the markers in M. dubia. Direct Seq-to-SSR approach was adopted and using an in-house Perl script, 426,390 SSR markers identified. For validation, selected 151 markers, of which 50 were genomic markers chosen randomly, and 101 were genic markers identified through BLAST2GO. Amplification was observed in all loci, and 81.4% generated high-quality, reproducible amplicons of the expected size. Out of 50 genomic markers, we used ten highly polymorphic markers to assess genetic diversity among 75 genotypes from three populations. One hundred fourteen alleles were recorded, with a moderate level of diversity and a positive fixation index. Twenty-nine genic markers representing 13 enzymes showing polymorphism for wood stiffness were selected for diversity assessment of 24 genotypes (12 genotypes each with high and low-stress wave velocity). The product size ranged from 87 to 279, covering the majority of the genome. Cluster and structure analysis segregated ~ 80% of the genotypes based on the trait. This is the first report of the development of genic markers from a genomic survey and has proved efficient in differentiating genotypes based on the trait. The markers developed in this study will be useful for genetic mapping, diversity estimation, marker-assisted selection for desired traits and breeding for wood traits in M. dubia.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02858-w.

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.

De novo transcriptome of the whole-body of the gastropod mollusk Philomycus bilineatus, a pest with medical potential in China

Philomycus bilineatus is a highly common gastropod mollusk pest in China and is also utilized to treat infectious diseases. However, no genomic resources are available for this non-model species. In the present study, the transcriptomic analysis of P. bilineatus was completed. After sequencing using the next generation sequencing technology, 9.11 Gb of clean reads were obtained, which led to the assembly and annotation of 145,523 transcripts and 125,690 unigenes. Unigenes were functionally classified using Gene Ontology (GO), euKaryotic Ortholog Groups of proteins (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG). A total of 27,554 unigenes were assigned into 55 GO terms, 13,989 unigenes were differentiated into 26 KOG categories, and 16,368 unigenes were assigned to 229 KEGG pathways. Furthermore, 16,614 simple sequence repeats (SSRs), 38 olfactory genes, and 40 antimicrobial peptide/protein genes were identified. The transcriptome profile of P. bilineatus will provide a valuable genomic resource for further study, will promote the development of new pest management strategies through interference of chemosensory communication, and will support potential medicinal uses of this species.

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.

Comparative Transcriptome Analysis of Celery Leaf Blades Identified an R2R3-MYB Transcription Factor that Regulates Apigenin Metabolism

Apigenin has been proven to possess many pharmacological properties, but the mechanism of regulation of apigenin biosynthesis in plants remains unclear. Apigenin is the main flavonoid in celery and is mainly accumulated in the middle stage of leaf blade development. In this study, comparative transcriptomic analysis revealed a large number of structural genes and transcription factor genes that may be involved in the apigenin metabolic pathway. On the basis of the apigenin content in different celery accessions, an R2R3-MYB transcription factor gene, named AgMYB1, was isolated from the high apigenin celery accession C014. Bioinformatics analysis indicated that AgMYB1 may be involved in flavonoid metabolism. AgMYB1 expression showed a positive relation with the expression of the apigenin accumulation marker gene FNSI and with the apigenin content in different celery tissues. Moreover, overexpression and antisense expression of AgMYB1 in transgenic celery plants significantly increased and reduced the expression of apigenin biosynthetic genes and the apigenin content, respectively. These findings suggest that AgMYB1 is involved in positive regulation of apigenin metabolism in celery.

The transcriptome analysis of Protaetia brevitarsis Lewis larvae

Larvae of the pest Protaetia brevitarsis are used to treat infections in traditional Chinese medicine. However, genomic information about this non-model species is currently lacking. To better understand the fundamental biology of this non-model species, its transcriptome was obtained using next generation sequencing and then analyzed. A total of 7.62 Gb of clean reads were obtained, which were assembled into 169,087 transcripts corresponding to 142,000 annotated unigenes. These unigenes were functionally classified according to Gene Ontology (GO), euKaryotic Ortholog Groups of proteins (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations. A total of 41,921 unigenes were assigned to 56 GO terms, 21,454 unigenes were divided among 26 KOG categories, and 16,368 unigenes were assigned to 32 KEGG pathways. In addition, 19,144 simple sequence repeats (SSRs) were identified. Furthermore, several kinds of natural antimicrobial peptides and proteins, 4 histones with potential antimicrobial activity, and 41 potential antimicrobial peptide sequences were identified. These data are the first reported whole transcriptome sequence of P. brevitarsis larvae, which represents a valuable genomic resource for studying this species, thus promoting the utilization of its medical potential.

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.

Transcriptome analysis of the Chinese giant salamander (Andrias davidianus) using RNA-sequencing

The Chinese giant salamander (Andrias davidianus) is an economically important animal on academic value. However, the genomic information of this species has been less studied. In our study, the transcripts of A. davidianus were obtained by RNA-seq to conduct a transcriptomic analysis. In total 132,912 unigenes were generated with an average length of 690 bp and N50 of 1263 bp by de novo assembly using Trinity software. Using a sequence similarity search against the nine public databases (CDD, KOG, NR, NT, PFAM, Swiss-prot, TrEMBL, GO and KEGG databases), a total of 24,049, 18,406, 36,711, 15,858, 20,500, 27,515, 36,705, 28,879 and 10,958 unigenes were annotated in databases, respectively. Of these, 6323 unigenes were annotated in all database and 39,672 unigenes were annotated in at least one database. Blasted with KEGG pathway, 10,958 unigenes were annotated, and it was divided into 343 categories according to different pathways. In addition, we also identified 29,790 SSRs. This study provided a valuable resource for understanding transcriptomic information of A. davidianus and laid a foundation for further research on functional gene cloning, genomics, genetic diversity analysis and molecular marker exploitation in A. davidianus.

De novo assembly, characterization and annotation for the transcriptome of Sphaeroma terebrans and microsatellite marker discovery

Sphaeroma terebrans, an economically and ecologically important marine wood-boring isopod, is mainly distributed in tropical and subtropical mangroves. Nevertheless, available genomic and transcriptomic information for this isopod is extremely deficient. Here, we first performed the assembly of S. terebrans transcriptome by Illumina sequencing. A total of 51,092 high-quality unigenes with an average length of 641 bp were obtained and 19,915 unigenes were successfully annotated in four public databases. Only 9932 out of 19915 unigenes were commonly annotated by all four databases. In addition, 9609 unigenes were categorized into 54 function categories of Gene Ontology (GO), and 14,512 unigenes were successfully grouped into 25 functional categories of the EuKaryotic Orthologous Groups (KOG) database. Moreover, 11,507 unigenes were assigned to 228 pathways of the Kyoto Encyclopedia of Genes and Genomes (KEGG). Out of 51,092 unigenes, a total of 4257 different microsatellites with motifs range from di- to hexa-nucleotide were identified from 3324 unigene sequences. Among 64 primer pairs selected for validation, 35 were successful in PCR amplification and 13 exhibited obvious repeat polymorphisms in the wild population of S. terebrans in Dongzhaigang (Hainan Island). The transcriptome dataset and the identified microsatellite markers (SSRs) will provide abundant information for researches on the discovery of new genes, metabolic mechanism and genetic diversity of S. terebrans.

Transcriptome profiling of fully open flowers in a frost-tolerant almond genotype in response to freezing stress

Spring frost is a major limiting abiotic stress for the cultivation of almonds [Prunus dulcis (Mill.)] in Mediterranean areas or the Middle East. Spring frost, in particular, damages almond fully open flowers, resulting to significant reduction in yield. Little is known about the genetic factors expressed after frost stress in Prunus spp. as well as in almond fully open flowers. Here, we provide the molecular signature of pistils of fully open flowers from a frost-tolerant almond genotype. The level of frost tolerance in this genotype was determined for all three flowering stages and was confirmed by comparing it to two other cultivars using several physiological analyses. Afterwards, comprehensive expression profiling of genes expressed in fully open flowers was performed after being exposed to frost temperatures (during post-thaw period). Clean reads, 27,104,070 and 32,730,772, were obtained for non-frost-treated and frost-treated (FT) libraries, respectively. A total of 62.24 Mb was assembled, generating 50,896 unigenes and 66,906 transcripts. Therefore, 863 upregulated genes and 555 downregulated genes were identified in the FT library. Functional annotation showed that most of the upregulated genes were related to various biological processes involved in responding to abiotic stress. For the first time, a highly expressed cold-shock protein was identified in the reproductive organ of fruit trees. The expression of six genes was validated by RT-PCR. As the first comprehensive analysis of open flowers in a frost-tolerant almond genotype, this study represents a key step toward the molecular breeding of fruit tree species for frost tolerance.

Development of hop transcriptome to support research into host-viroid interactions

Viroids, the smallest known pathogens, unable to encode any proteins, can cause severe diseases in their host plants. One of the proposed mechanisms of their pathogenicity includes silencing the host's genes via viroid-derived small RNAs, which are products of the host's immune response to the viroid's double stranded RNA. Humulus lupulus (hop) plants are hosts to several viroids; two of them, HLVd and CBCVd, are interesting models for studying host-viroid interactions, due to the symptomless infection of the former and severe stunting disease caused by the latter. To study these interactions, we constructed a deep hop NGS transcriptome based on 35 Gb paired-end sequencing data assembled into over 74 Mb of contigs. These transcripts were used for in-silico prediction of target transcripts of vd-sRNA of the two aforementioned viroids, using two different software tools. Prediction models revealed that 1062 and 1387 hop transcripts share nucleotide similarities with HLVd- and CBCVd-derived small RNAs, respectively, so they could be silenced in an RNA interference process. Furthermore, we selected 17 transcripts from 4 groups of targets involved in the metabolism of plant hormones, small RNA biogenesis, transcripts with high complementarity with viroid-derived small RNAs and transcripts targeted by CBCVd-derived small RNAs with high cellular concentrations. Their expression was monitored by reverse transcription quantitative PCR performed using leaf, flower and cone samples. Additionally, the expression of 5 pathogenesis related genes was monitored. Expression analysis confirmed high expression levels of four pathogenesis related genes in leaves of HLVd and CBCVd infected hop plants. Expression fluctuations were observed for the majority of targets, with possible evidence of downregulation of GATA transcription factor by CBCVd- and of linoleate 13S-lipoxygenase by HLVd-derived small RNAs. These results provide a deep transcriptome of hop and the first insights into complex viroid-hop plant interactions.

Characterization of the transcriptome and EST-SSR development in Boea clarkeana, a desiccation-tolerant plant endemic to China

BACKGROUND

Desiccation-tolerant (DT) plants can recover full metabolic competence upon rehydration after losing most of their cellular water (>95%) for extended periods of time. Functional genomic approaches such as transcriptome sequencing can help us understand how DT plants survive and respond to dehydration, which has great significance for plant biology and improving the drought tolerance of crops. Boea clarkeana Hemsl. (Gesneriaceae) is a DT dicotyledonous herb. Its genomic sequences characteristics remain unknown. Based on transcriptomic analyses, polymorphic EST-SSR (simple sequence repeats in expressed sequence tags) molecular primers can be designed, which will greatly facilitate further investigations of the population genetics and demographic histories of DT plants.

METHODS

In the present study, we used the platform Illumina HiSeq™2000 and de novo assembly technology to obtain leaf transcriptomes of B. clarkeana and conducted a BLASTX alignment of the sequencing data and protein databases for sequence classification and annotation. Then, based on the sequence information, the EST-SSR markers were developed, and the functional annotation of ESTs containing polymorphic SSRs were obtained through BLASTX.

RESULTS

A total of 91,449 unigenes were generated from the leaf cDNA library of B. clarkeana. Based on a sequence similarity search with a known protein database, 72,087 unigenes were annotated. Among the annotated unigenes, a total of 71,170 unigenes showed significant similarity to the known proteins of 463 popular model species in the Nr database, and 59,962 unigenes and 32,336 unigenes were assigned to Gene Ontology (GO) classifications and Cluster of Orthologous Groups (COG), respectively. In addition, 44,924 unigenes were mapped in 128 KEGG pathways. Furthermore, a total of 7,610 unigenes with 8,563 microsatellites were found. Seventy-four primer pairs were selected from 436 primer pairs designed for polymorphism validation. SSRs with higher polymorphism rates were concentrated on dinucleotides, pentanucleotides and hexanucleotides. Finally, 17 pairs with stable, highly polymorphic loci were selected for polymorphism screening. There was a total of 65 alleles, with 2-6 alleles at each locus. Primarily due to the unique biological characteristics of plants, the HE (0-0.196), HO (0.082-0.14) and PIC (0-0.155) per locus were very low. The functional annotation distribution centered on ESTs containing di- and tri-nucleotide SSRs, and the ESTs containing primers BC2, BC4 and BC12 were annotated to vegetative dehydration/desiccation pathways.

DISCUSSION

This work is the first genetic study of B. clarkeana as a new plant resource of DT genes. A substantial number of transcriptome sequences were generated in this study. These sequences are valuable resources for gene annotation and discovery as well as molecular marker development. These sequences could also provide a valuable basis for future molecular studies of B. clarkeana.

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.

Genetic Transformation and Genomic Resources for Next-Generation Precise Genome Engineering in Vegetable Crops

In the frame of modern agriculture facing the predicted increase of population and general environmental changes, the securement of high quality food remains a major challenge to deal with. Vegetable crops include a large number of species, characterized by multiple geographical origins, large genetic variability and diverse reproductive features. Due to their nutritional value, they have an important place in human diet. In recent years, many crop genomes have been sequenced permitting the identification of genes and superior alleles associated with desirable traits. Furthermore, innovative biotechnological approaches allow to take a step forward towards the development of new improved cultivars harboring precise genome modifications. Sequence-based knowledge coupled with advanced biotechnologies is supporting the widespread application of new plant breeding techniques to enhance the success in modification and transfer of useful alleles into target varieties. Clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 system, zinc-finger nucleases, and transcription activator-like effector nucleases represent the main methods available for plant genome engineering through targeted modifications. Such technologies, however, require efficient transformation protocols as well as extensive genomic resources and accurate knowledge before they can be efficiently exploited in practical breeding programs. In this review, we revise the state of the art in relation to availability of such scientific and technological resources in various groups of vegetables, describe genome editing results obtained so far and discuss the implications for future applications.

De novo assembly, functional annotation, and analysis of the giant reed (Arundo donax L.) leaf transcriptome provide tools for the development of a biofuel feedstock

BACKGROUND

Arundo donax has attracted renewed interest as a potential candidate energy crop for use in biomass-to-liquid fuel conversion processes and biorefineries. This is due to its high productivity, adaptability to marginal land conditions, and suitability for biofuel and biomaterial production. Despite its importance, the genomic resources currently available for supporting the improvement of this species are still limited.

RESULTS

We used RNA sequencing (RNA-Seq) to de novo assemble and characterize the A. donax leaf transcriptome. The sequencing generated 1249 million clean reads that were assembled using single-k-mer and multi-k-mer approaches into 62,596 unique sequences (unitranscripts) with an N50 of 1134 bp. TransDecoder and Trinotate software suites were used to obtain putative coding sequences and annotate them by mapping to UniProtKB/Swiss-Prot and UniRef90 databases, searching for known transcripts, proteins, protein domains, and signal peptides. Furthermore, the unitranscripts were annotated by mapping them to the NCBI non-redundant, GO and KEGG pathway databases using Blast2GO. The transcriptome was also characterized by BLAST searches to investigate homologous transcripts of key genes involved in important metabolic pathways, such as lignin, cellulose, purine, and thiamine biosynthesis and carbon fixation. Moreover, a set of homologous transcripts of key genes involved in stomatal development and of genes coding for stress-associated proteins (SAPs) were identified. Additionally, 8364 simple sequence repeat (SSR) markers were identified and surveyed. SSRs appeared more abundant in non-coding regions (63.18%) than in coding regions (36.82%). This SSR dataset represents the first marker catalogue of A. donax. 53 SSRs (PolySSRs) were then predicted to be polymorphic between ecotype-specific assemblies, suggesting genetic variability in the studied ecotypes.

CONCLUSIONS

This study provides the first publicly available leaf transcriptome for the A. donax bioenergy crop. The functional annotation and characterization of the transcriptome will be highly useful for providing insight into the molecular mechanisms underlying its extreme adaptability. The identification of homologous transcripts involved in key metabolic pathways offers a platform for directing future efforts in genetic improvement of this species. Finally, the identified SSRs will facilitate the harnessing of untapped genetic diversity. This transcriptome should be of value to ongoing functional genomics and genetic studies in this crop of paramount economic importance.

Transcriptome Analyses Reveal Candidate Pod Shattering-Associated Genes Involved in the Pod Ventral Sutures of Common Vetch (Vicia sativa L.)

The seed dispersion caused by pod shattering is a form of propagation used by many wild species. Loss of seeds from pod shattering is frequent in the common vetch (Vicia sativa L.), an important self-pollinating annual forage legume. However, pod shattering is one of the most important defects that limits the reproduction of the vetch in the field and the usage as a leguminous forage crop. To better understand the vetch pod shattering mechanism, we used high-throughput RNA sequencing to assess the global changes in the transcriptomes of the pod ventral sutures of shattering-susceptible and shattering-resistant vetch accessions screened from 541 vetch germplasms. A total of 1,285 significantly differentially expressed unigenes (DEGs) were detected, including 575 up-regulated unigenes and 710 down-regulated unigenes. Analyses of Gene Ontology and KEGG metabolic enrichment pathways of 1,285 DEGs indicated that 22 DEGs encoding cell wall modifications and hydrolases associated with pod shattering were highly expressed in shattering-susceptible accessions. These genes were mainly enriched in "hydrolase activity," "cytoplasm," and "carbohydrate metabolic process" systems. These cell wall modifications and hydrolases genes included β-glucosidase and endo-polygalacturonase, which work together to break down the glycosidic bonds of pectin and cellulose, and to promote the dissolution and disappearance of the cell wall in the ventral suture of the pod and make the pod more susceptible to shattering. We demonstrated the differences in gene transcription levels between the shattering-susceptible and shattering-resistant vetch accessions for the first time and our results provided valuable information for the identifying and characterizing of pod shattering regulation networks in vetch. This information may facilitate the future identification of pod shattering-related genes and their underlying molecular mechanisms in the common vetch.

Genome-Wide Transcriptional Excavation of Dipsacus asperoides Unmasked both Cryptic Asperosaponin Biosynthetic Genes and SSR Markers

BACKGROUND

Dipsacus asperoides is a traditional Chinese medicinal crop. The root is generally used as a medicine and is frequently prescribed by Chinese doctors for the treatment of back pain, limb paralysis, flutter trauma, tendon injuries, and fractures. With the rapid development of bioinformatics, research has been focused on this species at the gene or molecular level. For purpose of fleshing out genome information about D. asperoides, in this paper we conducted transcriptome analysis of this species.

PRINCIPAL FINDINGS

To date, many genes encoding enzymes involved in the biosynthesis of triterpenoid saponins in D.asperoides have not been elucidated. Illumina paired-end sequencing was employed to probe D. asperoides's various enzymes associated with the relevant mesostate. A total of 30, 832,805 clean reads and de novo spliced 43,243 unigenes were obtained. Of all unigenes, only 8.27% (3578) were successfully annotated in total of seven public databases: Nr, Nt, Swiss-Prot, GO, KOG, KEGG, and Pfam, which might be attributed to the poor studies on D. asperoides. The candidate genes encoding enzymes involved in triterpenoid saponin biosynthesis were identified and experimentally verified by reverse transcription qPCR, encompassing nine cytochrome P450s and 17 UDP-glucosyltransferases. Specifically, unearthly putative genes involved in the glycosylation of hederagenin were acquired. Simultaneously, 4490 SSRs from 43,243 examined sequences were determined via bioinformatics analysis.

CONCLUSION

This study represents the first report on the use of the Illumina sequence platform on this crop at the transcriptome level. Our findings of candidate genes encoding enzymes involved in Dipsacus saponin VI biosynthes is provide novel information in efforts to further understand the triterpenoid metabolic pathway on this species. The initial genetics resources in this study will contribute significantly to the genetic breeding program of D. asperoides, and are beneficial for clinical diagnosis and treatment.

Transcriptome Sequencing and Development of Genic SSR Markers of an Endangered Chinese Endemic Genus Dipteronia Oliver (Aceraceae)

Dipteronia Oliver (Aceraceae) is an endangered Chinese endemic genus consisting of two living species, Dipteronia sinensis and Dipteronia dyeriana. However, studies on the population genetics and evolutionary analyses of Dipteronia have been hindered by limited genomic resources and genetic markers. Here, the generation, de novo assembly and annotation of transcriptome datasets, and a large set of microsatellite or simple sequence repeat (SSR) markers derived from Dipteronia have been described. After Illumina pair-end sequencing, approximately 93.2 million reads were generated and assembled to yield a total of 99,358 unigenes. A majority of these unigenes (53%, 52,789) had at least one blast hit against the public protein databases. Further, 12,377 SSR loci were detected and 4179 primer pairs were designed for experimental validation. Of these 4179 primer pairs, 435 primer pairs were randomly selected to test polymorphism. Our results show that products from 132 primer pairs were polymorphic, in which 97 polymorphic SSR markers were further selected to analyze the genetic diversity of 10 natural populations of Dipteronia. The identification of SSR markers during our research will provide the much valuable data for population genetic analyses and evolutionary studies in Dipteronia.

De novo Transcriptome Assembly of a Chinese Locoweed (Oxytropis ochrocephala) Species Provides Insights into Genes Associated with Drought, Salinity, and Cold Tolerance

BACKGROUND

Locoweeds (toxic Oxytropis and Astraglus species), containing the toxic agent swainsonine, pose serious threats to animal husbandry on grasslands in both China and the US. Some locoweeds have evolved adaptations in order to resist various stress conditions such as drought, salt and cold. As a result they replace other plants in their communities and become an ecological problem. Currently very limited genetic information of locoweeds is available and this hinders our understanding in the molecular basis of their environmental plasticity, and the interaction between locoweeds and their symbiotic swainsonine producing endophytes. Next-generation sequencing provides a means of obtaining transcriptomic sequences in a timely manner, which is particularly useful for non-model plants. In this study, we performed transcriptome sequencing of Oxytropis ochrocephala plants followed by a de nove assembly. Our primary aim was to provide an enriched pool of genetic sequences of an Oxytropis sp. for further locoweed research.

RESULTS

Transcriptomes of four different O. ochrocephala samples, from control (CK) plants, and those that had experienced either drought (20% PEG), salt (150 mM NaCl) or cold (4°C) stress were sequenced using an Illumina Hiseq 2000 platform. From 232,209,506 clean reads 23,220,950,600 (~23 G nucleotides), 182,430 transcripts and 88,942 unigenes were retrieved, with an N50 value of 1237. Differential expression analysis revealed putative genes encoding heat shock proteins (HSPs) and late embryogenesis abundant (LEA) proteins, enzymes in secondary metabolite and plant hormone biosyntheses, and transcription factors which are involved in stress tolerance in O. ochrocephala. In order to validate our sequencing results, we further analyzed the expression profiles of nine genes by quantitative real-time PCR. Finally, we discuss the possible mechanism of O. ochrocephala's adaptations to stress environment.

CONCLUSION

Our transcriptome sequencing data present useful genetic information of a locoweed species. This genetic information will underpin further research in elucidating the environmental acclimation mechanism in locoweeds and the endophyte-plant association.

De novo assembly and functional annotation of Myrciaria dubia fruit transcriptome reveals multiple metabolic pathways for L-ascorbic acid biosynthesis

Background

Myrciaria dubia is an Amazonian fruit shrub that produces numerous bioactive phytochemicals, but is best known by its high L-ascorbic acid (AsA) content in fruits. Pronounced variation in AsA content has been observed both within and among individuals, but the genetic factors responsible for this variation are largely unknown. The goals of this research, therefore, were to assemble, characterize, and annotate the fruit transcriptome of M. dubia in order to reconstruct metabolic pathways and determine if multiple pathways contribute to AsA biosynthesis.

Results

In total 24,551,882 high-quality sequence reads were de novo assembled into 70,048 unigenes (mean length = 1150 bp, N50 = 1775 bp). Assembled sequences were annotated using BLASTX against public databases such as TAIR, GR-protein, FB, MGI, RGD, ZFIN, SGN, WB, TIGR_CMR, and JCVI-CMR with 75.2 % of unigenes having annotations. Of the three core GO annotation categories, biological processes comprised 53.6 % of the total assigned annotations, whereas cellular components and molecular functions comprised 23.3 and 23.1 %, respectively. Based on the KEGG pathway assignment of the functionally annotated transcripts, five metabolic pathways for AsA biosynthesis were identified: animal-like pathway, myo-inositol pathway, L-gulose pathway, D-mannose/L-galactose pathway, and uronic acid pathway. All transcripts coding enzymes involved in the ascorbate-glutathione cycle were also identified. Finally, we used the assembly to identified 6314 genic microsatellites and 23,481 high quality SNPs.

Conclusions

This study describes the first next-generation sequencing effort and transcriptome annotation of a non-model Amazonian plant that is relevant for AsA production and other bioactive phytochemicals. Genes encoding key enzymes were successfully identified and metabolic pathways involved in biosynthesis of AsA, anthocyanins, and other metabolic pathways have been reconstructed. The identification of these genes and pathways is in agreement with the empirically observed capability of M. dubia to synthesize and accumulate AsA and other important molecules, and adds to our current knowledge of the molecular biology and biochemistry of their production in plants. By providing insights into the mechanisms underpinning these metabolic processes, these results can be used to direct efforts to genetically manipulate this organism in order to enhance the production of these bioactive phytochemicals.

The accumulation of AsA precursor and discovery of genes associated with their biosynthesis and metabolism in M. dubia is intriguing and worthy of further investigation. The sequences and pathways produced here present the genetic framework required for further studies. Quantitative transcriptomics in concert with studies of the genome, proteome, and metabolome under conditions that stimulate production and accumulation of AsA and their precursors are needed to provide a more comprehensive view of how these pathways for AsA metabolism are regulated and linked in this species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2225-6) contains supplementary material, which is available to authorized users.

De novo sequencing and assembly analysis of transcriptome in the Sodom apple (Calotropis gigantea)

Background

The Sodom apple (Calotropis gigantea), a member of the Asclepiadaceae family, is a large evergreen shrub native to continental Asia and northern Africa. As an important medicinal shrub and a fiber resource plant, there is an urgent need for developing molecular markers to facilitate breeding and genetic improvement of varieties.

Results

In this study, using the Illumina high throughput sequencing technique we obtained about 45 million paired end sequencing reads, De novo assembled and generated a total of 133,634 transcripts with a mean of 1837.47 bp in length. Based on protein homology searches against available databases, a total of 21,851 unigenes were functionally annotated. In particular, many transcripts that encode for putative proteins involved in fiber and secondary metabolite biosynthesis were identified and analyzed. Key fiber genes identified were validated experimentally through Real-Time PCR technique. Various transcription factors involved in regulating plant response to abiotic stress were also identified. In addition, based on the unigene sequences assembled, 11,623 microsatellites loci were detected, which provide very useful resources for developing microsatellite molecular markers.

Conclusion

This study is the first report on transcriptome information in the Calotropis species and provides rich gene transcript resources for conducting further studies on understanding the molecular basis of fiber and secondary metabolite biosynthesis, serving the genetic improvement and resource utilization in Calotropis plants.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1908-3) contains supplementary material, which is available to authorized users.

RNA-seq Reveals Complicated Transcriptomic Responses to Drought Stress in a Nonmodel Tropic Plant, Bombax ceiba L

High-throughput transcriptome provides an unbiased approach for understanding the genetic basis and gene functions in response to different conditions. Here we sequenced RNA-seq libraries derived from a Bombax ceiba L. system under a controlled experiment. As a known medicinal and ornamental plant, B. ceiba grows mainly in hot-dry monsoon rainforests in Southeast Asia and Australia. Due to the specific growth environment, it has evolved a unique system that enables a physiologic response to drought stress. To date, few studies have characterized the genome-wide features of drought endurance in B. ceiba. In this study, we first attempted to characterize and identify the most differentially expressed genes and associated functional pathways under drought treatment and normal condition. Using RNA-seq technology, we generated the first transcriptome of B. ceiba and identified 59 differentially expressed genes with greater than 1,000-fold changes under two conditions. The set of upregulated genes implicates interplay among various pathways: plants growth, ubiquitin-mediated proteolysis, polysaccharides hydrolyzation, oxidative phosphorylation and photosynthesis, etc. In contrast, genes associated with stem growth, cell division, fruit ripening senescence, disease resistance, and proline synthesis are repressed. Notably, key genes of high RPKM levels in drought are AUX1, JAZ, and psbS, which are known to regulate the growth of plants, the resistance against abiotic stress, and the photosynthesis process. Furthermore, 16,656 microsatellite markers and 3,071 single-nucleotide polymorphisms (SNPs) were predicted by in silico methods. The identification and functional annotation of differentially expressed genes, microsatellites, and SNPs represent a major step forward and would serve as a valuable resource for understanding the complexity underlying drought endurance and adaptation in B. ceiba.

De Novo Characterization of Flower Bud Transcriptomes and the Development of EST-SSR Markers for the Endangered Tree Tapiscia sinensis

Tapiscia sinensis Oliv (Tapisciaceae) is an endangered species native to China famous for its androdioecious breeding system. However, there is a lack of genomic and transcriptome data on this species. In this study, the Tapiscia sinensis transcriptomes from two types of sex flower buds were sequenced. A total of 97,431,176 clean reads were assembled into 52,169 unigenes with an average length of 1116 bp. Through similarity comparison with known protein databases, 36,662 unigenes (70.27%) were annotated. A total of 10,002 (19.17%) unigenes were assigned to 124 pathways using the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database. Additionally, 10,371 simple sequence repeats (SSRs) were identified in 8608 unigenes, with 16,317 pairs of primers designed for applications. 150 pairs of primers were chosen for further validation, and the 68 pairs (45.5%) were able to produce clear polymorphic bands. Six polymorphic SSR markers were used to Bayesian clustering analysis of 51 T. sinensis individuals. This is the first report to provide transcriptome information and to develop large-scale SSR molecular markers for T. sinensis. This study provides a valuable resource for conservation genetics and functional genomics research on T. sinensis for future work.

De novo Assembly, Characterization of Immature Seed Transcriptome and Development of Genic-SSR Markers in Black Gram [Vigna mungo (L.) Hepper]

Black gram [V. mungo (L.) Hepper] is an important legume crop extensively grown in south and south-east Asia, where it is a major source of dietary protein for its predominantly vegetarian population. However, lack of genomic information and markers has become a limitation for genetic improvement of this crop. Here, we report the transcriptome sequencing of the immature seeds of black gram cv. TU94-2, by Illumina paired end sequencing technology to generate transcriptome sequences for gene discovery and genic-SSR marker development. A total of 17.2 million paired-end reads were generated and 48,291 transcript contigs (TCS) were assembled with an average length of 443 bp. Based on sequence similarity search, 33,766 TCS showed significant similarity to known proteins. Among these, only 29,564 TCS were annotated with gene ontology (GO) functional categories. A total number of 138 unique KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways were identified, of which majority of TCS are grouped into purine metabolism (678) followed by pyrimidine metabolism (263). A total of 48,291 TCS were searched for SSRs and 1,840 SSRs were identified in 1,572 TCS with an average frequency of one SSR per 11.9 kb. The tri-nucleotide repeats were most abundant (35%) followed by di-nucleotide repeats (32%). PCR primer pairs were successfully designed for 933 SSR loci. Sequences analyses indicate that about 64.4% and 35.6% of the SSR motifs were present in the coding sequences (CDS) and untranslated regions (UTRs) respectively. Tri-nucleotide repeats (57.3%) were preferentially present in the CDS. The rate of successful amplification and polymorphism were investigated using selected primers among 18 black gram accessions. Genic-SSR markers developed from the Illumina paired end sequencing of black gram immature seed transcriptome will provide a valuable resource for genetic diversity, evolution, linkage mapping, comparative genomics and marker-assisted selection in black gram.

De Novo Assembly and Transcriptome Analysis of Wheat with Male Sterility Induced by the Chemical Hybridizing Agent SQ-1

Wheat (Triticum aestivum L.), one of the world's most important food crops, is a strictly autogamous (self-pollinating) species with exclusively perfect flowers. Male sterility induced by chemical hybridizing agents has increasingly attracted attention as a tool for hybrid seed production in wheat; however, the molecular mechanisms of male sterility induced by the agent SQ-1 remain poorly understood due to limited whole transcriptome data. Therefore, a comparative analysis of wheat anther transcriptomes for male fertile wheat and SQ-1-induced male sterile wheat was carried out using next-generation sequencing technology. In all, 42,634,123 sequence reads were generated and were assembled into 82,356 high-quality unigenes with an average length of 724 bp. Of these, 1,088 unigenes were significantly differentially expressed in the fertile and sterile wheat anthers, including 643 up-regulated unigenes and 445 down-regulated unigenes. The differentially expressed unigenes with functional annotations were mapped onto 60 pathways using the Kyoto Encyclopedia of Genes and Genomes database. They were mainly involved in coding for the components of ribosomes, photosynthesis, respiration, purine and pyrimidine metabolism, amino acid metabolism, glutathione metabolism, RNA transport and signal transduction, reactive oxygen species metabolism, mRNA surveillance pathways, protein processing in the endoplasmic reticulum, protein export, and ubiquitin-mediated proteolysis. This study is the first to provide a systematic overview comparing wheat anther transcriptomes of male fertile wheat with those of SQ-1-induced male sterile wheat and is a valuable source of data for future research in SQ-1-induced wheat male sterility.

De novo transcriptome assembly of Ipomoea nil using Illumina sequencing for gene discovery and SSR marker identification

Ipomoea nil is widely used as an ornamental plant due to its abundance of flower color, but the limited transcriptome and genomic data hinder research on it. Using illumina platform, transcriptome profiling of I. nil was performed through high-throughput sequencing, which was proven to be a rapid and cost-effective means to characterize gene content. Our goal is to use the resulting information to facilitate the relevant research on flowering and flower color formation in I. nil. In total, 268 million unique illumina RNA-Seq reads were produced and used in the transcriptome assembly. These reads were assembled into 220,117 contigs, of which 137,307 contigs were annotated using the GO and KEGG database. Based on the result of functional annotations, a total of 89,781 contigs were assigned 455,335 GO term annotations. Meanwhile, 17,418 contigs were identified with pathway annotation and they were functionally assigned to 144 KEGG pathways. Our transcriptome revealed at least 55 contigs as probably flowering-related genes in I. nil, and we also identified 25 contigs that encode key enzymes in the phenylpropanoid biosynthesis pathway. Based on the analysis relating to gene expression profiles, in the phenylpropanoid biosynthesis pathway of I. nil, the repression of lignin biosynthesis might lead to the redirection of the metabolic flux into anthocyanin biosynthesis. This may be the most likely reason that I. nil has high anthocyanins content, especially in its flowers. Additionally, 15,537 simple sequence repeats (SSRs) were detected using the MISA software, and these SSRs will undoubtedly benefit future breeding work. Moreover, the information uncovered in this study will also serve as a valuable resource for understanding the flowering and flower color formation mechanisms in I. nil.

RNA-Seq analysis and gene discovery of Andrias davidianus using Illumina short read sequencing

The Chinese giant salamander, Andrias davidianus, is an important species in the course of evolution; however, there is insufficient genomic data in public databases for understanding its immunologic mechanisms. High-throughput transcriptome sequencing is necessary to generate an enormous number of transcript sequences from A. davidianus for gene discovery. In this study, we generated more than 40 million reads from samples of spleen and skin tissue using the Illumina paired-end sequencing technology. De novo assembly yielded 87,297 transcripts with a mean length of 734 base pairs (bp). Based on the sequence similarities, searching with known proteins, 38,916 genes were identified. Gene enrichment analysis determined that 981 transcripts were assigned to the immune system. Tissue-specific expression analysis indicated that 443 of transcripts were specifically expressed in the spleen and skin. Among these transcripts, 147 transcripts were found to be involved in immune responses and inflammatory reactions, such as fucolectin, β-defensins and lymphotoxin beta. Eight tissue-specific genes were selected for validation using real time reverse transcription quantitative PCR (qRT-PCR). The results showed that these genes were significantly more expressed in spleen and skin than in other tissues, suggesting that these genes have vital roles in the immune response. This work provides a comprehensive genomic sequence resource for A. davidianus and lays the foundation for future research on the immunologic and disease resistance mechanisms of A. davidianus and other amphibians.

De novo transcriptome assembly of pummelo and molecular marker development

Pummelo (Citrus grandis) is an important fruit crop worldwide because of its nutritional value. To accelerate the pummelo breeding program, it is essential to obtain extensive genetic information and develop relative molecular markers. Here, we obtained a 12-Gb transcriptome dataset of pummelo through a mixture of RNA from seven tissues using Illumina pair-end sequencing, assembled into 57,212 unigenes with an average length of 1010 bp. The annotation and classification results showed that a total of 39,584 unigenes had similar hits to the known proteins of four public databases, and 31,501 were classified into 55 Gene Ontology (GO) functional sub-categories. The search for putative molecular markers among 57,212 unigenes identified 10,276 simple sequence repeats (SSRs) and 64,720 single nucleotide polymorphisms (SNPs). High-quality primers of 1174 SSR loci were designed, of which 88.16% were localized to nine chromosomes of sweet orange. Of 100 SSR primers that were randomly selected for testing, 87 successfully amplified clear banding patterns. Of these primers, 29 with a mean PIC (polymorphic information content) value of 0.52 were effectively applied for phylogenetic analysis. Of the 20 SNP primers, 14 primers, including 54 potential SNPs, yielded target amplifications, and 46 loci were verified via Sanger sequencing. This new dataset will be a valuable resource for molecular biology studies of pummelo and provides reliable information regarding SNP and SSR marker development, thus expediting the breeding program of pummelo.

De novo assembly, transcriptome characterization, lignin accumulation, and anatomic characteristics: novel insights into lignin biosynthesis during celery leaf development

Celery of the family Apiaceae is a biennial herb that is cultivated and consumed worldwide. Lignin is essential for cell wall structural integrity, stem strength, water transport, mechanical support, and plant pathogen defense. This study discussed the mechanism of lignin formation at different stages of celery development. The transcriptome profile, lignin distribution, anatomical characteristics, and expression profile of leaves at three stages were analyzed. Regulating lignin synthesis in celery growth development has a significant economic value. Celery leaves at three stages were collected, and Illumina paired-end sequencing technology was used to analyze large-scale transcriptome sequences. From Stage 1 to 3, the collenchyma and vascular bundles in the petioles and leaf blades thickened and expanded, whereas the phloem and the xylem extensively developed. Spongy and palisade mesophyll tissues further developed and were tightly arranged. Lignin accumulation increased in the petioles and the mesophyll (palisade and spongy), and the xylem showed strong lignification. Lignin accumulation in different tissues and at different stages of celery development coincides with the anatomic characteristics and transcript levels of genes involved in lignin biosynthesis. Identifying the genes that encode lignin biosynthesis-related enzymes accompanied by lignin distribution may help elucidate the regulatory mechanisms of lignin biosynthesis in celery.

De novo transcriptome assembly, gene annotation, marker development, and miRNA potential target genes validation under abiotic stresses in Oenanthe javanica

Oenanthe javanica is an aquatic perennial herb with known medicinal properties and an edible vegetable with high vitamin and mineral content. The understanding of the biology of O. javanica is limited by the absence of information on its genome, transcriptome, and small RNA. In this study, transcriptome sequencing and small RNA sequencing were performed to annotate function genes, develop SSR markers and analyze potential target genes of miRNAs in O. javanica. All reads with total nucleotides number of 1,440,321,408 bp were assembled into 58,072 transcripts and 40,208 unigenes. A total of 1,233 SSRs were identified from O. javanica. Generated unigenes were aligned against seven databases and annotated with functions. A total of 29 potential targets were predicted. Expression of 10 miRNAs and their corresponding target genes under abiotic stresses (heat, cold, salinity, and drought) was validated. All ten miRNAs were confirmed to response to abiotic stresses. A pair of miRNA and its target gene was found. This study can serve as a valuable resource for future studies on O. javanica, which may focus on novel gene discovery, SSR development, gene mapping, and miRNA-affected processes and pathways. This can promote the development of the useful medicinal properties of O. javanica in medical science.

De novo transcriptome sequencing of the snail Echinolittorina malaccana: identification of genes responsive to thermal stress and development of genetic markers for population studies

Echinolittorina snails inhabit the upper intertidal rocky shore and face strong selection pressures from thermal extremes and fluctuations. Revealing the molecular processes of adaptive significance is greatly obstructed by the scarcity of genomic resource for these taxa. Here, we reported the first comprehensive transcriptome dataset for the genus Echinolittorina. Using Illumina HiSeq 2000 platform, about 52 M and 54 M paired-end clean reads were, respectively, generated for the control and heat-stressed libraries. Totally, 115,211 unique transcript fragments (unigenes) were assembled, with an average length of 453 bp and a N50 size of 492 bp. Approximately one third of the unigenes could be annotated according to their homology matches against the Nr, Swiss-Prot, COG, or KEGG databases, and they were found to represent 23,098 non-redundant genes. Gene expression comparison revealed that 1,267 and 6,663 annotated genes were, respectively, up- and downregulated with at least twofold changes upon heat stress. Gene Ontology and KEGG pathway analyses indicated that there were overrepresented amount of genes enriched in a broad spectrum of biological processes and pathways, including those associated with cytoskeleton organization, developmental regulation, signaling transduction, infection, and cardiac function. In addition, a transcriptome-wide search for polymorphic loci yielded a total of 11,228 simple sequence repeats (SSRs) from 9,938 unigenes and 138,631 single nucleotide polymorphism (SNP) and insertion/deletion (INDEL) sites among 22,770 unigenes. The large number of transcript sequences acquired, the biological pathways identified, and the candidate microsatellite and SNP/INDEL loci discovered in the study will serve as valuable resources for further investigations of genetic differentiation and thermal adaptation among populations.

De novo transcriptome sequence assembly and identification of AP2/ERF transcription factor related to abiotic stress in parsley (Petroselinum crispum)

Parsley is an important biennial Apiaceae species that is widely cultivated as herb, spice, and vegetable. Previous studies on parsley principally focused on its physiological and biochemical properties, including phenolic compound and volatile oil contents. However, little is known about the molecular and genetic properties of parsley. In this study, 23,686,707 high-quality reads were obtained and assembled into 81,852 transcripts and 50,161 unigenes for the first time. Functional annotation showed that 30,516 unigenes had sequence similarity to known genes. In addition, 3,244 putative simple sequence repeats were detected in curly parsley. Finally, 1,569 of the identified unigenes belonged to 58 transcription factor families. Various abiotic stresses have a strong detrimental effect on the yield and quality of parsley. AP2/ERF transcription factors have important functions in plant development, hormonal regulation, and abiotic response. A total of 88 putative AP2/ERF factors were identified from the transcriptome sequence of parsley. Seven AP2/ERF transcription factors were selected in this study to analyze the expression profiles of parsley under different abiotic stresses. Our data provide a potentially valuable resource that can be used for intensive parsley research.

De novo assembly of the perennial ryegrass transcriptome using an RNA-Seq strategy

BACKGROUND

Perennial ryegrass is a highly heterozygous outbreeding grass species used for turf and forage production. Heterozygosity can affect de-Bruijn graph assembly making de novo transcriptome assembly of species such as perennial ryegrass challenging. Creating a reference transcriptome from a homozygous perennial ryegrass genotype can circumvent the challenge of heterozygosity. The goals of this study were to perform RNA-sequencing on multiple tissues from a highly inbred genotype to develop a reference transcriptome. This was complemented with RNA-sequencing of a highly heterozygous genotype for SNP calling.

RESULT

De novo transcriptome assembly of the inbred genotype created 185,833 transcripts with an average length of 830 base pairs. Within the inbred reference transcriptome 78,560 predicted open reading frames were found of which 24,434 were predicted as complete. Functional annotation found 50,890 transcripts with a BLASTp hit from the Swiss-Prot non-redundant database, 58,941 transcripts with a Pfam protein domain and 1,151 transcripts encoding putative secreted peptides. To evaluate the reference transcriptome we targeted the high-affinity K+ transporter gene family and found multiple orthologs. Using the longest unique open reading frames as the reference sequence, 64,242 single nucleotide polymorphisms were found. One thousand sixty one open reading frames from the inbred genotype contained heterozygous sites, confirming the high degree of homozygosity.

CONCLUSION

Our study has developed an annotated, comprehensive transcriptome reference for perennial ryegrass that can aid in determining genetic variation, expression analysis, genome annotation, and gene mapping.

De novo assembly, gene annotation, and simple sequence repeat marker development using Illumina paired-end transcriptome sequences in the pearl oyster Pinctada maxima

We analyzed the mantle transcriptome of pearl oyster Pinctada maxima and developed EST-SSR markers using Illumina HiSeq 2000 paired-end sequencing technology. A total of 49,500,748 raw reads were generated. De novo assembly generated 108,704 unigenes with an average length of 407 bp. Sequence similarity search with known proteins or nucleotides revealed that 30,200 (27.78%) and 25,824 (23.76%) consensus sequences were homologous with the sequences in the non-redundant protein and Swiss-Prot databases, respectively, and that 19,701 (18.12%) of these unigenes were possibly involved in approximately 234 known signaling pathways in the Kyoto Encyclopedia of Genes and Genomes database. Ninety one biomineralization-related unigenes were detected. In a cultured stock, 1764 simple sequence repeats were identified and 56 primer pairs were randomly selected and tested. The rate of successful amplification was 68.3%. The developed molecular markers are helpful for further studies on genetic linkage analysis, gene localization, and quantitative trait loci mapping.